Next Article in Journal
Transcriptomic Profiling of Human Limbus-Derived Stromal/Mesenchymal Stem Cells—Novel Mechanistic Insights into the Pathways Involved in Corneal Wound Healing
Previous Article in Journal
A Novel In Vivo Model for Multiplexed Analysis of Callosal Connections upon Cortical Damage
Previous Article in Special Issue
A Novel Synthetic Antibody Library with Complementarity-Determining Region Diversities Designed for an Improved Amplification Profile
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Potential Insect Growth Inhibitor against Aedes aegypti: A Bioinformatics Approach

by
Glauber V. Da Costa
1,2,3,
Moysés F. A. Neto
4,
Alicia K. P. Da Silva
3,
Ester M. F. De Sá
3,
Luanne C. F. Cancela
3,
Jeanina S. Vega
3,
Cássio M. Lobato
2,3,
Juliana P. Zuliani
5,
José M. Espejo-Román
6,
Joaquín M. Campos
6,
Franco H. A. Leite
4 and
Cleydson B. R. Santos
1,2,6,*
1
Graduate Program in Network in Pharmaceutical Innovation, Federal University of Amapá, Macapá 68902-280, AP, Brazil
2
Laboratory of Modeling and Computational Chemistry, Department of Biological and Health Sciences, Federal University of Amapá, Macapá 68902-280, AP, Brazil
3
Laboratory of Biotechnology in Natural Products, Department of Biological and Health Sciences, Federal University of Amapá, Macapá 68902-280, AP, Brazil
4
Laboratory Molecular Modeling, State University of Feira de Santana, Feira de Santana 44036-900, BA, Brazil
5
Laboratory Cellular Immunology Applied to Health, Oswaldo Cruz Foundation, FIOCRUZ Rondônia, Porto Velho 78912-000, RO, Brazil
6
Department of Pharmaceutical and Organic Chemistry, Faculty of Pharmacy, Institute of Biosanitary Research ibs, University of Granada, 18071 Granada, Spain
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(15), 8218; https://doi.org/10.3390/ijms23158218
Submission received: 19 May 2022 / Revised: 30 June 2022 / Accepted: 11 July 2022 / Published: 26 July 2022
(This article belongs to the Special Issue Application of In Silico Techniques in Drug Design)

Abstract

:
Aedes aegypti is the main vector that transmits viral diseases such as dengue, hemorrhagic dengue, urban yellow fever, zika, and chikungunya. Worldwide, many cases of dengue have been reported in recent years, showing significant growth. The best way to manage diseases transmitted by Aedes aegypti is to control the vector with insecticides, which have already been shown to be toxic to humans; moreover, insects have developed resistance. Thus, the development of new insecticides is considered an emergency. One way to achieve this goal is to apply computational methods based on ligands and target information. In this study, sixteen compounds with acceptable insecticidal activities, with 100% larvicidal activity at low concentrations (2.0 to 0.001 mg·L−1), were selected from the literature. These compounds were used to build up and validate pharmacophore models. Pharmacophore model 6 (AUC = 0.78; BEDROC = 0.6) was used to filter 4793 compounds from the subset of lead-like compounds from the ZINC database; 4142 compounds (dG < 0 kcal/mol) were then aligned to the active site of the juvenile hormone receptor Aedes aegypti (PDB: 5V13), 2240 compounds (LE < −0.40 kcal/mol) were prioritized for molecular docking from the construction of a chitin deacetylase model of Aedes aegypti by the homology modeling of the Bombyx mori species (PDB: 5ZNT), which aligned 1959 compounds (dG < 0 kcal/mol), and 20 compounds (LE < −0.4 kcal/mol) were predicted for pharmacokinetic and toxicological prediction in silico (Preadmet, SwissADMET, and eMolTox programs). Finally, the theoretical routes of compounds M01, M02, M03, M04, and M05 were proposed. Compounds M01M05 were selected, showing significant differences in pharmacokinetic and toxicological parameters in relation to positive controls and interaction with catalytic residues among key protein sites reported in the literature. For this reason, the molecules investigated here are dual inhibitors of the enzymes chitin synthase and juvenile hormonal protein from insects and humans, characterizing them as potential insecticides against the Aedes aegypti mosquito.

1. Introduction

Mosquitoes are one of the largest groups of insects (pests) with medical importance in the global control of disease transmission, and the highlight is the Aedes aegypti mosquito, the main transmitter of arboviruses such as dengue, dengue hemorrhagic fever, urban yellow fever, and chikungunya [1]. Until lately, the alarming growth of cases of these diseases in different regions of the world was restricted to emerging countries, now becoming a global concern. Currently, there is only one registered and widely distributed vaccine for yellow fever [2]. The way to reduce the spread of these viral pathologies is through the control of the transmitting mosquito; the main means is the use of chemical insecticides in the different stages of the mosquito’s development.
A class of these insecticides is called insect growth inhibitors (IGIs), which are used to control insect pests in agriculture, forestry, stored products, and public health situations [3]. IGIs act in the early stages of development (larval and pupal) and on two main neurohormones (neuropeptides), the ecdysteroids (molting hormones) and the juvenile sesquiterpenoid hormones (JHs), which have the function of controlling molting and metamorphosis [4].
However, in post-embryonic development, both hormones are related in adulthood to the regulation of reproductive maturation, release, transport, excretion, and metabolism as the main source for chitin synthesis and insect cuticle sclerotization [5,6,7]. A cuticle and attached epidermis are one of the largest and most important organs in the insect body, the function of which is to protect against physical damage, pathogen penetration, and desiccation. Associated with growth hormones, insects need to form new cuticles and periodically eliminate old ones; this involves many biochemical reactions, such as chitinization [8,9,10,11].
Chitin is a linear amino polysaccharide polymer composed of units of β-1,4-N-acetyl-d-glucosamine (GlcNAc). They are found in the cuticles of insects, crustaceans, and fungi and have the constitution of chains of α-chitin, considered the most abundant form of chitin since its form of organization is in a very stable antiparallel way. This form of sclerotization, depending on the levels of organization of the epidermis, creates structural and architectural differences where several cuticular proteins can be cross-linked with orthodiphenols and their quinones. In this way, the epidermis, composed of a single layer of underlying epidermal cells, is generated and can originate not only the differences within a species of insect but also among insects [12,13].
Eight enzymatic steps are necessary to convert the disaccharide precursor trehalose to chitin in the final step. The UDP-N-acetyl-glucosamine units are finally converted into chitin polymers by the enzyme chitin synthase, and this fact occurs at the limit of two topological spaces, from the cytosol in the lumen of the vesicles or outside the cell [10,13].
Therefore, there is great interest in insecticides that act specifically in the development of insects since their use on a large scale has provided resistance and also presented adverse reactions in other organisms, especially those with chitin (crustaceans). Therefore, through computational methods, it is possible to design compounds that meet the requirements of pesticides, as they are specifically developed for mosquito control.
The present work aims to search for new chemical entities with potential insecticidal activity with the dual inhibition of juvenile hormone protein and chitin deacetylases of Aedes aegypti, employing molecular modeling approaches by pharmacophore and docking-based virtual screening. In addition, compounds were developed through the pharmacokinetic characteristics (ADME) of Aedes aegypti physiology and toxicological properties directed at other organisms, mainly humans. The general outline that summarizes the methodological steps in this document is shown in Scheme 1 (see more details in the Materials and Methods section) [14].

2. Results and Discussion

2.1. Pharmacophore Model Generation and Evaluation

Ligand-based approaches can be useful in the prioritization of hit compounds with the same stereo-electronic features of known active inhibitors (e.g., pharmacophore models). Thus, the employment of potent compounds in pharmacophore model building is directly related to a reliable 3D representation of stereo-electronic features and pharmacophore quality [15,16,17].
First, all models presented strain energy of less than two orders of magnitude (<100 kcal/mol) and were all considered acceptable, complying with this condition. The non-dominant analysis provided by the Pareto index, the score of which is different from zero, would determine if there are dominant parameters that could affect the statistical quality of the models; this also did not occur [18]. Regarding the steric overlap for the model (Sterics), pharmacophore agreement (H_bond) and agreement between the characteristics of the compound used in the construction of the pharmacophore model in relation to the pharmacophore requirements (Mol_qry) revealed that they were statistically equivalent (Pareto = 0), and, therefore, this metric was not sufficient to select a useful pharmacophore model (see Table 1). In order to circumvent this limitation, the ability to differentiate actives compounds from false positives (decoys) was employed by the area under the ROC curve [15,19].
A test dataset with 7 inhibitors (chitin synthase—Ki < 0.001 mg·L−1 and juvenile hormone—Ki < 0.001 mg·L−1) and 350 decoys was employed in ROC curves built to calculate the respective areas under the curve (AUC-ROC) of each pharmacophore model (QFIT value; 0–100). An AUC equal to 1.0 would be found in a model with the capacity to recognize all actives before decoys, and an AUC less than 0.5 would be associated with models worse than a random selection. An AUC greater than 0.70 is associated with a moderate predictive ability [19]. According to the AUC calculation, all pharmacophore models had the ability to differentiate active compounds from false positives (AUC < 0.70; Figure 1).
Although all pharmacophore models had good predictive ability (AUC > 0.7), ROC curve analysis is a continuous metric, and thus, it is not able to classify active compounds in their early phases [20]. Thus, the Boltzmann-enhanced discrimination of the ROC curve (BEDROC) was employed to probe the pharmacophore models’ predictions to classify active compounds before decoys (Table 2).
Pharmacophore models with BEDROC > 0.5 are considered better than a random assay and useful for the early and ordered recovery of active compounds [21,22,23]. Although two pharmacophore models, models 06 and 09, showed satisfactory BEDROC values, pharmacophore model 6 was selected as it had the highest value for pharmacophore-based virtual screening.
Pharmacophore model 6, which was used in virtual screening, aims to find molecules that have in their structure the stereo-electronic requirements that are similar to the present model: two hydrophobic centers (HYs), two hydrogen bond donors (HBDs), and two hydrogen bond acceptors (HBAs). Compounds that have similarities to the physicochemical parameters of model 6 were identified since the program used was GALAHAD™ (Tripos Inc., St. Louis, MO, USA) for the alignment of very flexible compounds and different chemotypes. The adjustments were performed using a genetic flexible algorithm to align the compounds, allowing different conformations in the active site and defined by their conformers from the adjustments of the dihedral angles that were evaluated through a function of the scoring [15,24,25,26,27]. Structure–activity relationship studies for the most favorable recognition of macromolecules should consist of the aromatic presence, and the hydrophobic characteristics of hydrogen acceptance are an essential requirement for the recognition of chitin synthase [28,29]. For juvenile hormone protein, the donor and acceptor hydrogen bonds are essential parameters for binding to the protein [30,31].
In order to characterize potent inhibitors with the pharmacophore model, a potent compound and a weak compound were superimposed on pharmacophore model 6. Compounds were analyzed experimentally by Sun et al. [32], and the compound with a value of 0.001 mg·L−1 was considered potent; the weak compound was the insect growth inhibitor flucycloxuron, used as a positive control (0.05 mg·L−1) (Figure 2).
The objective of the screening was to obtain a compound that has a dual action on both chitin synthase and juvenile hormone receptors. Pharmacophore model 6 allowed us to reproduce the requirements for chitin synthase and the recognition of the juvenile hormone receptor; therefore, it is useful to prioritize the potential inhibitors of Aedes aegypti insect growth regulators.

2.2. Pharmacophore-Based Virtual Screening

In total, 214,446 compounds were prepared in multiple protonation states and multiple tautomeric forms and were collected in structure data format (SDF), which is available in multiple conformations from the Sigma-AldrichTM subset of the ZINC [33,34] database (https://zinc15.docking.org/, accessed on 12 June 2021). The database of previously prepared compounds will be able to predict compounds with better properties, such as lipophilicity and water solubility [35,36]. For each molecule in the database, up to 300 conformers were generated using the MMFF94 force field implemented with OMEGA software (Open Eye Scientific Software, Santa Fe, NM, USA, http://www.eyesopen.com, accessed on 12 June 2021) [37], and the partial atomic charges were calculated by means of the empirical Gasteiger-Hückel method. Compounds that comply with the previous conditions are better suited for molecular docking [38]. Then, 99,623 compounds overlapped pharmacopeia model 06 with partial stereo requirements (0.47 < QFIT > 98.12), with the QFIT score being considered as a parameter. A score of 100 demonstrates that a better fit between the pharmacophore model and the database has occurred; when this score is zero, no overlap has occurred [39]. As a pharmacophore-based virtual screening result, 4793 compounds (QFIT > 58.94) were prioritized with QFIT between 58.95 and 98.12 [23,39,40]. However, this approach is not useful for evaluating the spatial requirements of the binding site. Thus, molecular docking was applied to the filtered compounds from the pharmacophore-based virtual screening.

2.3. Construction of the Bombyx-Aedes 3D Model and Molecular Docking

The macromolecule deposited in the PDB-encoded 5ZNT is a chitin deacetylase 1 protein. Although the CDA1 gene of the Bombyx mori organism, without a mutation, has no crystallized ligand, Liu [41] determined the structural and biochemical characteristics of the catalytic mechanism experimentally since the enzyme is classified as a hydrolase; its catalytic activity has been recorded (EC: 3.5.1.41): [(1 → 4)-N-acetyl-β-d-glucosaminyl] (n) + n H2O = n acetate + chitosan [42].
The primary structure of the Aedes aegypti species was searched with the descriptives (Aedes aegypti + chitin deacetylase 1 protein) at https://www.ncbi.nlm.nih.gov/protein, accessed on 27 June 2021, and we found that the protein AAEL003419-PA with 535 amino acids was used for the alignment through BLAST. From similar sequences found, the most compatible organisms (with the highest percentage of identification, above 80%, and consultation coverage above 70%) were selected [43].
In the following steps, built in the MODELLER 10.1 program, the comparison, alignment, and construction of the Bombyx-Aedes model and its validation were performed. The model generated by the program was aligned in the Pymol 2.0.7 program (Warren Lyford DeLano, The PyMOL Molecular Graphics System, Schrödinger, LLC, New York, NY, USA) as shown in Figure 3, in which the RMSD alignment overlap was 0.396, from which the choice of the best atomic model with a lower score is considered adequate in its structural composition and reliability according to Rahman et al. [44].
The Bombyx-Aedes model presented a percentage of more than 90% of residues in favorable regions, demonstrating a good stereochemical quality, as analyzed by the Ramachandram plot and shown in Figure 4. The Ramachandran plot of the Bombyx-Aedes model shows the distribution of the torsion angles Φ and Ψ of the main chain of the model’s built protein [45]. The Ramachandram plot is divided into different zones: (a) those where the color intensity is least intense because of the fact that there are no interatomic clashes; (b) those where there are moderate clashes, and (c) those where the clashes are extremely severe, with the most intense coloration [46]. A Global Model Quality Estimate GMQE value of 0.5 was found, indicating a model with higher reliability.
Carbohydrate esterases of family 4 (family CE-4) are metalloenzymes; predominantly, they use Co2+ and Zn2+ ions as divalent metal cofactors in vitro. These enzymes are constituted by a structure formed by a triad of zinc and His–His–Asp bonds, acting as a catalytic domain and an acid and base. It has two oxygenated acetates that coordinate zinc in a distorted octahedral shape [42,47,48,49]. According to Dong [50], a recombinant form of peptidoglycan GlcNAc-N-acetylase from S. pneumoniae (PgdASm) binds to a zinc ion within its active site and in vivo and demonstrates that the protein PgdASm can be fully activated by the presence of the zinc ion. Since any modifications to the CE4 motifs can provide a means to change the specificity of the substrate, a preference for the metal or its reaction mechanism represents a step towards the identification of mechanism-based inhibitors for this important class of enzymes [51].

2.3.1. Selection, Validation, and Virtual Screening for Potential Dual Inhibitors

The molecular adjustment was performed using the AutoDock Vina program, and all parameter values were adjusted throughout the process to define receptors, exhaustiveness, center_xyz, and size_xyz. To generate the result in binding affinity values, the AutoDock Vina program foresees the connection of the lowest energy state at the time of the formation of the protein complex (receptor) with the ligand, with the objective of a correlation between the values of the functions and the affinity of the ligand [52,53,54,55]. Among the 4793 compounds prioritized in the virtual pharmacophore-based screening of the Zinc15 database, 4142 compounds (ΔG < 0 kcal/mol) were aligned to the binding site of the crystalline structure of the juvenile hormone protein macromolecules (code PDB: 5V13) of Aedes aegypti; its complexed binding site has the reference ligand JH3 (methyl (2E,6E)-9-[(2R)-3,3-dimethyloxiran-2-yl]-3,7-dimethylnona-2,6-dienoate). Therefore, 2240 (LE <−0.40 kcal/mol) compounds were chosen since the efficiency of the ligand selects a more concordant compound for the choice of a principal compound [56,57,58].
Molecular docking was performed using a validation protocol to determine the pose (conformation + orientation) of the ligands (JH3) within the enzyme active sites (PDB ID 5V13). The root mean square deviation (RMSD) between the reference ligand and the experimental JH3 ligand (1.19 Å) was calculated. The RMSD analysis was determined by considering the most appropriate position of the initial structure around the X-ray crystallographic complex through an experiment of the most stable position so that there is no structural change between the protein and the reference ligand when complexed with the macromolecule under study. The best molecular fit is determined by an RMSD less than or equal to 1.5 [59,60,61,62,63], as shown in Figure 5.
Kenny et al. [64] stated that it does not seem possible to quantify the absolute binding efficiency for compounds objectively. The efficiency can still be defined in a relative way by scaling the affinity differences by the corresponding molecular size differences, and the differences were prioritized for molecular docking for the chitin deacetylase model A. aegypti.
Molecular docking to the binding site of the crystalline structure of the macromolecules of the juvenile hormone protein (PDB code: 5V13) was carried out on all compounds, and, after filtering the ligand efficiently, 2240 compounds went on to the next stage of molecular docking with the model Bombyx-Aedes enzyme insect chitin deacetylase built from the Aedes aegypti mosquito.

2.3.2. Interaction Analysis of Inhibitors and Juvenile Hormone (PDB 5V13)

The juvenile hormone protein presents in its binding site the main interaction with the epoxy group and forms a hydrogen bond with the phenolic hydroxy group of Tyr129; the rest of the isoprenoid chain is surrounded by hydrophobic side chains, including those of Phe144, Tyr64, Trp53, Val65, Val68, Leu72, Leu74, Val51, and Tyr33 [37,65].
When we analyzed our ligands, important interactions were carried out with the main amino acid residues of the binding site of the Aedes aegypti juvenile hormone macromolecule (COD: 5V13). All selected ligands (M01M05) and the reference ligand complexed with the macromolecule (JH3) interacted with residues Tyr64 (except M03), Trp53, Trp50 (except M04), Try33 (except M03 and M05), and Val51 (except M03 and M04). M03 was the only one that interacted with the amino acid residues Val65 and Val68, as shown in Figure 6.

2.3.3. Interaction Analysis of Inhibitors and Enzyme Model Insect Chitin Deacetylase from Aedes aegypti

The selection of compounds took place in a double inhibition of macromolecules, where the parameter to which they were associated performed more satisfactorily in filtering the choice of the most appropriate ligand that would, in particular, function as both an inhibitor of the juvenile hormone and a growth regulator of the Aedes aegypti mosquito. In view of this approach, 2240 (LE < −0.40 kcal/mol) compounds were selected in the first stage with the juvenile hormone macromolecule. For the next stage, there was, again, a selection of the molecular docking based on the insect model Aedes aegypti chitin deacetylase, constructed from the enzyme insect chitin deacetylase deposited in the PDB (PDB code: 5ZNT); 1959 compounds (ΔG < 0 kcal/mol) were aligned and 20 compounds (LE < −0.4 kcal/mol) were prioritized for pharmacokinetic and toxicological prediction.
When we analyzed the interactions described by Liu [41], the 5ZNT protein, with the presence of the amino acid residues observed experimentally (Asp205, His 261, and Arg306; Figure 7A), and the interactions in the Bombyx-Aedes model (Figure 7B; Asp201 and 202, His257 and 262 and Arg302) maintained the same regions around zinc. This demonstrates the conservation of important amino acids for carrying out catalytic activity and, in Figure 7C, the overlap of the 5ZNT protein with the Bombyx-Aedes model, with the presence of amino acid residues closer to the metal (Asp206, His261, and 265). Figure 7D shows the overlapping of the Bombyx-Aedes model protein with the 5ZNT protein, with the presence of amino acid residues closer to the metal (His257 and 261).
In addition, the chitin deacetylase enzyme catalyzes the removal of acetyl groups from chitin and modifies this polymer during its synthesis and reorganization. The enzyme has a structure of two regions: a subdomain of NodB homology conserved in CE4 and C-terminal loops. The homology domain of NodB is a barrel (β/α)7 (residues 162–471) composed of seven parallel β strips arranged in a barrel surrounded by six α helices. The active site and the substrate connection gap are located through the upper center of the barrel (β/α)7, which contains a metal bond triad conserved throughout the CE4 family, namely, a zinc ion coordinated by Asp206 and Asp205, His261 and His265, in addition to Arg306 and canonical Arg306, responsible for stabilizing the catalytic base Asp205 [41,42].
A long crack in the substrate was also observed on the surface of the insect chitin deacetylases (CDAs) from the silkmoth Bombyx mori (BmCDA1-CAD), with only one hydrophobic residue exposed to the solvent (Tyr242) [41] that can help in the bonding of the substrate. From these interactions, the program is able to score the best connection with the most promising compounds.
When analyzing the ligands with the best scores, interactions of all ligands (M01M05) and zinc with the amino acid residues His256, His261 (except M01), and Tyr305 (except zinc) were observed. Only compounds M01 and M05 interacted with the amino acid residue Asp302, and compounds M01 and M03 were the only ones that interacted with His430 and Leu306. Compounds M01, M04, and M05 interacted with Ala303, as shown in Figure 8.
After analysis, 20 compounds with the best score values were selected for pharmaco-kinetic and toxicological screening in search of potential insecticide growth regulators. However, this approach is not useful for evaluating the permeability of compounds via the biological barriers of mosquitoes and the toxicological potential in humans. Thus, in silico pharmacokinetic and toxicological prediction was applied to the filtered compounds from the docking-based virtual screening.

2.4. In Silico Pharmacokinetic and Toxicological Properties

The development of a potential insecticide that can exercise its possible mechanisms of action by inhibiting the growth of insects must have a concentration that is sufficient to remain bioactive with its effective action. For this to occur, the insecticide needs to cross various barriers, such as the insect’s cuticle or trachea [37], as demonstrated by commercial insecticides. Their possible mechanisms of action include the inhibition of the enzyme chitin synthase located in the cuticle or the juvenile hormone protein or the insect ecdysone (diflubenzuron and buprofezin) located in the insect’s central nervous system (CNS).
Accepting that the insect’s physiology is a key point in the development of a new insecticide, the in silico pharmacokinetic properties that have been developed mainly include properties inherent to humans. However, the planning of an insecticide is based on some properties that are shared by other live organisms since the physiology is similar in some physicochemical properties [66,67,68,69]. In this study, four positive controls (diflubenzuron, DFB, teflubenzuron, TEF, flucycloxuron, FCX, and buprofezin, BUP) were used. These are insecticides that are on the market and have a proven action in inhibiting the growth of insects; they are used as references in relation to physical and chemical properties.
Costa et al. [37] described the possible ways of better penetration and action of an insecticide to gain access to the central nervous system at the adult stage of the insect, which is through the cuticle, trachea, and integument. Simon-Delson [70] and Denecker et al. [71] stated that it is common for insecticides to penetrate through the route of ingestion by the insect. However, another way is through its endocuticle in the larval period, when this tissue is in the development phase, with a large production of chitin synthase and growth hormones [72,73,74,75].
However, depending on the route that the insecticide takes, the insect’s physiological system will have several means to reduce the effect of an insecticide: for example, low penetration of the insecticide, where it does not reach sufficient concentration to exert its effective action, and where the insect has a high excretion of several insecticides, characteristic of its resistance to xenobiotics, as reported in the literature [1,76].
The planning of the insecticidal potential was developed with the screening of 20 compounds from the previous steps, of which five compounds were chosen; they had a high power of penetration through the insect pathways and sufficient quantity to inhibit the enzyme chitin synthase or the protein of the juvenile hormone or ecdysteroid hormone [77], which are the possible target receptors.
Figure S1 (supplementary material) shows the positive controls and compounds that possibly have these properties, as can be seen on the bioavailability radar generated by the SwissADME program, in which the pink area represents the ideal range for each property (lipophilicity: XLOGP3 between −0.7 and + 5.0, size: MW between 150 and 500 g/mol, polarity: TPSA between 20 and 130 Å2, solubility: log S not exceeding 6, saturation: sp3 hybridization carbon fraction >0.25, and flexibility: <9 rotating links) [78,79].
We checked the radar bioavailability graph, which has six axes of important bioavailability properties for the absorption of the insect’s body wall and CNS. Of the positive controls DFB, TEF, FCX, and BUP and compounds M01M05, as seen in Figure S1 (supplementary material) with possible insecticidal activity, BPU and M03 were the only ones that did not violate any of the parameters, the pharmacokinetic properties of which remained within the range of the graph. The positive control FCX, on the other hand, violated the polarity parameter, making it a very bad insecticide to achieve high concentrations in the target receptors.
However, DFB, TEF (positive controls), M01, M02, M04, and M05, as per Figure S1 (supplementary material), violated the saturation property with an sp3 hybridization fraction greater than 0.25, demonstrating that the design of the compound substituents could increase the complementarity of receptors/ligands, making it possible that additional protein–ligand interactions would not be accessible to a flat aromatic ring and, therefore, would improve potency and selectivity for a given target, which should mitigate the off-target effects [80].
It is clear from this data that a greater saturation would allow higher complexity of the compounds, including the presence of stereocenters. Such compounds will have access to more available chemical space, making possible the complementarity of the spatial connections of the active sites of the target protein [80], criteria obeyed by M01, M02, M04, and M05 (Figure 9). However, saturation will also affect physical properties: saturation has been correlated with solubility; the more saturated the compound, the greater the likelihood of solubility [80].
Thus, as the increase in log Po/w leads to a decrease in aqueous solubility, analyzing another method of solubility that influences the capture, distribution, transport, and, eventually, the bioactivity of a drug at the site of its actions could facilitate drug development, mainly in the ease of handling and formulation [81,82]. The FCX compound again differs in value, as expected, and cannot cross the blood–brain barrier easily. According to Table 3, the positive controls DFB and TEF are poorly soluble, whereas BUP and all planar compounds are soluble in an aqueous system. To ratify our results regarding the penetration of the blood–brain barrier and absorption by the insect cuticle, we used the data presented by the Preadmet tool, which makes biological predictions that show the same relationship of the compound’s passage to reach the studied enzymes in sufficient concentrations to exercise its effect. The TEF and FCX compounds did not cross the blood–brain barrier due to their log Po/w and S results (see Table 3).
Therefore, according to the literature, the greater the hydrophobic property, the greater the logP value. The parameters were analyzed to determine whether the insecticide can overcome the cuticle and blood–brain barrier. Another parameter of biological activity that is considered a good model for analyzing the blood–brain barrier is the apparent cell permeability property, in nm/s, of two renal cell lines from Madin-Darby dogs (MDCK) [83,84,85,86,87], as can be seen in Table 3, following our virtual screening based on hieratic ligands, where we selected, according to the criteria above, five compounds in the data set that had values within the parameters described by the literature.
The pharmacokinetic properties are directly related to the proper planning of the insecticidal potential since different substituents on a given chemical entity are used as references to potentially create more efficient compounds according to their hydrophobic, steric, and electronic characteristics. Therefore, when analyzing our compounds, we take into account the octanol–water partition coefficient to define the hydrophobic character of a substituent [88] on different platforms that are condensed by the SwissADME program [79].
The SwissADME program for better visualization allows the BOILED-Egg method of estimated brain or instinct permeation, an intuitive method to simultaneously predict two main parameters of the ADME, that is, passages across membranes and brain access (BBB). Although conceptually very simple, as it has only two physical–chemical descriptors (WLOGP and TPSA, for lipophilicity and apparent polarity), this classification model was constructed with extreme care for statistical significance and robustness [78,79,89].
The selection of compounds M01M05 was based on the structural simplicity of the derivatives and the fact that it allows the study of their effects on the cuticle and CNS. Since solubility is one of the main physicochemical properties for the development of a drug, the partition coefficient between n-octanol and water in its logarithmic form (log Po/w) could determine the assessment of the pharmacological properties of absorption, distribution, metabolism, excretion, and toxicity (ADMET) of a candidate chemical lead for its initial selection. It establishes benchmarks against which compounds synthesized during lead optimization can be evaluated [90,91].
Log Po/w, processed by the SwissADME program, as can be seen in Table 3, uses five multiple computational methods (iLOGP, XLOGP3, WLOGP, MLOGP, SILICOS-IT) for estimating and increasing the accuracy predictors of this physical–chemical property and generating the consensus log Po/w. This parameter is the arithmetic means of the values predicted by the five proposed methods based on free solvation energies in n-octanol and water, calculated by the generalized-born model, and can give an estimation of the solvent-accessible surface (GB/SA) [92,93,94,95,96]. As to the parameters established in the literature and by positive controls, only the FCX control presented a value above 5 (consensus log Po/w of 5.37); however, this compound can overcome the gastrointestinal barrier, which in the insect’s structure would be the passage through the cuticle. However, it cannot overcome the blood–brain barrier in sufficient concentrations due to its atomic nature, with three aromatic rings, together with a cyclopropane, making it possible to increase the log Po/w [97,98,99].
From the result of Table 3, four compounds that allowed the best entry into the absorption of the central nervous system (CNS) through the blood–brain barrier and one that had excellent absorption properties through the membrane were selected, as shown in the graph (see Figure S2 in the supplementary material), in which the white and yellow regions (yolk) are equivalent to the area most likely to be absorbed by the cell and most likely to access the brain (penetration), respectively. In agreement with Geldenhuys [90], who, in his study, developed a computational approach for the brain uptake of compounds in an insect model system, we concluded that in silico database filtering models of small compounds used in drug discovery were able to express the correlation between the BBB models. When analyzing the graph of the three growth-inhibiting insecticides on the market, DFB and BPU, together with compounds M01, M02, M03, and M04, are found in the yellow region (yolk) [89]. On the other hand, the insecticides TEF and FCX, and compound M05 are in the white region of HIA (gastrointestinal absorption), corroborating previous data. When performing the in silico experiment in the SwissADMET program for the 20 compounds, the program did not generate results for compounds ZINC95851157 and ZINC12447446. They cannot be seen in the graph because their results are out of reach, which can be explained by the fact that they have deficient pharmacokinetic properties, similar to the BCR-ABL-1 inhibitor. They are predicted with an accuracy of 83% (true negatives are out of range) [78].
This is a clear indication that these compounds are highly likely to be absorbed by the gastrointestinal tract and permeate the insect’s brain; the compounds TEF, FCX, and M05 did not cross the blood–brain barrier in a high concentration, in spite of our plan for the growth inhibitors’ possible mechanism by inhibiting the enzyme chitin synthase present in the cuticle. Another analysis parameter in the graph is the possibility of active efflux from the cell to expel the compound that passes through the barriers that are substrates or non-substrates of the permeability glycoproteins PGP+ (substrates) and PGP (not substrates). Blue and red dots symbolize compounds; none of the compounds were PGP+, shown as blue and red dots on the graph, and none of the compounds were PGP+. Probably they prevent the intracellular accumulation of toxic compounds, as they are eliminated outside the cell through active efflux through biological membranes [91,92,93].
The glycoprotein-P of the Aedes aegypti mosquito shows similarity in its identity to other species of insects and to humans [94]. It was observed in the literature that the baseline ATPase activity of ATP-binding cassette transporters in larvae increased insecticide efflux in the presence of organophosphate larvicides (temephos); a growth inhibitor (diflubenzuron) showed a significant increase in the toxicity of these insecticides [95]. However, Figure S2 (supplementary material) shows that all positive controls and the studied compounds were PGP since this depends on the species of insects, the structural components of the compounds, and how they are analyzed [94,96,97,98].
To design an insecticide that is developed based on certain mechanisms of action, be it the chitin synthase enzyme or the Aedes aegytpi juvenile hormone protein, it is expected that the insecticide requires a minimum concentration for the desired effect. A toxicity assessment in humans is essential for any insecticide [99]. Shown in Table 4, toxicological predictions were made using the eMolTox server (http://xundrug.cn/moltox, accessed on 30 June 2021).
EMolTox is a web server for predicting potential toxicity associated with the structure of a compound; it predicts the compound with high confidence using the aggregate predictor method (ACP). It was elaborated based on literature data from a total of 174 in vitro/in vivo experimental data sets related to toxicology. The data were used to build the model through physical–chemical descriptors and Mondrian’s shaped prediction, which offers similar prediction accuracy for both active and inactive sets, where the validity values are around 0.90, at a significance level of 0.1. This toxicological web server predicts the molecular toxicity of the target molecule based on the mechanism of action in different animal models and their organs [100,101,102,103,104,105,106,107,108,109,110].
The toxicological web server eMolTox presented the following data, of which the potential injury alerts for genotoxicity were for the positive controls FCX and BUP, in addition to injury to the liver and endocrine, whose actions are in some pathway signaling that is similar to foreheads in vitro. As for the DFB and TEF positive controls, the alerts for injuries were in the liver, nervous system, immune system, and kidney, and their actions were in modulators; all confidence levels were higher than 0.981. On the other hand, only the M04 compound was on alert for skin injury, together with alerts in the liver for the action of blocking the (OATP1) B1 transporter, which is the same for the M02 and M03 compounds. The organic anion transporter polypeptide (OATP) 1B is a detoxifying limiting protein that mediates the uptake and elimination of many endogenous and xenobiotic compounds (drugs or insecticides) through the sinusoidal membrane of liver cells [111,112,113,114,115]; M02, M03, and M04 are possibly blocking this transporter, making them molecular structures to be eliminated by living beings. M01 and M05 did not show any actions or injuries (Table 4).
The possible toxicity of insect growth inhibitors specific to several species has been analyzed in the literature. It can be found in the Scientific Opinion of the Panel on Plant Protection Products and Their Residues (PPR Panel), at the request of the PRASA Unit of the European Food Safety Authority (EFSA), on the genotoxic and carcinogenic potential of buprofezin in the context of human risk assessment. It was evaluated in a two-year study that, in mice and rats, neither the nature nor the incidence of tumors was affected by its administration [116]. However, with teflubenzuron, the main target organ of which is the liver, after short or long-term exposure in species of rats, mice, or dogs, hepatotoxicity was observed in an 18-month oncogenicity study in mice. Benign hepatocellular adenomas, observed only in this study with mice, were not considered relevant to human exposure [117].
In one study, BUP and organophosphates were tested for chromosomal aberrations in somatic and germ cells in mice, induced by high doses. The growth inhibitor was a relatively safer compound than conventional organophosphate insecticides [118]. It is undoubtedly common that the toxicity of BUPs is related to their metabolites because they are degraded by hydrolysis and photodegradation, generating metabolites 2,6-difluorobenzamide, 4-chlorophenylurea, 4-chloroacetanilide, 4-chloroaniline, and N-methyl-4-chloroaniline, of which the last three are classified as mutagenic [119,120]. It was observed in blue-green algae (Plectonema boryanum) that when the phenylurea derivatives were exposed and became highly resistant, the main metabolite was p-chlorophenylurea [1,119,121]. None of our M01M05 compounds (Figure 9) were planned without the chlorine atom, and, therefore, they have promising insecticide potential without the mutagenic effect.
The compound that showed excellent results for the standards established by the program when analyzing Table 4 was compound M05, which did not present any toxicity. Regarding hepatoxicity [122,123] and cytotoxicity, all positive controls were alerted as active, while M02 and M03 compounds were alerted as a transportation block in the liver; together with M04 (alerted for skin sensitivity), these compounds could give rise to a few adverse reactions. Teflubenzuron (TEF) is used in the delousing of salmon farms; analyses in non-target species such as lobster (Homarus gammarus) and shrimp (Palaemon elegans) have shown exoskeletal abnormalities (deformities of the claws or rigid legs); deformities of juvenile lobsters have been observed in several organs [124,125,126].

2.5. Structure–Activity Relationship of the Promising Molecule

Insect growth inhibitors possess the basic chemical structure of benzoylphenylurean compounds, which consist of three regions for SAR analysis: the benzoyl ring (A), the urea bridge (B), and the aniline (C), as shown in Figure 10. The best-known insecticide is diflunzuron, which has a structure composed of benzoyl urea, similar to sulfonylureas used as oral antidiabetics. Vila-Carriles et al. [127] found that the substitution of aromatic or lipophilic heterocyclic groups increased the potency of compounds such as glibenclamide by 100–1000 times and that the increased affinity was the result of interactions with two overlapping binding sites that were analogous to BPU. This can lead to new insecticides through the screening of bioassays, and in this study, the compound LY 131215, which does not contain a urea bridge, had the best result concerning insecticidal activity [128].
After the high-performance pharmacophoric screening of the studied compounds, only M01 presented urea in its molecular structure. M01 and its derivatives have 19 studies and patents that are mainly related to antibacterial tests, nitric oxide, DNA inhibition, neuroinflammatory treatment and autoimmune diseases, ribonucleotide reductase, tyrosinase inhibitors, and sickle cell anemia treatments [116,117,129,130,131,132,133,134].
Another widely studied compound, M05, has 26 studies related mainly to histone deacetylase, urease, metallo-β-lactamase inhibitor FEZ-1, antitumor activity, DNA synthesis, hydroxamic β-glucosyl-phenolic acid inhibitor, ribonucleotide reductase, 5-lipoxygenase inhibitors, and activity against trypanosome [59,135,136,137,138,139,140,141,142,143]. On the other hand, the M02 compound has only one study analyzing antibacterial activity [144]; in relation to M03 and M04 compounds, there are no biological studies; none of our compound cells were tested against insect-growth-inhibiting activity to transform them into a new insecticide after additional tests in vitro and in vivo.

2.6. Prediction of Synthetic Accessibility (SA) and Theoretical Synthetic Routes

The selected compounds M01M05, after hierarchical virtual screening, were submitted for evaluation of synthetic accessibility (Table 5); all five compounds were easily synthesized, with scores above 89 [145,146,147]. According to classification by score, the value of 100 is the maximum synthetic accessibility; that is, the molecule is more easily synthesized if it has a higher score (see Table 5). In this study, we propose the theoretical synthetic route of promising compounds, as can be seen in Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15 and Figure 16.

2.6.1. Theoretical Synthetic Route (Compound M01)

Compound M01 can be prepared according to the research carried out by B. Buettelman et al. [148]. Conversion of butyraldehyde I with hydroxylamine hydrochloride (NH2OH·HCl) in EtOH and H2O as solvents and using NaOH as a base will allow us to obtain the final butyraldehyde oxime—M01 (Figure 11).

2.6.2. Theoretical Synthetic Route (Compound M02)

We propose a theoretical synthetic route for M02 based on the preparation of oxime derivative V as an intermediate, with IV and NH2OH·HCl as initial reactants and NaOAc and EtOH/H2O as solvents under reflux conditions [149]. By a reduction reaction with sodium cyanoborohydride (NaBH3CN) of intermediate V, final hydroxylamine M02 will be obtained [150] (Figure 12).
Furthermore, we designed a secondary strategy based on the previous work of C. Berini et al. [151]. Firstly, nitrone VII could be furnished by condensation between hydroxylamine VI and acetaldehyde with magnesium sulfate in DCM. Eventually, the nucleophilic addition of pyrrole to intermediate VII using HCl in MeOH would originate the target compound (Figure 13).

2.6.3. Theoretical Synthetic Route (Compound M03)

Oxime derivative M03 can be synthesized by a condensation reaction (Figure 14) between aldehyde VIII and hydroxylamine hydrochloride in a 30% MeOH aqueous solution [152].
In order to improve the yield, we suggest the use of a green chemistry procedure by J. Yu et al. [153]. The final compound M03 will be formed through an aerobic oxidation of amine IX with (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO), InCl3, and acetaldoxime as catalysts, using toluene as a solvent (Figure 15).

2.6.4. Theoretical Synthetic Route (Compound M04)

We present a synthetic route to synthesize compound M04 (Figure 16) based on the formation of iodobenzaldehyde XI from iodobenzyl alcohol X using pyridinium dichromate as an oxidizing agent and DCM and Et2O as solvents. The condensation reaction between intermediate XI and hydroxylamine hydrochloride and the subsequent reduction with sodium cyanoborohydride will furnish the desired N-(iodobenzyl) hydroxylamine M04 [154].

2.6.5. Theoretical Synthetic Route (Compound M05)

We proffer a theoretical synthetic route for compound M05 using intermediate XIV as the key substrate. The initial strategy to obtain XIV consists of two steps [155]. Firstly, the coupling reaction of acid derivative XII with o-(tetrahydro-2H-pyran-2-yl) hydroxylamine (NH2OTHP) in DMF with 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide (EDC) and N,N-diisopropylethylamine (DIPEA) will provide derivative XIII. Then, the cleavage of the THP-protecting group in acidic conditions will produce benzamide XIV. The second possibility of synthesizing the key intermediate can be through palladium-catalyzed hydrogenation [156], removing the benzyl group of starting material XV. Finally, the target compound M05 can be furnished from intermediate XIV by a hydrogenation reaction using platinum/ceria-zirconia mixed oxide (Pt/CeO2-ZrO2) as a catalytic system [157] (Figure 17).

3. Materials and Methods

3.1. Dataset

A dataset of 14 benzoylphenylurean analogs was chosen to determine the sets for training and testing as they were designed and synthesized from the reference structure, with modifications at specific points to generate the model from a single base structure (Figure 18); the dataset had experimental data with excellent larvicidal activity [32], potentially for the chitin synthase and the juvenile hormone receptor.
GALAHADTM (Tripos Inc.), was employed to build ten pharmacophore models from a training set (see Figure 19) of previously related inhibition activity to chitin synthase and juvenile hormone receptors [32,120].
First, they were drawn using MarvinTM Sketch 16.9.5 software (CHEMAXON, Budapest, Hungary, https://www.chemaxon.com/, accessed on 15 March 2021), from which we selected the most reliable tautomers (pH = 6.5) for all compounds [158]. Next, 2D structures were converted to a 3D format by CONCORD implemented in the SYBYL-XTM 2.0 package (Tripos Inc.). All structures were energy-minimized using a conjugate gradient (CG) with a convergence criterion of 0.001 kcal/mol and a Tripos force field (dielectric constant ε = 80.0 and maximum number of iterations = 50,000) [159]. Partial atomic charges were calculated using the Gasteiger-Hückel method [160], available on SYBYL-X 2.0TM (Tripos Inc.).
The study of larvicidal activity was carried out by Sun [32] et al., from which the fourteen compounds with 100% larvicidal activity at low concentrations (2 to 0.001 mg·L−1) were selected, seven of which were used to determine the training set and seven for the test set. They were used randomly in the construction of the pharmacophore models, while the other compounds with the same criteria were used in the validation of the pharmacophore model.

3.2. Pharmacophore Model Generation

The training datasets were flexibly aligned to each other in order to build hypermolecular alignments that map common pharmacophore features. The genetic algorithm employed in this step starts with 45 conformations (population size) of each benzoylphenylurean derivative and evolves through a maximum of 90 generations by means of default genetic operators (mutation rate = 0.2 and crossover rate = 0.2) available in the GALAHADTM module in SYBYL-XTM 2.0 [161] (Tripos Inc.).

3.3. Pharmacophore Model Evaluation

Pharmacophore models with the lowest energy (<100.0 kcal/mol) were selected and analyzed for the results of the Pareto score. To determine the best pharmacophore model and its discriminatory power, the models were used to build decoys that were removed from the DUD-E server, March (http://dude.docking.org/, accessed on 15 March 2021), [162] and executed in the SigmaPlotTM program v. 12.0 (SYSTAT, Leland Wilkinson, Chicago, IL, USA) by calculating the area under the receiver operating characteristics curve (AUC-ROC curve); it was possible to visualize the recognition of assets and decoys for better validation of the models through the true-positive rate versus the false-positive rate.
The Boltzmann-enhanced discrimination of the ROC curve (BEDROC) was employed to discriminate models with AUC > 0.7 according to early recognition of actives using the ROCKER server [163] (www.jyu.fi/rocker, accessed on 15 March 2021). Additionally, BEDROC is oriented by a search algorithm (α = 16.10) that corresponds to the initial 10% of the database. The pharmacophore model that demonstrated early recognition (BEDROC > 0.5) was employed in the virtual screening step.

3.4. Pharmacophore-Based Virtual Screening

The best pharmacophore model was automatically converted to the UNITYTM 3D query that was employed to virtually screen the Sigma-Aldrich subset from the ZINC database [34] (https://zinc15.docking.org/, accessed on 15 March 2021). Compounds were ranked according to QFIT values, and those with values higher than mean plus two times standard deviation (Equation (1)) were selected for molecular docking.
X = x + 2×σ
  • X = QFIT value
  • x = Average
  • σ = Standard deviation
Equation (1): Mathematical equation of the mean plus two times the standard deviation.

3.5. Homology Modeling, Molecular Docking, and Virtual Screening for Potential Dual Modulators

The primary sequence of Aedes aegypti species was collected in the National Center for Biotechnology Information (NCBI), with the descriptives (Aedes aegypti + chitin deacetylase 1 protein) at https://www.ncbi.nlm.nih.gov/protein, accessed on 27 June 2021. We found that the protein AAEL003419-PA was used as a primary protein; the BLASTTM, https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 27 June 2021 [164] performed and compared the primary biological sequences and found the most similar and statistically significant ones in the PDB (Protein Data Bank).
The most similar macromolecules were introduced and analyzed in the MODELLER 10.1 program, which, with the primary structure of Aedes aegypti, performed the comparison, alignment, construction, and validation of the model, called Bombyx-Aedes.
Using the SWISS-MODEL server [165,166,167], structure assessment (https://swissmodel.expasy.org/assess, accessed on 27 June 2021) was employed to analyze the model created for A. aegypti chitin deacetylase using the global model quality estimation (GMQE) score on the accuracy of the constructed model and the model used [168]; the prioritized model had a GMQE score of <0.50 [169]. Furthermore, Ramachandran plots were used to evaluate the energetically favored regions of the model [170]. In the Pymol program, the 3D alignment between the structures was performed, which collected the coordinates of the alpha-carbons of each of the residues in this analysis.
The conformational search and scoring evaluation were performed using AutoDock Vina 1.1.2 (Scripps Research Institute, Olson Laboratory, La Jolla, CA, USA). The ability to generate satisfactory solutions to scoring functions was probed according to the root mean square deviation value (RMSD < 2Å) [171,172].
The compounds prioritized in pharmacophore-based virtual screening were employed in molecular docking, and in silico approaches of the Aedes aegypti juvenile hormone receptor (PDB: 5V13) were performed according to previous studies [173,174]. After the molecular docking, ligand efficiency (LE) was calculated through the binding affinity of compounds employed in docking studies (Equation (2)), and the compounds with LE < −0.4 kcal/mol were prioritized to dockig with the A. aegypti chitin deacetylase receptor model.
LE = −ΔG/N
  • LE = Ligand efficiency
  • ΔG = Binding affinity
  • N = Number of heavy atoms
Equation (2). Mathematical equation of ligand efficiency.
The ligand efficiency (LE) was calculated for compounds employed in docking with the A. aegypti chitin deacetylase receptor, and compounds with LE < −0.4 kcal/mol were prioritized for pharmacokinetic and toxicological analysis.

3.6. In Silico Pharmacokinetic and Toxicological Properties

For the analysis of pharmacokinetic properties, the best compounds from the previous step were selected, and we used the Preadmet (https://preadmet.bmdrc.kr/, accessed on 27 June 2021) and SwissADME tools [79,175] and the proficient method of cerebral or intestinal permeation Estimation D (BOILED-Egg) [78] (http://www.swissadme.ch/index.php, accessed on 27 June 2021). The analysis of the toxicity profile of the compounds was evaluated using the online platform eMolTox (http://xundrug.cn/moltox, accessed on 27 June 2021).

3.7. Synthetic Accessibility (SA) Prediction

The prediction of the synthetic accessibility (SA) of the compounds selected by the hierarchical virtual screening was performed using the AMBIT program (http://ambit.sourceforge.net/reactor.html, accessed on 10 July 2021), which implements the tautomer generations by the algorithm method of combinatorics, exhaustively calculating all tautomeric forms of a given organic compound by molecular and stereochemical complexity and considering the presence of fused and bridged systems to suggest the canonical tautomer by an energy base of classification by different structural and topological characteristics combined with an additive scheme using weighted molecular descriptors. By generating a score ranging from 0 to 100 (maximum synthetic accessibility), the molecule is more easily synthesized when the score is closer to 100 [176,177,178,179].

4. Conclusions

This study has established the planning of insect-growth-inhibiting compounds with potential insecticidal activity that have the possible mechanism of action of dual inhibition of the juvenile hormone protein; the study has also established a chitin deacetylase model of Aedes aegypti through molecular modeling approaches via pharmacophore and virtual screening based on molecular docking by predicting molecular target binding. In addition, compounds were developed through the pharmacokinetic characteristics (ADME) of the physiology of the Aedes aegypti mosquito, and toxicological properties directed at other organisms, mainly humans, as well as the synthetic theoretical routes of compounds M01, M02, M03, M04, and M05 have been proposed.
Finally, compounds M01M05 were selected, showing parameters similar and superior to those observed for positive controls and also significant differences in the interaction with key residues in each catalytic protein site reported in the literature. Accordingly, the molecules investigated here are dual inhibitors of the enzymes chitin synthase and juvenile hormonal protein (from insects), characterizing them as potential insecticides against the mosquito Aedes aegypti. With these promising results, our research group intends to carry out biological tests in vivo with the compounds obtained in this work to confirm the in silico predictions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23158218/s1.

Author Contributions

Conceptualization, G.V.D.C., J.M.E.-R., F.H.A.L. and C.B.R.S.; Data curation, A.K.P.D.S. and C.M.L.; Formal analysis, M.F.A.N., A.K.P.D.S., E.M.F.D.S., L.C.F.C., C.M.L., J.M.C. and C.B.R.S.; Funding acquisition, G.V.D.C., J.P.Z., J.M.C. and C.B.R.S.; Investigation, A.K.P.D.S., E.M.F.D.S., L.C.F.C., J.S.V., C.M.L., J.M.E.-R. and F.H.A.L.; Methodology, G.V.D.C., M.F.A.N., E.M.F.D.S., L.C.F.C., J.S.V., J.M.E.-R., J.M.C. and F.H.A.L.; Project administration, C.B.R.S.; Software, M.F.A.N.; Supervision, F.H.A.L. and C.B.R.S.; Validation, J.P.Z.; Visualization, A.K.P.D.S., J.S.V., J.P.Z. and J.M.C.; Writing—original draft, G.V.D.C., M.F.A.N., J.M.E.-R. and F.H.A.L.; Writing—review and editing, J.P.Z. and C.B.R.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We gratefully acknowledge the computational support provided by the Laboratory of Modeling and Computational Chemistry (LMQC) and the Graduate Program in-Network in Pharmaceutical Innovation of the Department of Biological Sciences at the Federal University of Amapá (UNIFAP/Macapá-Brazil) and thank the Laboratory of Molecular Modeling of the State University of Feira de Santana (Bahia-Brazil). We would also like to thank the Laboratory of Cellular Immunology Applied to Health of the Oswaldo Cruz Foundation (FIOCRUZ); the Department of Pharmaceutical and Organic Chemistry, Faculty of Pharmacy of the University of Granada (Spain); and the Researcher Assistance Program—PAPESQ/UNIFAP for financial support.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Georghiou, G.P.; Saito, T. Pest Resistance to Pesticides; Springer Science & Business Media: Boston, MA, USA, 1983; ISBN 9788578110796. [Google Scholar]
  2. Monath, T.P. Yellow Fever Vaccine. Expert Rev. Vaccines 2005, 4, 553–574. [Google Scholar] [CrossRef] [PubMed]
  3. Mondal, K.A.M.S.H.; Parween, S. Insect Growth Regulators and Their Potential in the Management of Stored-Product Insect Pests. Integr. Pest Manag. Rev. 2000, 5, 255–295. [Google Scholar] [CrossRef]
  4. Gilbert, L.I. (Ed.) The Juvenile Hormones; PLENUM PRE.; Springer US: Boston, MA, USA, 1976; ISBN 978-1-4684-7949-2. [Google Scholar]
  5. Wing, K.D.; Sparks, T.C.; Lovell, V.M.; Levinson, S.O.; Hammock, B.D. The Distribution of Juvenile Hormone Esterase and Its Interrelationship with Other Proteins Influencing Juvenile Hormone Metabolism in the Cabbage Looper, Trichoplusia Ni. Insect Biochem. 1981, 11, 473–485. [Google Scholar] [CrossRef]
  6. Hoffmann, K.H.; Lorenz, M.W. Recent Advances in Hormones in Insect Pest Control. Phytoparasitica 1998, 26, 323–330. [Google Scholar] [CrossRef]
  7. Noriega, F.G. Nutritional Regulation of JH Synthesis: A Mechanism to Control Reproductive Maturation in Mosquitoes? Insect Biochem. Mol. Biol. 2004, 34, 687–693. [Google Scholar] [CrossRef]
  8. Hepburn, H.R.; Roberts, D.C. Stiffness and Tanning of Sclerites. J. Insect Physiol. 1975, 21, 1741–1746. [Google Scholar] [CrossRef]
  9. Londershausen, M.; Kammann, V.; Spindler-Barth, M.; Spindler, K.D.; Thomas, H. Chitin Synthesis in Insect Cell Lines. Insect Biochem. 1988, 18, 631–636. [Google Scholar] [CrossRef]
  10. Kramer, K.J.; Koga, D. Insect Chitin. Physical State, Synthesis, Degradation and Metabolic Regulation. Insect Biochem. 1986, 16, 851–877. [Google Scholar] [CrossRef]
  11. Andersen, S.O. Insect Cuticular Sclerotization: A Review. Insect Biochem. Mol. Biol. 2010, 40, 166–178. [Google Scholar] [CrossRef]
  12. Merzendorfer, H. Insect Chitin Synthases: A Review. J. Comp. Physiol. B Biochem. Syst. Environ. Physiol. 2006, 176, 1–15. [Google Scholar] [CrossRef]
  13. Doucet, D.; Retnakaran, A. Insect Chitin. Metabolism, Genomics and Pest Management, 1st ed.; Elsevier Ltd.: Amsterdam, The Netherlands, 2012; Volume 43, ISBN 9780123915009. [Google Scholar]
  14. Neto, M.F.d.A.; dos Santos, C.B.R.; Magalhães-Junior, J.T.; Leite, F.H.A. Identification of Novel Aedes Aegypti Odorant-Binding Protein 1 Modulators by Ligand and Structure-Based Approaches and Bioassays. J. Biomol. Struct. Dyn. 2020, 40, 117–129. [Google Scholar] [CrossRef] [PubMed]
  15. Shepphird, J.K.; Clark, R.D. A Marriage Made in Torsional Space: Using GALAHAD Models to Drive Pharmacophore Multiplet Searches. J. Comput.-Aided Mol. Des. 2006, 20, 763–771. [Google Scholar] [CrossRef] [PubMed]
  16. Long, W.; Liu, P.; Li, Q.; Xu, Y.; Gao, J. 3D-QSAR Studies on a Class of IKK-2 Inhibitors with GALAHAD Used to Develop Molecular Alignment Models. QSAR Comb. Sci. 2008, 27, 1113–1119. [Google Scholar] [CrossRef]
  17. Leite, F.H.A.; Froes, T.Q.; da Silva, S.G.; de Souza, E.I.M.; Vital-Fujii, D.G.; Trossini, G.H.G.; Pita, S.S.d.R.; Castilho, M.S. An Integrated Approach towards the Discovery of Novel Non-Nucleoside Leishmania Major Pteridine Reductase 1 Inhibitors. Eur. J. Med. Chem. 2017, 132, 322–332. [Google Scholar] [CrossRef] [PubMed]
  18. Nicolaou, C.A.; Brown, N. Multi-Objective Optimization Methods in Drug Design. Drug Discov. Today Technol. 2013, 10, e427–e435. [Google Scholar] [CrossRef]
  19. Bittencourt, J.A.H.M.; Neto, M.F.A.; Lacerda, P.S.; Bittencourt, R.C.V.S.; Silva, R.C.; Lobato, C.C.; Silva, L.B.; Leite, F.H.A.; Zuliani, J.P.; Rosa, J.M.C.; et al. In Silico Evaluation of Ibuprofen and Two Benzoylpropionic Acid Derivatives with Potential Anti-Inflammatory Activity. Molecules 2019, 24, 1476. [Google Scholar] [CrossRef] [Green Version]
  20. Zhao, W.; Hevener, K.E.; White, S.W.; Lee, R.E.; Boyett, J.M. A Statistical Framework to Evaluate Virtual Screening. BMC Bioinform. 2009, 13, 225. [Google Scholar] [CrossRef] [Green Version]
  21. Truchon, J.F.; Bayly, C.I. Evaluating Virtual Screening Methods: Good and Bad Metrics for the “Early Recognition” Problem. J. Chem. Inf. Model. 2007, 47, 488–508. [Google Scholar] [CrossRef]
  22. Chaput, L.; Martinez-Sanz, J.; Saettel, N.; Mouawad, L. Benchmark of Four Popular Virtual Screening Programs: Construction of the Active/Decoy Dataset Remains a Major Determinant of Measured Performance. J. Cheminform. 2016, 8, 56. [Google Scholar] [CrossRef] [Green Version]
  23. Mascarenhas, A.M.S.; de Almeida, R.B.M.; de Araujo Neto, M.F.; Mendes, G.O.; da Cruz, J.N.; dos Santos, C.B.R.; Botura, M.B.; Leite, F.H.A. Pharmacophore-Based Virtual Screening and Molecular Docking to Identify Promising Dual Inhibitors of Human Acetylcholinesterase and Butyrylcholinesterase. J. Biomol. Struct. Dyn. 2021, 39, 6021–6030. [Google Scholar] [CrossRef]
  24. Brooijmans, N.; Kuntz, I.D. Molecular Recognition and Docking Algorithms. Annu. Rev. Biophys. Biomol. Struct. 2003, 32, 335–373. [Google Scholar] [CrossRef] [PubMed]
  25. Warren, G.L.; Andrews, C.W.; Capelli, A.M.; Clarke, B.; LaLonde, J.; Lambert, M.H.; Lindvall, M.; Nevins, N.; Semus, S.F.; Senger, S.; et al. A Critical Assessment of Docking Programs and Scoring Functions. J. Med. Chem. 2006, 49, 5912–5931. [Google Scholar] [CrossRef] [PubMed]
  26. Yang, S.Y. Pharmacophore Modeling and Applications in Drug Discovery: Challenges and Recent Advances. Drug Discov. Today 2010, 15, 444–450. [Google Scholar] [CrossRef]
  27. Huang, S.Y.; Grinter, S.Z.; Zou, X. Scoring Functions and Their Evaluation Methods for Protein-Ligand Docking: Recent Advances and Future Directions. Phys. Chem. Chem. Phys. 2010, 12, 12899–12908. [Google Scholar] [CrossRef]
  28. Kamst, E.; Zegelaar-jaarsveld, K.; Marel, G.A. Van Der Chemical Synthesis of N-Acetylglucosamine Derivatives and Their Use as Glycosyl Acceptors by the Mesorhizobium Loti Chitin Oligosaccharide Synthase NodC. Carbohydr. Res. 1999, 321, 176–189. [Google Scholar] [CrossRef]
  29. Rohidas, R.; Sanjaykumar, P. CHEMISTRY Pharmacophore-Based 3D-QSAR Study of Fungal Chitin Synthase Inhibitors. Med. Chem. Res. 2013, 22, 1762–1772. [Google Scholar] [CrossRef]
  30. Liszekova, D.; Polakovic, M.; Ben, M.; Farkas, R. Molecular Determinants of Juvenile Hormone Action as Revealed by 3D QSAR Analysis in Drosophila. PLoS ONE 2009, 4, e6001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Latha, R.S.; Padmanabhan, R.V.J.; Comsia, Á.C.Á. CHEMISTRY 3D-QSAR Studies on the Biological Activity of Juvenile Hormone Mimetic Compounds for Culex Pipiens Larvae. Med. Chem. Res. 2013, 22, 5948–5960. [Google Scholar] [CrossRef]
  32. Sun, R.; Li, Y.; Lü, M.; Xiong, L.; Wang, Q. Synthesis, Larvicidal Activity, and SAR Studies of New Benzoylphenylureas Containing Oxime Ether and Oxime Ester Group. Bioorg. Med. Chem. Lett. 2010, 20, 4693–4699. [Google Scholar] [CrossRef]
  33. Irwin, J.J.; Shoichet, B.K. ZINC—A Free Database of Commercially Available Compounds for Virtual Screening. J. Chem. Inf. Model. 2005, 45, 177–182. [Google Scholar] [CrossRef] [Green Version]
  34. Sterling, T.; Irwin, J.J. ZINC 15—Ligand Discovery for Everyone. J. Chem. Inf. Model. 2015, 55, 2324–2337. [Google Scholar] [CrossRef] [PubMed]
  35. Kralj, S.; Jukič, M.; Bren, U. Commercial SARS-CoV-2 Targeted, Protease Inhibitor Focused and Protein–Protein Interaction Inhibitor Focused Molecular Libraries for Virtual Screening and Drug Design. Int. J. Mol. Sci. 2022, 23, 393. [Google Scholar] [CrossRef] [PubMed]
  36. Kolarič, A.; Jukič, M.; Bren, U. Novel Small-Molecule Inhibitors of the SARS-CoV-2 Spike Protein Binding to Neuropilin 1. Pharmaceuticals 2022, 15, 165. [Google Scholar] [CrossRef] [PubMed]
  37. Da Costa, G.V.; Ferreira, E.F.B.; de Ramos, R.S.; Da Silva, L.B.; De Sá, E.M.F.; Da Silva, A.K.P.; Lobato, C.M.; Souto, R.N.P.; da Silva de P, C.H.T.; Federico, L.B.; et al. Hierarchical Virtual Screening of Potential Insectides Inhibitors of Acetylcholinesterase and Juvenile Hormone from Temephos. Pharmaceuticals 2019, 12, 61. [Google Scholar] [CrossRef] [Green Version]
  38. Jukič, M.; Janežič, D.; Bren, U. Ensemble Docking Coupled to Linear Interaction Energy Calculations for Identification of Coronavirus Main Protease (3CLpro) Non-Covalent Small-Molecule Inhibitors. Molecules 2020, 25, 5808. [Google Scholar] [CrossRef]
  39. Pinto, V.d.S.; Araújo, J.S.C.; Silva, R.C.; da Costa, G.V.; Cruz, J.N.; Neto, M.F.D.A.; Campos, J.M.; Santos, C.B.R.; Leite, F.H.A.; Junior, M.C.S. In Silico Study to Identify New Antituberculosis Molecules from Natural Sources by Hierarchical Virtual Screening and Molecular Dynamics Simulations. Pharmaceuticals 2019, 12, 36. [Google Scholar] [CrossRef] [Green Version]
  40. Vyas, V.K.; Ghate, M.; Goel, A. Pharmacophore Modeling, Virtual Screening, Docking and in Silico ADMET Analysis of Protein Kinase B (PKB β) Inhibitors. J. Mol. Graph. Model. 2013, 42, 17–25. [Google Scholar] [CrossRef]
  41. Liu, L.; Zhou, Y.; Qu, M.; Qiu, Y.; Guo, X.; Zhang, Y.; Liu, T.; Yang, J.; Yang, Q. Structural and Biochemical Insights into the Catalytic Mechanisms of Two Insect Chitin Deacetylases of the Carbohydrate Esterase 4 Family. J. Biol. Chem. 2019, 294, 5774–5783. [Google Scholar] [CrossRef] [Green Version]
  42. Blair, D.E.; Hekmat, O.; Schüttelkopf, A.W.; Shrestha, B.; Tokuyasu, K.; Withers, S.G.; Van Aalten, D.M.F. Structure and Mechanism of Chitin Deacetylase from the Fungal Pathogen Colletotrichum Lindemuthianum. Biochemistry 2006, 45, 9416–9426. [Google Scholar] [CrossRef] [Green Version]
  43. Nene, V.; Wortman, J.R.; Lawson, D.; Haas, B.; Kodira, C.; Tu, Z.; Loftus, B.; Xi, Z.; Megy, K.; Grabherr, M.; et al. Genome Sequence of Aedes Aegypti, a Major Arbovirus Vector. Science 2007, 316, 1718–1723. [Google Scholar] [CrossRef] [Green Version]
  44. Rahman, T.; Nagar, A.; Duffy, E.B.; Okuda, K.; Silverman, N.; Harton, J.A. NLRP3 Sensing of Diverse Inflammatory Stimuli Requires Distinct Structural Features. Front. Immunol. 2020, 11, 7–10. [Google Scholar] [CrossRef] [PubMed]
  45. Hasan, N.; Oany, A.R.; Hirashima, A. Comparative Structure Analysis of Tyramine- β -Hydroxylase from Fruit Fly and ADME/T-Based Profiling of 1-Arylimidazole-2 (3 H) -Thiones as Potential Inhibitors. Int. J. Anal. Bio-Sci. 2016, 4, 10. [Google Scholar]
  46. Carugo, O.; Djinovic Carugo, K. Half a Century of Ramachandran Plots. Acta Crystallogr. Sect. D Biol. Crystallogr. 2013, 69, 1333–1341. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Urch, J.E.; Hurtado-Guerrero, R.; Brosson, D.; Liu, Z.; Eijsink, V.G.H.; Texier, C.; Van Aalten, D.M.F. Structural and Functional Characterization of a Putative Polysaccharide Deacetylase of the Human Parasite Encephalitozoon Cuniculi. Protein Sci. 2009, 18, 1197–1209. [Google Scholar] [CrossRef] [Green Version]
  48. Little, D.J.; Poloczek, J.; Whitney, J.C.; Robinson, H.; Nitz, M.; Howell, P.L. The Structure- and Metal-Dependent Activity of Escherichia Coli PgaB Provides Insight into the Partial De-N-Acetylation of Poly-β-1,6-N-Acetyl-d-Glucosamine. J. Biol. Chem. 2012, 287, 31126–31137. [Google Scholar] [CrossRef] [Green Version]
  49. Andrés, E.; Albesa-Jové, D.; Biarnés, X.; Moerschbacher, B.M.; Guerin, M.E.; Planas, A. Structural Basis of Chitin Oligosaccharide Deacetylation. Angew. Chem.-Int. Ed. 2014, 53, 6882–6887. [Google Scholar] [CrossRef]
  50. Dong, M.D.; Urch, J.E.; Ten Cate, J.M.; Rao, V.A.; Van Aalten, D.M.F.; Crielaard, W. Streptococcus Mutans SMU.623c Codes for a Functional, Metal-Dependent Polysaccharide Deacetylase That Modulates Interactions with Salivary Agglutinin. J. Bacteriol. 2009, 91, 394–402. [Google Scholar] [CrossRef] [Green Version]
  51. Taylor, E.J.; Gloster, T.M.; Turkenburg, J.P.; Vincent, F.; Brzozowski, A.M.; Dupont, C.; Shareck, F.; Centeno, M.S.J.; Prates, J.A.M.; Puchart, V.; et al. Structure and Activity of Two Metal Ion-Dependent Acetylxylan Esterases Involved in Plant Cell Wall Degradation Reveals a Close Similarity to Peptidoglycan Deacetylases. J. Biol. Chem. 2006, 281, 10968–10975. [Google Scholar] [CrossRef] [Green Version]
  52. Westhead, D.R.; Clark, D.E.; Murray, C.W. A Comparison of Heuristic Search Algorithms for Molecular Docking. J. Comput.-Aided Mol. Des. 1997, 11, 209–228. [Google Scholar] [CrossRef]
  53. Eldridge, M.D.; Murray, C.W.; Auton, T.R.; Paolini, G.V.; Mee, R.P. Empirical Scoring Functions: I. The Development of a Fast Empirical Scoring Function to Estimate the Binding Affinity of Ligands in Receptor Complexes. J. Comput.-Aided Mol. Des. 1997, 11, 425–445. [Google Scholar] [CrossRef]
  54. Morris, G.M.; Goodsell, D.S.; Halliday, R.S.; Huey, R.; Hart, W.E.; Belew, R.K.; Olson, A.J. Automated Docking Using a Lamarckian Genetic Algorithm and an Empirical Binding Free Energy Function. J. Comput. Chem. 1998, 19, 1639–1662. [Google Scholar] [CrossRef] [Green Version]
  55. Nguyen, N.T.; Nguyen, T.H.; Pham, T.N.H.; Huy, N.T.; Bay, M.V.; Pham, M.Q.; Nam, P.C.; Vu, V.V.; Ngo, S.T. Autodock Vina Adopts More Accurate Binding Poses but Autodock4 Forms Better Binding Affinity. J. Chem. Inf. Model. 2020, 60, 204–211. [Google Scholar] [CrossRef] [PubMed]
  56. Valdés-Tresanco, M.S.; Valdés-Tresanco, M.E.; Valiente, P.A.; Moreno, E. AMDock: A Versatile Graphical Tool for Assisting Molecular Docking with Autodock Vina and Autodock4. Biol. Direct 2020, 15, 12. [Google Scholar] [CrossRef] [PubMed]
  57. Costa Júnior, D.B.; Araújo, J.S.C.; de Mattos Oliveira, L.; Neri, F.S.M.; Moreira, P.O.L.; Taranto, A.G.; Fonseca, A.L.; de Pilla Varotti, F.; Leite, F.H.A. Identification of Novel Antiplasmodial Compound by Hierarquical Virtual Screening and in Vitro Assays. J. Biomol. Struct. Dyn. 2020, 39, 3378–3386. [Google Scholar] [CrossRef]
  58. Hopkins, A.L.; Keserü, G.M.; Leeson, P.D.; Rees, D.C.; Reynolds, C.H. The Role of Ligand Efficiency Metrics in Drug Discovery. Nat. Rev. Drug Discov. 2014, 13, 105–121. [Google Scholar] [CrossRef]
  59. Ul-Haq, Z.; Wadood, A.; Uddin, R. CoMFA and CoMSIA 3D-QSAR Analysis on Hydroxamic Acid Derivatives as Urease Inhibitors. J. Enzym. Inhib. Med. Chem. 2009, 24, 272–278. [Google Scholar] [CrossRef]
  60. Yusuf, D.; Davis, A.M.; Kleywegt, G.J.; Schmitt, S. An Alternative Method for the Evaluation of Docking Performance: RSR vs RMSD. J. Chem. Inf. Model. 2008, 48, 1411–1422. [Google Scholar] [CrossRef]
  61. Kores, K.; Lešnik, S.; Bren, U.; Janežič, D.; Konc, J. Discovery of Novel Potential Human Targets of Resveratrol by Inverse Molecular Docking. J. Chem. Inf. Model. 2019, 59, 2467–2478. [Google Scholar] [CrossRef]
  62. Furlan, V.; Bren, U. Insight into Inhibitory Mechanism of Pde4d by Dietary Polyphenols Using Molecular Dynamics Simulations and Free Energy Calculations. Biomolecules 2021, 11, 479. [Google Scholar] [CrossRef]
  63. Lešnik, S.; Bren, U. Mechanistic Insights into Biological Activities of Polyphenolic Compounds from Rosemary Obtained by Inverse Molecular Docking. Foods 2022, 11, 67. [Google Scholar] [CrossRef]
  64. Kenny, P.W. The Nature of Ligand Efficiency. J. Cheminform. 2019, 11, 8. [Google Scholar] [CrossRef] [PubMed]
  65. Kim, I.H.; Pham, V.; Jablonka, W.; Goodman, W.G.; Ribeiro, J.M.C.; Andersen, J.F. A Mosquito Hemolymph Odorant-Binding Protein Family Member Specifically Binds Juvenile Hormone. J. Biol. Chem. 2017, 292, 15329–15339. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Wu, L.; Aster, J.C.; Blacklow, S.C.; Lake, R.; Artavanis-Tsakonas, S.; Griffin, J.D. MAML1, a Human Homologue of Drosophila Mastermind, Is a Transcriptional Co-Activator for NOTCH Receptors. Nat. Genet. 2000, 26, 484–489. [Google Scholar] [CrossRef] [PubMed]
  67. Borges, N.M.; Sartori, G.R.; Ribeiro, J.F.R.; Rocha, J.R.; Martins, J.B.L.; Montanari, C.A.; Gargano, R. Similarity Search Combined with Docking and Molecular Dynamics for Novel HAChE Inhibitor Scaffolds. J. Mol. Model. 2018, 24, 41. [Google Scholar] [CrossRef]
  68. Lobato Rodrigues, A.B.; Martins, R.L.; Rabelo, É.D.M.; Tomazi, R.; Santos, L.L.; Brandão, L.B.; Faustino, C.G.; Ferreira Farias, A.L.; Dos Santos, C.B.R.; de Castro Cantuária, P.; et al. Development of nano-emulsions based on Ayapana triplinervis essential oil for the control of Aedes aegypti larvae. PLoS ONE 2021, 16, e0254225. [Google Scholar] [CrossRef]
  69. Sinha, S.; Pant, K.; Pal, M.; Sharma, P.; Anand, J.; Pant, B.; Verma, D.; Kumar, R.; Tripathi, Y.; Birla, K.; et al. Molecular Docking Studies of Phytocompounds in Common Weeds with AChE Proteins of Aphids and Beetles. J. Biol. Act. Prod. Nat. 2020, 10, 18–33. [Google Scholar] [CrossRef]
  70. Simon-Delso, N.; Amaral-Rogers, V.; Belzunces, L.P.; Bonmatin, J.M.; Chagnon, M.; Downs, C.; Furlan, L.; Gibbons, D.W.; Giorio, C.; Girolami, V.; et al. Systemic Insecticides (Neonicotinoids and Fipronil): Trends, Uses, Mode of Action and Metabolites. Environ. Sci. Pollut. Res. 2015, 22, 5–34. [Google Scholar] [CrossRef]
  71. Denecke, S.; Swevers, L.; Douris, V.; Vontas, J. How Do Oral Insecticidal Compounds Cross the Insect Midgut Epithelium? Insect Biochem. Mol. Biol. 2018, 103, 22–35. [Google Scholar] [CrossRef]
  72. Bennet-clark, H.C. Active Control of the Mechancal Ixsect Endocuticle Properties. Ins. Physiol. 1967, 8, 627–633. [Google Scholar] [CrossRef]
  73. Sowa, B.A. An in vitro system for the quantitative measurement of CHITIN synthesis in the cockroach: Inhibition by TH 6040 and polyoxin D. Insect Biochem. 1975, 5, 855–859. [Google Scholar] [CrossRef]
  74. Verloop, A.; Ferrell, C.D. Benzoylphenyl Ureas—A New Group of Larvicides Interfering with Chitin Deposition. Am. Chem. Soc. 1977, 19111, 237–270. [Google Scholar] [CrossRef]
  75. Tecle, B.; Ruzo, L.O.; Casida, J.E. Radiosynthesis of [Benzoyl-3,4,5-3H] diflubenzuron by a Route Applicable to Other High-Potency Insect Growth Regulators. J. Agric. Food Chem. 1988, 36, 178–180. [Google Scholar] [CrossRef]
  76. Panini, M.; Manicardi, G.C.; Moores, G.D.; Mazzoni, E. An Overview of the Main Pathways of Metabolic Resistance in Insects. Invertebr. Surviv. J. 2016, 13, 326–335. [Google Scholar]
  77. Sparks, T.C.; Nauen, R. IRAC: Mode of Action Classification and Insecticide Resistance Management. Pestic. Biochem. Physiol. 2015, 121, 122–128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Daina, A.; Zoete, V. A BOILED-Egg to Predict Gastrointestinal Absorption and Brain Penetration of Small Molecules. ChemMedChem 2016, 11, 1117–1121. [Google Scholar] [CrossRef] [Green Version]
  79. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A Free Web Tool to Evaluate Pharmacokinetics, Drug-Likeness and Medicinal Chemistry Friendliness of Small Molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef] [Green Version]
  80. Lovering, F.; Bikker, J.; Humblet, C. Escape from Flatland: Increasing Saturation as an Approach to Improving Clinical Success. J. Med. Chem. 2009, 52, 6752–6756. [Google Scholar] [CrossRef]
  81. Delaney, J.S. ESOL: Estimating Aqueous Solubility Directly from Molecular Structure. J. Chem. Inf. Comput. Sci. 2004, 44, 1000–1005. [Google Scholar] [CrossRef]
  82. Ali, J.; Camilleri, P.; Brown, M.B.; Hutt, A.J.; Kirton, S.B. Revisiting the General Solubility Equation: In Silico Prediction of Aqueous Solubility Incorporating the Effect of Topographical Polar Surface Area. J. Chem. Inf. Model. 2012, 52, 420–428. [Google Scholar] [CrossRef]
  83. Irvine, J.D.; Takahashi, L.; Lockhart, K.; Cheong, J.; Tolan, J.W.; Selick, H.E.; Grove, J.R. MDCK (Madin-Darby Canine Kidney) Cells: A Tool for Membrane Permeability Screening. J. Pharm. Sci. 1999, 88, 28–33. [Google Scholar] [CrossRef]
  84. Yamashita, F.; Hashida, M. In Silico Approaches for Predicting ADME Properties of Drugs. Drug Metab Pharm. 2004, 19, 327–338. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Ma, X.L.; Chen, C.; Yang, J. Predictive Model of Blood-Brain Barrier Penetration of Organic Compounds. Acta Pharmacol. Sin. 2005, 26, 500–512. [Google Scholar] [CrossRef] [Green Version]
  86. Voutchkova, A.M.; Ferris, L.A.; Zimmerman, J.B.; Anastas, P.T. Toward Molecular Design for Hazard Reduction-Fundamental Relationships between Chemical Properties and Toxicity. Tetrahedron 2010, 66, 1031–1039. [Google Scholar] [CrossRef]
  87. Regina, K.J.; Brown, S.M.; Kharasch, E.D.; Crafford, A.; Holtzman, M.J.; Campbell, S.D. P-Glycoprotein Is a Major Determinant of Norbuprenorphine Brain Exposure and Antinociception. J. Pharmacol. Exp. Ther. 2012, 343, 53–61. [Google Scholar] [CrossRef] [Green Version]
  88. Hansch, C.; Silipo, C. Quantitative Structure-Activity Relation of Reversible Dihydrofolate Reductase Inhibitors. Diaminotriazines. J. Med. Chem. 1974, 17, 661–667. [Google Scholar] [CrossRef] [PubMed]
  89. Babatunde, S.; Adamu, O.; Gideon, U.; Sani, S. QSAR Modeling, Molecular Docking and ADMET/Pharmacokinetic Studies: A Chemometrics Approach to Search for Novel Inhibitors of Norepinephrine Transporter as Potent Antipsychotic Drugs. J. Iran. Chem. Soc. 2020, 17, 1953–1966. [Google Scholar] [CrossRef] [Green Version]
  90. Geldenhuys, W.J.; Bloomquist, J.R. Development of an a Priori Computational Approach for Brain Uptake of Compounds in an Insect Model System. Bioorg. Med. Chem. Lett. 2021, 40, 127930. [Google Scholar] [CrossRef]
  91. Szakács, G.; Váradi, A.; Özvegy-Laczka, C.; Sarkadi, B. The Role of ABC Transporters in Drug Absorption, Distribution, Metabolism, Excretion and Toxicity (ADME-Tox). Drug Discov. Today 2008, 13, 379–393. [Google Scholar] [CrossRef]
  92. Poongavanam, V.; Haider, N.; Ecker, G.F. Fingerprint-Based in Silico Models for the Prediction of P-Glycoprotein Substrates and Inhibitors. Bioorg. Med. Chem. 2012, 20, 5388–5395. [Google Scholar] [CrossRef] [Green Version]
  93. Montanari, F.; Ecker, G.F. Prediction of Drug-ABC-Transporter Interaction-Recent Advances and Future Challenges. Adv. Drug Deliv. Rev. 2015, 86, 17–26. [Google Scholar] [CrossRef] [Green Version]
  94. Figueira-Mansur, J.; Ferreira-Pereira, A.; Mansur, J.F.; Franco, T.A.; Alvarenga, E.S.L.; Sorgine, M.H.F.; Neves, B.C.; Melo, A.C.A.; Leal, W.S.; Masuda, H.; et al. Silencing of P-Glycoprotein Increases Mortality in Temephos-Treated Aedes Aegypti Larvae. Insect Mol. Biol. 2013, 22, 648–658. [Google Scholar] [CrossRef] [PubMed]
  95. Yavaşoglu, S.İ.; Alkhaibari, A.M.; Şimşek, F.M. Investigation of the insecticide resistance levels and mechanisms of the Mediterranean and Aegean Aedes caspius populations of Turkey. Int. J. Pest Manag. 2022. [Google Scholar] [CrossRef]
  96. Porretta, D.; Gargani, M.; Bellini, R.; Medici, A.; Punelli, F.; Urbanelli, S. Defence Mechanisms against Insecticides Temephos and Diflubenzuron in the Mosquito Aedes Caspius: The P-Glycoprotein Efflux Pumps. Med. Vet. Entomol. 2008, 22, 48–54. [Google Scholar] [CrossRef] [PubMed]
  97. Lima, E.P.; Goulart, M.O.F.; Neto, M.L.R. Evaluation of the Role of ATP-Binding Cassette Transporters as a Defense Mechanism against Temephos in Populations of Aedes Aegypti. Mem. Do Inst. Oswaldo Cruz 2014, 109, 964–966. [Google Scholar] [CrossRef] [Green Version]
  98. Merzendorfer, H. ABC Transporters and Their Role in Protecting Insects from Pesticides and Their Metabolites, 1st ed.; Elsevier Inc.: Amsterdam, The Netherlands, 2014; Volume 46, ISBN 9780124170100. [Google Scholar]
  99. Clark, R.D. Predicting Mammalian Metabolism and Toxicity of Pesticides in Silico. Pest Manag. Sci. 2018, 74, 1992–2003. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  100. Ji, C.; Svensson, F.; Zoufir, A.; Bender, A. EMolTox: Prediction of Molecular Toxicity with Confidence. Bioinformatics 2018, 34, 2508–2509. [Google Scholar] [CrossRef] [Green Version]
  101. Mokhnache, K.; Soltani, E.-K.; Madoui, S.; Khither, H.; Karbab, A.; Charef, N.; Arrar, L. In Silico Toxicological, Anti-Tubercular Effect Evaluation And In Vitro Marine Pathogenic Bacteria Inhibition of N-[(3-Chloro-4-Nitro-Phenyl)Methyleneamino]Pyridine-4-Carboxamidine. J. Drug Deliv. Ther. 2019, 9, 23–27. [Google Scholar] [CrossRef]
  102. Morger, A.; Mathea, M.; Achenbach, J.H.; Wolf, A.; Buesen, R.; Schleifer, K.J.; Landsiedel, R.; Volkamer, A. KnowTox: Pipeline and Case Study for Confident Prediction of Potential Toxic Effects of Compounds in Early Phases of Development. J. Cheminform. 2020, 12, 24. [Google Scholar] [CrossRef] [Green Version]
  103. Mokhnache, K.; Karbab, A.; Soltani, E.K.; Bououden, W.; Ouhida, S.; Arrar, L.; Esteban, M.A.; Charef, N.; Mubarak, M.S. Synthesis, Characterization, Toxic Substructure Prediction, Hepatotoxicity Evaluation, Marine Pathogenic Bacteria Inhibition, and DFT Calculations of a New Hydrazone Derived from Isoniazid. J. Mol. Struct. 2020, 1221. [Google Scholar] [CrossRef]
  104. Abaukaka, Y.A.; Sanusi, S.; Ozigi, K.A.; Malo, F.U. Assessment of the Cytotoxic and Mutagenic Potential of Dichlorvos (DDVP) Using in Silico Classification Model; a Health Hazard Awareness in Nigeria. Environ. Health Toxicol. 2020, 35, e2020016. [Google Scholar] [CrossRef]
  105. Sosnin, S.; Karlov, D.; Tetko, I.V.; Fedorov, M.V. Comparative Study of Multitask Toxicity Modeling on a Broad Chemical Space. J. Chem. Inf. Model. 2019, 59, 1062–1072. [Google Scholar] [CrossRef] [Green Version]
  106. Hernandez, H.; Shivraj, L. In Silico Toxicity Prediction Using an Integrative Multimodel Approach. Res. Rep. 2020, 5, 2020–2022. [Google Scholar] [CrossRef]
  107. Tugcu, G.; Kırmızıbekmez, H.; Aydın, A. The Integrated Use of in Silico Methods for the Hepatotoxicity Potential of Piper Methysticum. Food Chem. Toxicol. 2020, 145, 111663. [Google Scholar] [CrossRef] [PubMed]
  108. Bhamare, P.; Umadoss, P.; Upmanyu, N.; Dubey, R. Identification, Isolation, Structural Characterisation, Synthesis and: In Silico Toxicity Prediction of the Alkaline Hydrolytic Degradation Product of Brivaracetam by Using LC-PDA, Preparative HPLC, LC/HESI/LTQ, FTIR, and 1H NMR. Anal. Methods 2020, 12, 1868–1881. [Google Scholar] [CrossRef]
  109. Rashid, M.I.; Andleeb, S.; Ali, A. Evaluation of Pyocyanin Induced Systemic Pathogenicity of Pseudomonas Aeruginosa. Pak. J. Pharm. Sci. 2020, 33, 915–922. [Google Scholar] [CrossRef]
  110. Webel, H.E.; Kimber, T.B.; Radetzki, S.; Neuenschwander, M.; Nazaré, M.; Volkamer, A. Revealing Cytotoxic Substructures in Molecules Using Deep Learning. J. Comput.-Aided Mol. Des. 2020, 34, 731–746. [Google Scholar] [CrossRef] [Green Version]
  111. Kim, R.B. Organic Anion-Transporting Polypeptide (OATP) Transporter Family and Drug Disposition. Eur. J. Clin. Investig. Suppl. 2003, 33, 1–5. [Google Scholar] [CrossRef]
  112. Marzolini, C.; Tirona, R.G.; Kim, R.B. Pharmacogenomics of the OATP and OAT Families. Pharmacogenomics 2004, 5, 273–282. [Google Scholar] [CrossRef]
  113. Smith, N.F.; Figg, W.D.; Sparreboom, A. Role of the Liver-Specific Transporters OATP1B1 and OATP1B3 in Governing Drug Elimination. Expert Opin. Drug Metab. Toxicol. 2005, 1, 429–445. [Google Scholar] [CrossRef]
  114. Alam, K.; Pahwa, S.; Wang, X.; Zhang, P.; Ding, K.; Alaa, H.; Li, L.; Yue, W.; City, O.; States, U.; et al. Downregulation of Organic Anion Transporting Polypeptide (OATP) 1B1 Transport Function by Lysosomotropic Drug Chloroquine: Implication in OATP-Mediated Drug-Drug Interactions. Mol. Pharmaceut. 2016, 13, 839–851. [Google Scholar] [CrossRef] [Green Version]
  115. Van De Steeg, E.; Stránecký, V.; Hartmannová, H.; Nosková, L.; Hřebíček, M.; Wagenaar, E.; Van Esch, A.; De Waart, D.R.; Oude Elferink, R.P.J.; Kenworthy, K.E.; et al. Complete OATP1B1 and OATP1B3 Deficiency Causes Human Rotor Syndrome by Interrupting Conjugated Bilirubin Reuptake into the Liver. J. Clin. Investig. 2012, 122, 519–528. [Google Scholar] [CrossRef] [PubMed]
  116. Vankayala, S.L.; Hargis, J.C.; Woodcock, H.L. How Does Catalase Release Nitric Oxide? A Computational Structure-Activity Relationship Study. J. Chem. Inf. Model. 2013, 53, 2951–2961. [Google Scholar] [CrossRef] [Green Version]
  117. Berrino, E.; Bozdag, M.; Del Prete, S.; Alasmary, F.A.S.; Alqahtani, L.S.; AlOthman, Z.; Capasso, C.; Supuran, C.T. Inhibition of α-, β-, γ-, and δ-Carbonic Anhydrases from Bacteria and Diatoms with N′-Aryl-N-Hydroxy-Ureas. J. Enzym. Inhib. Med. Chem. 2018, 33, 1194–1198. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  118. Fahmy, M.A.; Abdalla, E.F. Genotoxicity Evaluation of Buprofezin, Petroleum Oil and Profenofos in Somatic and Germ Cells of Male Mice. J. Appl. Toxicol. 1998, 18, 301–305. [Google Scholar] [CrossRef]
  119. Booth, D.F.M. Pesticides in Aquatic Environments. In Environmental Science Research; Springer Science & Business Media: Berlin/Heidelberg, Germany, 1977; p. 263. ISBN 978-1-4684-2870-4. [Google Scholar]
  120. Sun, R.; Liu, C.; Zhang, H.; Wang, Q. Benzoylurea Chitin Synthesis Inhibitors. J. Agric. Food Chem. 2015, 63, 6847–6865. [Google Scholar] [CrossRef] [PubMed]
  121. Chang, J.; Li, W.; Xu, P.; Guo, B.; Wang, Y.; Li, J.; Wang, H. The Tissue Distribution, Metabolism and Hepatotoxicity of Benzoylurea Pesticides in Male Eremias Argus after a Single Oral Administration. Chemosphere 2017, 183, 1–8. [Google Scholar] [CrossRef]
  122. Greene, N.; Fisk, L.; Naven, R.T.; Note, R.R.; Patel, M.L.; Pelletier, D.J. Developing Structure-Activity Relationships for the Prediction of Hepatotoxicity. Chem. Res. Toxicol. 2010, 23, 1215–1222. [Google Scholar] [CrossRef]
  123. Banerjee, P.; Dehnbostel, F.O.; Preissner, R. Prediction Is a Balancing Act: Importance of Sampling Methods to Balance Sensitivity and Specificity of Predictive Models Based on Imbalanced Chemical Data Sets. Front. Chem. 2018, 6, 362. [Google Scholar] [CrossRef]
  124. de Knecht, J.A.; van Herwijnen, R. Environmental Risk Limits for Deltamethrin; RIVM Letter report 601716015/2008; National Institute for Public Health and the Environment: Bilthoven, The Netherlands, 2008; pp. 1–58. [Google Scholar]
  125. Olsvik, P.A.; Samuelsen, O.B.; Agnalt, A.L.; Lunestad, B.T. Transcriptional Responses to Teflubenzuron Exposure in European Lobster (Homarus Gammarus). Aquat. Toxicol. 2015, 167, 143–156. [Google Scholar] [CrossRef]
  126. Samuelsen, O.B.; Lunestad, B.T.; Farestveit, E.; Grefsrud, E.S.; Hannisdal, R.; Holmelid, B.; Tjensvoll, T.; Agnalt, A.L. Mortality and Deformities in European Lobster (Homarus Gammarus) Juveniles Exposed to the Anti-Parasitic Drug Teflubenzuron. Aquat. Toxicol. 2014, 149, 8–15. [Google Scholar] [CrossRef]
  127. Vila-Carriles, W.H.; Zhao, G.; Bryan, J. Defining a Binding Pocket for Sulfonylureas in ATP-sensitive Potassium Channels. FASEB J. 2007, 21, 18–25. [Google Scholar] [CrossRef] [PubMed]
  128. Retnakaran, A. Effect of 3 New Moult-Inhibiting Insect Growth Regulators on the Spruce Budworm1. J. Econ. Entomol. 1980, 73, 520–524. [Google Scholar] [CrossRef]
  129. Mandloi, M.; Karmarkar, S.; Joshi, S.; Khadikar, P. V Inhibition of DNA Synthesis by Hydroxyureas and Related Compounds: Structure-Activity Relationships. Oxid. Commun. 2001, 24, 291–301. [Google Scholar]
  130. Das, P.K.; Caaveiro, J.M.M.; Luque, S.; Klibanov, A.M. Binding of Hydrophobic Hydroxamic Acids Enhances Peroxidase’s Stereoselectivity in Nonaqueous Sulfoxidations. J. Am. Chem. Soc. 2002, 124, 782–787. [Google Scholar] [CrossRef]
  131. Uesato, S.; Hashimoto, Y.; Nishino, M.; Nagaoka, Y.; Kuwajima, H. N-Substituted Hydroxyureas as Urease Inhibitors. Chem. Pharm. Bull. 2002, 50, 1280–1282. [Google Scholar] [CrossRef] [Green Version]
  132. Rohrman, B.A.; Mazziotti, D.A. Quantum Chemical Design of Hydroxyurea Derivatives for the Treatment of Sickle-Cell Anemia. J. Phys. Chem. B 2005, 109, 13392–13396. [Google Scholar] [CrossRef]
  133. Ryan, S.W. (12) United States Patent Date of Patent. Syst. Method Program. A Weigh. Scale Usinga Key Signal Enter A Program. Mode 2009, 1, 14. [Google Scholar]
  134. McCourt, M.P.; Greene, R.S.; Smith, C.; McClory, H.; Wirth, C.; Sumbler, B.; Emery, M.L.; Leskiw, B.D.; Krause, J.G.; McGahan, M.E.; et al. Computational Analysis of the Apoptotic Activity of a Small Library of Novel Substituted Hydroxy Ureas Designed as Potential Inhibitors of Ribonucleotide Reductase. Lett. Drug Des. Discov. 2011, 8, 580–584. [Google Scholar] [CrossRef]
  135. Elford, H.L.; Van’t Riet, B.; Wampler, G.L.; Lin, A.L.; Elford, R.M. Regulation of Ribonucleotide Reductase in Mammalian Cells by Chemotherapeutic Agents. Adv. Enzym. Regul. 1981, 19, 151–168. [Google Scholar] [CrossRef]
  136. Grady, R.W.; Bienen, E.J.; Clarkson, A.B. P-Alkyloxybenz Hydroxamic Acids, Effective Inhibitors of the trypanosome Glycerol-3-Phosphate Oxidase. Mol. Biochem. Parasitol. 1986, 19, 231–240. [Google Scholar] [CrossRef]
  137. Gupta, S.P.; Gupta, J.K. Quantitative Structure-Activity Relationship Study on Some 5-Lipoxygenase Inhibitors. J. Enzym. Inhib. Med. Chem. 1990, 3, 179–188. [Google Scholar] [CrossRef] [PubMed]
  138. Summers, J.B.; Mazdiyasni, H.; Holms, J.H.; Ratajczyk, J.D.; Stewart, A.O.; Dyer, R.D.; Carter, G.W.; Kim, K.H. Hydroxamic Acid Inhibitors of 5-Lipoxygenase: Quantitative Structure–Activity Relationships. J. Med. Chem. 1990, 33, 992–998. [Google Scholar] [CrossRef]
  139. Khadika, P.V.; Sharma, S.; Sharma, V.; Joshi, S.; Lukovits, I.; Kaveeshwar, M. A QSAR Study of the Effect of Benzohydroxamic Acids on DNA Synthesis. Bull. Des Sociétés Chim. Belg. 1997, 106, 767–772. [Google Scholar]
  140. Xuan, X.; Yifan, L. The Pattern Recognition Study on the Antitumor Activity of Aroyl Compounds. Chin. J. Med. Chem. 2001, 11, 17–20. [Google Scholar]
  141. Liénard, B.M.R.; Horsfall, L.E.; Galleni, M.; Frère, J.M.; Schofield, C.J. Inhibitors of the FEZ-1 Metallo-β-Lactamase. Bioorg. Med. Chem. Lett. 2007, 17, 964–968. [Google Scholar] [CrossRef] [PubMed]
  142. Bayat, A.; Fattahi, A. Investigation of the Scavenging Mechanism of Tyrosyl Radical by Hydroxybenzohydroxamic Acid Derivatives: A DFT Study. Comput. Theor. Chem. 2013, 1018, 35–44. [Google Scholar] [CrossRef]
  143. Lim, J.; Song, Y.; Jang, J.H.; Jeong, C.H.; Lee, S.; Park, B.; Seo, Y.H. Aspirin-Inspired Acetyl-Donating HDACs Inhibitors. Arch. Pharmacal Res. 2018, 41, 967–976. [Google Scholar] [CrossRef]
  144. Hania, M. Synthesis and antibacterial activity of 2-acetyl and 2-formylpyrrole oximes, semicarbazones and phenylhydrazones. Asian J. Chem. 2002, 14, 2. [Google Scholar]
  145. Ramos, R.S.; Borges, R.S.; de Souza, J.S.N.; Araujo, I.F.; Chaves, M.H.; Santos, C.B.R. Identification of Potential Antiviral Inhibitors from Hydroxychloroquine and 1,2,4,5-Tetraoxanes Analogues and Investigation of the Mechanism of Action in SARS-CoV-2. Int. J. Mol. Sci. 2022, 23, 1781. [Google Scholar] [CrossRef]
  146. Fukunishi, Y.; Kurosawa, T.; Mikami, Y.; Nakamura, H. Prediction of Synthetic Accessibility Based on Commercially Available Compound Databases. J. Chem. Inf. Model. 2014, 54, 3259–3267. [Google Scholar] [CrossRef]
  147. Gentile, F.; Yaacoub, J.C.; Gleave, J.; Fernandez, M.; Ton, A.-T.; Ban, F.; Stern, A.; Cherkasov, A. Artificial Intelligence–Enabled Virtual Screening of Ultra-Large Chemical Libraries with Deep Docking. Nat. Protoc. 2022, 17, 672–697. [Google Scholar] [CrossRef] [PubMed]
  148. Buettelmann, B.; Jakob-Roetne, R.; Knust, H.; Lucas, M.C.; Thomas, A. Preparation of Isoxazole-Pyridine Derivatives as GABA Modulators. PCT Int. Appl. WO 2010, 4, 65. [Google Scholar]
  149. Mo, X.; Morgan, T.D.R.; Ang, H.T.; Hall, D.G. Scope and Mechanism of a True Organocatalytic Beckmann Rearrangement with a Boronic Acid/Perfluoropinacol System under Ambient Conditions. J. Am. Chem. Soc. 2018, 140, 5264–5271. [Google Scholar] [CrossRef]
  150. Zeng, G.; Li, H.; Wei, Y.; Xuan, W.; Zhang, R.; Breden, L.E.; Wang, W.; Liang, F. Sen Engineering Iron Responses in Mammalian Cells by Signal-Induced Protein Proximity. ACS Synth. Biol. 2017, 6, 921–927. [Google Scholar] [CrossRef] [PubMed]
  151. Berini, C.; Minassian, F.; Pelloux-Léon, N.; Denis, J.N.; Vallée, Y.; Philouze, C. Efficient Stereoselective Nucleophilic Addition of Pyrroles to Chiral Nitrones. Org. Biomol. Chem. 2008, 6, 2574–2586. [Google Scholar] [CrossRef]
  152. Niu, C.; Yin, L.; Nie, L.F.; Dou, J.; Zhao, J.Y.; Li, G.; Aisa, H.A. Synthesis and Bioactivity of Novel Isoxazole Chalcone Derivatives on Tyrosinase and Melanin Synthesis in Murine B16 Cells for the Treatment of Vitiligo. Bioorg. Med. Chem. 2016, 24, 5440–5448. [Google Scholar] [CrossRef]
  153. Yu, J.; Cao, X.; Lu, M. A Novel and Efficient Catalytic System Including TEMPO/Acetaldoxime/InCl3 for Aerobic Oxidation of Primary Amines to Oximes. Tetrahedron Lett. 2014, 55, 5751–5755. [Google Scholar] [CrossRef]
  154. Viswanathan, V.S.; Ryan, M.J.; Dhruv, H.D.; Gill, S.; Eichhoff, O.M.; Seashore-Ludlow, B.; Kaffenberger, S.D.; Eaton, J.K.; Shimada, K.; Aguirre, A.J.; et al. Dependency of a therapy-resistant state of cancer cells on a lipid peroxidase pathway. Nature 2017, 547, 453–457. [Google Scholar] [CrossRef]
  155. Kim, J.H.; Ali, K.H.; Oh, Y.J.; Seo, Y.H. Design, synthesis, and biological evaluation of histone deacetylase inhibitor with novel salicylamide zinc binding group. Medicine 2022, 101, e29049. [Google Scholar] [CrossRef]
  156. Orlowska, E.; Roller, A.; Wiesinger, H.; Pignitter, M.; Jirsa, F.; Krachler, R.; Kandioller, W.; Keppler, B.K. Benzoic Hydroxamate-Based Iron Complexes as Model Compounds for Humic Substances: Synthesis, Characterization and Algal Growth Experiments. RSC Adv. 2016, 6, 40238–40249. [Google Scholar] [CrossRef] [Green Version]
  157. Redina, E.A.; Vikanova, K.V.; Kapustin, G.I.; Mishin, I.V.; Tkachenko, O.P.; Kustov, L.M. Selective Room-Temperature Hydrogenation of Carbonyl Compounds under Atmospheric Pressure over Platinum Nanoparticles Supported on Ceria-Zirconia Mixed Oxide. Eur. J. Org. Chem. 2019, 2019, 4159–4170. [Google Scholar] [CrossRef]
  158. Mayer, R.T.; Chen, A.C.; Deloach, J.R.; Mayer, R.T.; Chen, A.C.; Deloach, J.R. Characterization of A Chitin Synthase from the Stable Fly, Stomoxys Calcitrans (L.). Insect Biochem. 1980, 10, 549–556. [Google Scholar] [CrossRef]
  159. Clark, M.; Cramer, R.D.; Van Opdenbosch, N. Validation of the General Purpose Tripos 5.2 Force Field. J. Comput. Chem. 1989, 10, 982–1012. [Google Scholar] [CrossRef]
  160. Gasteiger, J.; Marsili, M. Iterative Partial Equalization of Orbital Electronegativity-a Rapid Access to Atomic Charges. Tetrahedron 1980, 36, 3219–3228. [Google Scholar] [CrossRef]
  161. Richmond, N.J.; Abrams, C.A.; Wolohan, P.R.N.; Abrahamian, E.; Willett, P.; Clark, R.D. GALAHAD: 1. Pharmacophore Identification by Hypermolecular Alignment of Ligands in 3D. J. Comput.-Aided Mol. Des. 2006, 20, 567–587. [Google Scholar] [CrossRef]
  162. Mysinger, M.M.; Carchia, M.; Irwin, J.J.; Shoichet, B.K. Directory of Useful Decoys, Enhanced (DUD-E): Better Ligands and Decoys for Better Benchmarking. J. Med. Chem. 2012, 55, 6582–6594. [Google Scholar] [CrossRef] [PubMed]
  163. Lätti, S.; Niinivehmas, S.; Pentikäinen, O.T. Rocker: Open Source, Easy-to-Use Tool for AUC and Enrichment Calculations and ROC Visualization. J. Cheminform. 2016, 8, 45. [Google Scholar] [CrossRef] [Green Version]
  164. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. BLAST+: Architecture and Applications. BMC Bioinform. 2009, 9, 421. [Google Scholar] [CrossRef] [Green Version]
  165. Guex, N.; Peitsch, M.C.; Schwede, T. Celebrating 30 Years Automated Comparative Protein Structure Modeling with SWISS-MODEL and Swiss-PdbViewer: A Historical Perspective. Electrophoresis 2009, 30, 162–173. [Google Scholar] [CrossRef]
  166. Waterhouse, A.; Bertoni, M.; Bienert, S.; Studer, G.; Tauriello, G.; Gumienny, R.; Heer, F.T.; De Beer, T.A.P.; Rempfer, C.; Bordoli, L.; et al. SWISS-MODEL: Homology Modelling of Protein Structures and Complexes. Nucleic Acids Res. 2018, 46, 296–303. [Google Scholar] [CrossRef] [Green Version]
  167. Bienert, S.; Waterhouse, A.; De Beer, T.A.P.; Tauriello, G.; Studer, G.; Bordoli, L.; Schwede, T. The SWISS-MODEL Repository—New Features and Functionality. Nucleic Acids Res. 2017, 45, 313–319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  168. Studer, G.; Rempfer, C.; Waterhouse, A.M.; Haas, J.; Schwede, T. QMEANDisCo–Distance Constraints Applied on Model Quality Estimation. Struct. Bioinform. 2019, 36, 1765–1771. [Google Scholar] [CrossRef] [PubMed]
  169. Benkert, P.; Biasini, M.; Schwede, T. Toward the Estimation of the Absolute Quality of Individual Protein Structure Models. Struct. Bioinform. 2011, 27, 343–350. [Google Scholar] [CrossRef]
  170. Wang, G.; Dunbrack, R.L., Jr. Bioinformatics Applications note pisces: A Protein Sequence Culling Server. Bioinformatics 2003, 19, 1589–1591. [Google Scholar] [CrossRef] [Green Version]
  171. Cole, J.C.; Murray, C.W.; Nissink, J.W.M.; Taylor, R.D.; Taylor, R. Comparing Protein-Ligand Docking Programs Is Difficult. Proteins: Struct. Funct. Genet. 2005, 60, 325–332. [Google Scholar] [CrossRef]
  172. Ramírez, D.; Caballero, J. Is It Reliable to Take the Molecular Docking Top Scoring Position as the Best Solution without Considering Available Structural Data? Molecules 2018, 23, 1038. [Google Scholar] [CrossRef] [Green Version]
  173. Ramos, R.d.S.; Costa, J.d.S.; Macêdo, W.J.d.C.; Rodrigues, A.B.L.; Silva, C.; Santos, C.; da Costa, G.V.; Taft, C.A.; Silva, R.C.; Souto, R.N.P.; et al. Identification of Potential Inhibitors from Pyriproxyfen with Insecticidal Activity by Virtual Screening. Pharmaceuticals 2019, 12, 20. [Google Scholar] [CrossRef] [Green Version]
  174. Leão, R.P.; Cruz, J.V.; da Costa, G.V.; Cruz, J.N.; Ferreira, E.F.B.; Silva, R.C.; de Lima, L.R.; Borges, R.S.; Dos Santos, G.B.; Santos, C.B.R. Identification of New Rofecoxib-Based Cyclooxygenase-2 Inhibitors: A Bioinformatics Approach. Pharmaceuticals 2020, 13, 209. [Google Scholar] [CrossRef]
  175. Costa, R.A.; Rocha, E.C.M.; Silva, R.C.; Gonçalves, A.S.; Santos, C.B.R.; Brasil, D.S.B. A Computational Approach Applied to the Study of Potential Allosteric Inhibitors Protease NS2B/NS3 from Dengue Virus. Molecules 2022, 27, 20. [Google Scholar] [CrossRef]
  176. Boda, K.; Seidel, T.; Gasteiger, J. Structure and Reaction Based Evaluation of Synthetic Accessibility. J. Comput.-Aided Mol. Des. 2007, 21, 311–325. [Google Scholar] [CrossRef]
  177. Kochev, N.T.; Paskaleva, V.H.; Jeliazkova, N. Ambit-Tautomer: An Open Source Tool for Tautomer Generation. Mol. Inform. 2013, 32, 481–504. [Google Scholar] [CrossRef] [PubMed]
  178. Ertl, P.; Schuffenhauer, A. Estimation of Synthetic Accessibility Score of Drug-like Molecules Based on Molecular Complexity and Fragment Contributions. J. Cheminform. 2009, 1, 8–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  179. Sun, J.; Jeliazkova, N.; Chupakin, V.; Golib-Dzib, J.F.; Engkvist, O.; Carlsson, L.; Wegner, J.; Ceulemans, H.; Georgiev, I.; Jeliazkov, V.; et al. ExCAPE-DB: An Integrated Large Scale Dataset Facilitating Big Data Analysis in Chemogenomics. J. Cheminform. 2017, 9, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Scheme 1. Methodological scheme. First stage (pharmacophore models) > second stage (comparative modeling) and third stage (hierarchical virtual screening—model and coupling).
Scheme 1. Methodological scheme. First stage (pharmacophore models) > second stage (comparative modeling) and third stage (hierarchical virtual screening—model and coupling).
Ijms 23 08218 sch001
Figure 1. ROC curve for pharmacophore models. Below the diagonal line, poor models are represented for random selection (AUC < 0.5).
Figure 1. ROC curve for pharmacophore models. Below the diagonal line, poor models are represented for random selection (AUC < 0.5).
Ijms 23 08218 g001
Figure 2. (A)—potent inhibitor (0.001 mg·L−1; QFIT = 84.94); (B)—flucycloxuron (0.05 mg·L−1; QFIT = 58.66), superposed to model pharmacophore 06. Spheres represent H-bond donors (green), H-bond acceptors (magenta) and hydrophobic groups (cyan). The size of spheres varies according to the tolerance radius calculated using GALAHADTM. All the distances are measured in angstroms.
Figure 2. (A)—potent inhibitor (0.001 mg·L−1; QFIT = 84.94); (B)—flucycloxuron (0.05 mg·L−1; QFIT = 58.66), superposed to model pharmacophore 06. Spheres represent H-bond donors (green), H-bond acceptors (magenta) and hydrophobic groups (cyan). The size of spheres varies according to the tolerance radius calculated using GALAHADTM. All the distances are measured in angstroms.
Ijms 23 08218 g002
Figure 3. Overlap between alpha-carbons of macromolecule 5ZNT (green) and the Bombyx-Aedes model (red); RMSD value 0.396.
Figure 3. Overlap between alpha-carbons of macromolecule 5ZNT (green) and the Bombyx-Aedes model (red); RMSD value 0.396.
Ijms 23 08218 g003
Figure 4. Ramachandran, the Bombyx-Aedes Model terrain. The most favorable regions are shown in dark green, neutral in light green, and non-favorable in white; each amino acid residue is marked with different colors, generated by the structure assessment program (https://swissmodel.expasy.org/assess/KGJVtX/01, accessed on 27 June 2021).
Figure 4. Ramachandran, the Bombyx-Aedes Model terrain. The most favorable regions are shown in dark green, neutral in light green, and non-favorable in white; each amino acid residue is marked with different colors, generated by the structure assessment program (https://swissmodel.expasy.org/assess/KGJVtX/01, accessed on 27 June 2021).
Ijms 23 08218 g004
Figure 5. Overlays of crystallographic ligand poses in JH3 (in red) with the calculated pose (in blue).
Figure 5. Overlays of crystallographic ligand poses in JH3 (in red) with the calculated pose (in blue).
Ijms 23 08218 g005
Figure 6. Interaction of selected molecules, in yellow, with the amino acid residues, in green, of the macromolecule 5V13 (juvenile hormone), in green. (A) M01 (ZINC13514543); (B) M02 (ZINC14413017); (C) M03 (ZINC05680928); (D) M04 (ZINC04691148); (E) M05 (ZINC95851150); (F) JH3.
Figure 6. Interaction of selected molecules, in yellow, with the amino acid residues, in green, of the macromolecule 5V13 (juvenile hormone), in green. (A) M01 (ZINC13514543); (B) M02 (ZINC14413017); (C) M03 (ZINC05680928); (D) M04 (ZINC04691148); (E) M05 (ZINC95851150); (F) JH3.
Ijms 23 08218 g006
Figure 7. Interactions in the catalytic region around the metal zinc in the 5ZNT protein (A) and the Bombyx-Aedes model (B) and the overlapping of the amino acid residues closest to the zinc metal in the 5ZNT protein (C) and the Bombyx-Aedes model (D).
Figure 7. Interactions in the catalytic region around the metal zinc in the 5ZNT protein (A) and the Bombyx-Aedes model (B) and the overlapping of the amino acid residues closest to the zinc metal in the 5ZNT protein (C) and the Bombyx-Aedes model (D).
Ijms 23 08218 g007
Figure 8. Interaction between the selected molecules, in yellow, with the amino acid residues, in green, of the macromolecule Bombyx-Aedes model, in green. (A) M01 (ZINC13514543); (B) M02 (ZINC14413017); (C) M03 (ZINC05680928); (D) M04 (ZINC04691148); (E) M05 (ZINC95851150); (F) Zinc Ray.
Figure 8. Interaction between the selected molecules, in yellow, with the amino acid residues, in green, of the macromolecule Bombyx-Aedes model, in green. (A) M01 (ZINC13514543); (B) M02 (ZINC14413017); (C) M03 (ZINC05680928); (D) M04 (ZINC04691148); (E) M05 (ZINC95851150); (F) Zinc Ray.
Ijms 23 08218 g008
Figure 9. Promising compounds (A) M01 (ZINC13514543); (B) M02 (ZINC14413017); (C) M03 (ZINC05680928); (D) M04 (ZINC04691148); (E) M05 (ZINC95851150). Source: authors.
Figure 9. Promising compounds (A) M01 (ZINC13514543); (B) M02 (ZINC14413017); (C) M03 (ZINC05680928); (D) M04 (ZINC04691148); (E) M05 (ZINC95851150). Source: authors.
Ijms 23 08218 g009
Figure 10. The basic chemical structure of benzoylphenylurean compounds consists of three regions for SAR analysis: the benzoyl ring (A), the urea bridge (B), and the aniline (C).
Figure 10. The basic chemical structure of benzoylphenylurean compounds consists of three regions for SAR analysis: the benzoyl ring (A), the urea bridge (B), and the aniline (C).
Ijms 23 08218 g010
Figure 11. Synthetic route for the preparation of M01. DCM (dichloromethane). Starting material I is commercially available.
Figure 11. Synthetic route for the preparation of M01. DCM (dichloromethane). Starting material I is commercially available.
Ijms 23 08218 g011
Figure 12. Theoretical synthetic route for the preparation of M02. Starting material IV is commercially available.
Figure 12. Theoretical synthetic route for the preparation of M02. Starting material IV is commercially available.
Ijms 23 08218 g012
Figure 13. Alternative synthetic route for the preparation of M02. Starting material VI is commercially available.
Figure 13. Alternative synthetic route for the preparation of M02. Starting material VI is commercially available.
Ijms 23 08218 g013
Figure 14. Traditional method for the preparation of compound M03. Starting material VIII is commercially available.
Figure 14. Traditional method for the preparation of compound M03. Starting material VIII is commercially available.
Ijms 23 08218 g014
Figure 15. Alternative synthetic route for the preparation of compound M03. Starting material IX is commercially available.
Figure 15. Alternative synthetic route for the preparation of compound M03. Starting material IX is commercially available.
Ijms 23 08218 g015
Figure 16. Synthetic route for the preparation of MO4. THF (tetrahydrofuran), Potassium hydroxide, Starting material X is commercially available.
Figure 16. Synthetic route for the preparation of MO4. THF (tetrahydrofuran), Potassium hydroxide, Starting material X is commercially available.
Ijms 23 08218 g016
Figure 17. Theoretical synthetic route for the preparation of M05. Starting materials XII and XV are commercially available.
Figure 17. Theoretical synthetic route for the preparation of M05. Starting materials XII and XV are commercially available.
Ijms 23 08218 g017
Figure 18. The 2D structure of benzoylphenylurean derivatives. R1 and R2 are radicals, and X = methyl or halogens.
Figure 18. The 2D structure of benzoylphenylurean derivatives. R1 and R2 are radicals, and X = methyl or halogens.
Ijms 23 08218 g018
Figure 19. Structures of selected training set compounds, with canonical smiles and concentrations ([ ]) with 100% larvicidal activity (mg·L−1). (A)—Compound 01; (B)—Compound 05; (C)—Compound 06; (D)—Compound 07; (E)—Compound 12; (F)—Compound 13; (G)—Compound 16.
Figure 19. Structures of selected training set compounds, with canonical smiles and concentrations ([ ]) with 100% larvicidal activity (mg·L−1). (A)—Compound 01; (B)—Compound 05; (C)—Compound 06; (D)—Compound 07; (E)—Compound 12; (F)—Compound 13; (G)—Compound 16.
Ijms 23 08218 g019
Table 1. Internal statistical parameters for pharmacophore models provided by GALAHAD™.
Table 1. Internal statistical parameters for pharmacophore models provided by GALAHAD™.
ModelENERGY
(kcal/mol)
STERICSH-BONDMOL_QRY
015.552662.60627.30365.56
026.832502.50669.00378.08
036.622724.20646.60358.13
046.882888.40619.20356.44
056.352496.40650.00362.48
066.732876.90610.40350.63
077.872936.50641.60357.15
086.312604.20647.90343.24
097.032638.60674.00335.78
106.822802.00611.70346.41
Table 2. BEDROC value of pharmacophore models, considering α = 16.10.
Table 2. BEDROC value of pharmacophore models, considering α = 16.10.
ModelsBEDROC
010.36
020.30
040.41
050.23
060.60
070.49
090.52
100.44
Table 3. Pharmacokinetic properties of selected molecules.
Table 3. Pharmacokinetic properties of selected molecules.
In Silico Pharmacokinetic Properties
CODSwissADMEPREADMET
LipophilicityWater Solubility
Consensus log Po/w 1 X Log SBBB 2Caco-2 3MDCK 4
DFB3.32−5.062.2520.470.53
TEF4.43−5.944.6623.090.07
FCX5.37−8.232.2925.120.05
BPU3.35−4.570.1356.972.61
M010.55−1.460.5320.77247.45
M020.90−1.551.7013.7031.05
M031.46−2.130.613.20395.90
M042.42−3.431.6119.290.46
M050.28−1.040.3519.063.21
1 Log Po/w, partition coefficients for octanol/water); 2 BBB, in silico prediction of brain/blood partition coefficient; 3 Caco-2, in silico prediction of colon adenocarcinoma cell; 4 MDCK, in silico prediction of a Madin-Darby canine kidney cell.
Table 4. Toxicological properties of selected molecules in the toxicity model—eMolTox.
Table 4. Toxicological properties of selected molecules in the toxicity model—eMolTox.
CODActionInjuryConfidence
FCXAgonist of liver X receptor betaLiver0.993
Antagonist of the farnesoid-X receptor (FXR) signaling pathway0.995
Disruptors of the mitochondrial membrane potential0.981
Block Bile Salt Export Pump0.991
Differential cytotoxicity against isogenic chicken DT40 cell lines with known DNA damage response pathways—Rad54/Ku70 mutant cell lineGenotoxicity0.992
Differential cytotoxicity (isogenic chicken DT40 cell lines)Genotoxicity0.986
DFBActivator of the aryl hydrocarbon receptor (AhR) signaling pathwayLiver0.997
Modulator of P2X purinoceptor 7Nervous system, immune, kidney0.985
TEFModulator of P2X purinoceptor 7Nervous system, immune, kidney0.993
BUPDifferential cytotoxicity against isogenic chicken DT40 cell lines with known DNA damage response pathways—Rad54/Ku70 mutant cell lineGenotoxicity0.992
Antagonist of the estrogen receptor alpha (ER-alpha) signaling pathwayEndocrine0.996
Activators of the human pregnane X receptor (PXR) signaling pathwayLiver0.996
Agonist of the constitutive androstane receptor (CAR) signaling pathway0.992
M01We do not find any toxic action from data-driven models.
M02Block OATP1B1 TransporterLiver0.993
M03Block OATP1B1 TransporterLiver0.982
M04Skin SensitizationSkin0.999
Block OATP1B1 TransporterLiver0.999
M05We do not find any toxic action from data-driven models.
Table 5. Prediction of synthetic accessibility (SA).
Table 5. Prediction of synthetic accessibility (SA).
CODCOMPOUNDSSA
M01ZINC1351454396.096
M02ZINC1441301792.731
M03ZINC0568092889.776
M04ZINC0469114892.951
M05ZINC9585115091.439
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Da Costa, G.V.; Neto, M.F.A.; Da Silva, A.K.P.; De Sá, E.M.F.; Cancela, L.C.F.; Vega, J.S.; Lobato, C.M.; Zuliani, J.P.; Espejo-Román, J.M.; Campos, J.M.; et al. Identification of Potential Insect Growth Inhibitor against Aedes aegypti: A Bioinformatics Approach. Int. J. Mol. Sci. 2022, 23, 8218. https://doi.org/10.3390/ijms23158218

AMA Style

Da Costa GV, Neto MFA, Da Silva AKP, De Sá EMF, Cancela LCF, Vega JS, Lobato CM, Zuliani JP, Espejo-Román JM, Campos JM, et al. Identification of Potential Insect Growth Inhibitor against Aedes aegypti: A Bioinformatics Approach. International Journal of Molecular Sciences. 2022; 23(15):8218. https://doi.org/10.3390/ijms23158218

Chicago/Turabian Style

Da Costa, Glauber V., Moysés F. A. Neto, Alicia K. P. Da Silva, Ester M. F. De Sá, Luanne C. F. Cancela, Jeanina S. Vega, Cássio M. Lobato, Juliana P. Zuliani, José M. Espejo-Román, Joaquín M. Campos, and et al. 2022. "Identification of Potential Insect Growth Inhibitor against Aedes aegypti: A Bioinformatics Approach" International Journal of Molecular Sciences 23, no. 15: 8218. https://doi.org/10.3390/ijms23158218

APA Style

Da Costa, G. V., Neto, M. F. A., Da Silva, A. K. P., De Sá, E. M. F., Cancela, L. C. F., Vega, J. S., Lobato, C. M., Zuliani, J. P., Espejo-Román, J. M., Campos, J. M., Leite, F. H. A., & Santos, C. B. R. (2022). Identification of Potential Insect Growth Inhibitor against Aedes aegypti: A Bioinformatics Approach. International Journal of Molecular Sciences, 23(15), 8218. https://doi.org/10.3390/ijms23158218

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop