Quantitative structure-activity relationship (QSAR) for insecticides: development of predictive in vivo insecticide activity models

Predictionof Environmental FateandToxicityofInsecticidesUsing Multi-Target QSAR Approach

Ecotoxicological risk assessments form the foundation of regulatory decisions for industrial chemicals used in various sectors. In this study, a multi-target-QSAR model established by a backpropagation neural network trained with the Levenberg-Marquardt (LM) algorithm was used to construct a statistically robust and easily interpretable Mt-QSAR model with high external predictability for the simultaneous prediction of the environmental fate in form of octanol-water partition coefficient (LogP), (BCF) and acute oral toxicity in mammals and birds (LD50rat ) and (LD50bird ) for a wide range of chemical structural classes of insecticides. Principal component analysis was performed on descriptors selected by the SW-MLR method, and the selected PCs were used for constructing the SW-MLR-PCA-ANN model. The developed well-trained model (RMSE=0.83, MPE=0.004, CCC=0.82, IIC=0.78, R2 =0.69) was statistically robust as indicated by the external validation parameters (RMSE=0.93, MPE=0.008, CCC=0.77, IIC=0.68, R2 =0.61). The AD of the developed Mt-QSAR model was also defined to identify the most reliable predictions. Finally, the missing values in the dataset for the aforementioned targets were predicted using the constructed Mt-QSAR model. The proposed approach can be used for simultaneous prediction of the environmental fate of new insecticides, especially ones that haven't been tested yet.

In Silico Resources of Drug-Likeness as a Mirror: What Are We Lacking in Pesticide-Likeness?

Unfavorable bioavailability is an important aspect underlying the failure of drug candidates. Computational approaches for evaluating drug-likeness can minimize these risks. Over the past decades, computational approaches for evaluating drug-likeness have sped up the process of drug development and were also quickly derived to pesticide-likeness. As a result of many critical differences between drugs and pesticides, many kinds of methods for drug-likeness cannot be used for pesticide-likeness. Therefore, it is crucial to comprehensively compare and analyze the differences between drug-likeness and pesticide-likeness, which may provide a basis for solving the problems encountered during the evaluation of pesticide-likeness. Here, we systematically collected the recent advances of drug-likeness and pesticide-likeness and compared their characteristics. We also evaluated the current lack of studies on pesticide-likeness, the molecular descriptors and parameters adopted, the pesticide-likeness model on pesticide target organisms, and comprehensive analysis tools. This work may guide researchers to use appropriate methods for developing pesticide-likeness models. It may also aid non-specialists to understand some important concepts in drug-likeness and pesticide-likeness.

Design New Compound of Meisoindigo Derivative as Anti Breast Cancer Based on QSAR Approach

Design new compounds of Meisoindigo using the QSAR approach with semiempirical methods have been successfully made with better biological activity as anti-breast cancer results. The research method used to design the new compound of meisoindigo derivatives uses the semiempirical methods. The research procedure divide into tree step, molecular modeling, QSAR equation analysis, and design of new compounds. The PM3 method was chosen as a better method because it has model results that are more representative of physicochemical aspects. The selection of the best model is made by multilinear regression statistical analysis. New compounds derived from the design are expected to bind to the Cyclin-Dependent Kinase 4 (CDK4) enzyme that helps in the cell cycle to prevent cell division. Based on the design, the compound of (E)-2-(2-oxo-1-(2-(trifluoromethyl) benzyl)indolin-3-ylidene)-N-(quinoline-7-yl) acetamide choose as a new compound with more better biological activity (log 1/IC50 = 6.992) than before (log 1/IC50 = 5.823)

Computational chemistry study of toxicity of some m-tolyl acetate derivatives insecticides and molecular design of structurally related products

Molecular descriptors (including quantum chemical, topological and physicochemical descriptors) were calculated for five m-tolyl derivatives insecticides [namely, carbosulfan (CBS), carbofuran (CBF), isoprocab (IFP), methiocarb (MTC) and isocarbophos (ICP)]. Calculated quantum chemical parameters included the total energy, the electronic energy, the binding energy, the core-core repulsion energy, the heat of formation, the dipole moment and the frontier molecular orbital energies. All the calculated quantum chemical parameters (except dipole moment) exhibited strong correlation with the experimental LD₅₀ values of the studied insecticides (at various Hamiltonians). Calculated topological parameters included the molecular topological index (MTI), polar surface area (PSA), total connectivity (TC), total valence connectivity (TVC), Wiener index (WI), topological diameter (TD) and Balaban index (BI). However, only MTI, PSA, WI and BI exhibited excellent correlation with the toxicological activity of the insecticides. Also among all the calculated physicochemical parameters [logP, surface area (SA), surface volume (SV), hydration energy (E_Hydr), polarizability (PLZ) and refractivity (RFT)], only SV, E_Hydr, PLZ and RFT were useful in establishing quantitative structure activity relationship (QSAR). Application of QSAR indicated that the calculated theoretical LD₅₀ values for the studied insecticides displayed excellent correlation with experimentally derived LD₅₀ values. However, best results were obtained from quantum chemical descriptors under modified neglect of atomic overlap (MNDO). The toxicity profile of the insecticides also correlated strongly with ionization energy, electron affinity, global softness and global harness. Reactive sites of each of the insecticides were established using Fukui function, Huckel charges and HOMO/LUMO diagrams. Six new molecules were proposed and their theoretical activities were estimated. The proposed molecules included 2-methyl-2-(methylthio)-2,3-dihydrobenzofuran-7-yl methylcarbamate, O-methyl O-2-((methylaminooxy)carbonyl)phenyl phosphoramidothioate, 2-((methylaminooxy)carbonyl)phenyl methylcarbamate, 2-(1-(methylthio)ethyl)phenyl methylcarbamate, N-methyl-O-(2-(methylthiooxy) benzoyl) hydroxyl amine and 4-methyl naphthalen-2-yl methylcarbamate. Some of the proposed molecules exhibited negative values of LD₅₀ (indicating extreme toxicity) while two of them exhibited values that are comparable to existing insecticides.

Characterization of acetylcholinesterase from elm left beetle, Xanthogaleruca luteola and QSAR of temephos derivatives against its activity

Insect acetylcholinesterase (AChE) is the principal target for organophosphate (OP) and carbamate (CB) insecticides. In this research, an AChE from third instar larvae of elm left beetle, Xanthogaleruca luteola was purified by affinity chromatography. The enzyme was purified 75.29-fold with a total yield of 8.51%. As shown on denaturing SDS-PAGE, the molecular mass of purified AChE was 70kDa. The enzyme demonstrated maximum activity at pH7 and 35°C. Furthermore, a series of temephos (Tem) derivatives with the general structure of P(O)XP(O) (1-44) were prepared, synthesized and characterized by ³¹P, ¹³C, ¹H NMR and FT-IR spectral techniques. The toxicity of 36 new Tem derivatives was screened on the third instar larvae and the compound compound 1,2 cyclohexane-N,N'-bis(N,N'-piperidine phosphoramidate) exhibited the highest insecticidal potential. The method of kinetic analysis is applied in order to obtain the maximum velocity (V_max), the Michaelis constant (K_m) and the parameters characterizing the inhibition type for inhibitors with >75% mortality in preliminary bioassay. The inhibition mechanism was mixed and inhibitory constant (K_i) was calculated as 4.70μM^-1min^-1 for this compound. Quantitative structure-activity relationship (QSAR) equations of these compounds indicated that the electron orbital energy has major effect on insecticidal properties.

Acetylcholinesterase of the sand fly, Phlebotomus papatasi (Scopoli): construction, expression and biochemical properties of the G119S orthologous mutant

Background

Phlebotomus papatasi vectors zoonotic cutaneous leishmaniasis. Previous expression of recombinant P. papatasi acetylcholinesterase (PpAChE1) revealed 85% amino acid sequence identity to mosquito AChE and identified synthetic carbamates that effectively inhibited PpAChE1 with improved specificity for arthropod AChEs compared to mammalian AChEs. We hypothesized that the G119S mutation causing high level resistance to organophosphate insecticides in mosquitoes may occur in PpAChE1 and may reduce sensitivity to inhibition. We report construction, expression, and biochemical properties of rPpAChE1 containing the G119S orthologous mutation.

Methods

Targeted mutagenesis introduced the G119S orthologous substitution in PpAChE1 cDNA. Recombinant PpAChE1 enzymes containing or lacking the G119S mutation were expressed in the baculoviral system. Biochemical assays were conducted to determine altered catalytic properties and inhibitor sensitivity resulting from the G119S substitution. A molecular homology model was constructed to examine the modeled structural interference with docking of inhibitors of different classes. Genetic tests were conducted to determine if the G119S orthologous codon existed in polymorphic form in a laboratory colony of P. papatasi.

Results

Recombinant PpAChE1 containing the G119S substitution exhibited altered biochemical properties, and reduced inhibition by compounds that bind to the acylation site on the enzyme (with the exception of eserine). Less resistance was directed against bivalent or peripheral site inhibitors, in good agreement with modeled inhibitor docking. Eserine appeared to be a special case capable of inhibition in the absence of covalent binding at the acylation site. Genetic tests did not detect the G119S mutation in a laboratory colony of P. papatasi but did reveal that the G119S codon existed in polymorphic form (GGA + GGC).

Conclusions

The finding of G119S codon polymorphism in a laboratory colony of P. papatasi suggests that a single nucleotide transversion (GGC → AGC) may readily occur, causing rapid development of resistance to organophosphate and phenyl-substituted carbamate insecticides under strong selection. Careful management of pesticide use in IPM programs is important to prevent or mitigate development and fixation of the G119S mutation in susceptible pest populations. Availability of recombinant AChEs enables identification of novel inhibitory ligands with improved efficacy and specificity for AChEs of arthropod pests.

Synthesis and quantitative structure–activity relationship (QSAR) studies of novel rosin-based diamide insecticides

The quantitative structure–activity relationship (QSAR) of two series of rosin-based diamides with insecticidal activity against P. xylostella was studied.

Acetylcholinesterase of the sand fly, Phlebotomus papatasi (Scopoli): cDNA sequence, baculovirus expression, and biochemical properties

Background

Millions of people and domestic animals around the world are affected by leishmaniasis, a disease caused by various species of flagellated protozoans in the genus Leishmania that are transmitted by several sand fly species. Insecticides are widely used for sand fly population control to try to reduce or interrupt Leishmania transmission. Zoonotic cutaneous leishmaniasis caused by L. major is vectored mainly by Phlebotomus papatasi (Scopoli) in Asia and Africa. Organophosphates comprise a class of insecticides used for sand fly control, which act through the inhibition of acetylcholinesterase (AChE) in the central nervous system. Point mutations producing an altered, insensitive AChE are a major mechanism of organophosphate resistance in insects and preliminary evidence for organophosphate-insensitive AChE has been reported in sand flies. This report describes the identification of complementary DNA for an AChE in P. papatasi and the biochemical characterization of recombinant P. papatasi AChE.

Methods

A P. papatasi Israeli strain laboratory colony was utilized to prepare total RNA utilized as template for RT-PCR amplification and sequencing of cDNA encoding acetylcholinesterase 1 using gene specific primers and 3’-5’-RACE. The cDNA was cloned into pBlueBac4.5/V5-His TOPO, and expressed by baculovirus in Sf21 insect cells in serum-free medium. Recombinant P. papatasi acetylcholinesterase was biochemically characterized using a modified Ellman’s assay in microplates.

Results

A 2309 nucleotide sequence of PpAChE1 cDNA [GenBank: JQ922267] of P. papatasi from a laboratory colony susceptible to insecticides is reported with 73-83% nucleotide identity to acetylcholinesterase mRNA sequences of Culex tritaeniorhynchus and Lutzomyia longipalpis, respectively. The P. papatasi cDNA ORF encoded a 710-amino acid protein [GenBank: AFP20868] exhibiting 85% amino acid identity with acetylcholinesterases of Cx. pipiens, Aedes aegypti, and 92% amino acid identity for L. longipalpis. Recombinant P. papatasi AChE1 was expressed in the baculovirus system and characterized as an insect acetylcholinesterase with substrate preference for acetylthiocholine and inhibition at high substrate concentration. Enzyme activity was strongly inhibited by eserine, BW284c51, malaoxon, and paraoxon, and was insensitive to the butyrylcholinesterase inhibitors ethopropazine and iso-OMPA.

Conclusions

Results presented here enable the screening and identification of PpAChE mutations resulting in the genotype for insensitive PpAChE. Use of the recombinant P. papatasi AChE1 will facilitate rapid in vitro screening to identify novel PpAChE inhibitors, and comparative studies on biochemical kinetics of inhibition.

Rational Design of Novel Anti-microtubule Agent (9-Azido-Noscapine) from Quantitative Structure Activity Relationship (QSAR) Evaluation of Noscapinoids

An anticough medicine, noscapine [(S)-3-((R)4-methoxy-6-methyl-5,6,7,8-tetrahydro-[1,3]dioxolo[4,5-g]isoquinolin-5-yl)-6,7-dimethoxyiso-benzofuran-1(3H)-one], was discovered in the authors’ laboratory as a novel type of tubulin-binding agent that mitigates polymerization dynamics of microtubule polymers without changing overall subunit-polymer equilibrium. To obtain systematic insight into the relationship between the structural framework of noscapine scaffold and its antitumor activity, the authors synthesized strategic derivatives (including two new ones in this article). The IC50 values of these analogs vary from 1.2 to 56.0 µM in human acute lymphoblastic leukemia cells (CEM). Geometrical optimization was performed using semiempirical quantum chemical calculations at the 3-21G* level. Structures were in agreement with nuclear magnetic resonance analysis of molecular flexibility in solution and crystal structures. A genetic function approximation algorithm of variable selection was used to generate the quantitative structure activity relationship (QSAR) model. The robustness of the QSAR model ( R2 = 0.942) was analyzed by values of the internal cross-validated regression coefficient ( R2LOO = 0.815) for the training set and determination coefficient ( R2test = 0.817) for the test set. Validation was achieved by rational design of further novel and potent antitumor noscapinoid, 9-azido-noscapine, and reduced 9-azido-noscapine. The experimentally determined value of pIC50 for both the compounds (5.585 M) turned out to be very close to predicted pIC50 (5.731 and 5.710 M).

Perspectives on Non-Animal Alternatives for Assessing Sensitization Potential in Allergic Contact Dermatitis

Skin sensitization remains a major environmental and occupational health hazard. Animal models have been used as the gold standard method of choice for estimating chemical sensitization potential. However, a growing international drive and consensus for minimizing animal usage have prompted the development of in vitro methods to assess chemical sensitivity. In this paper, we examine existing approaches including in silico models, cell and tissue based assays for distinguishing between sensitizers and irritants. The in silico approaches that have been discussed include Quantitative Structure Activity Relationships (QSAR) and QSAR based expert models that correlate chemical molecular structure with biological activity and mechanism based read-across models that incorporate compound electrophilicity. The cell and tissue based assays rely on an assortment of mono and co-culture cell systems in conjunction with 3D skin models. Given the complexity of allergen induced immune responses, and the limited ability of existing systems to capture the entire gamut of cellular and molecular events associated with these responses, we also introduce a microfabricated platform that can capture all the key steps involved in allergic contact sensitivity. Finally, we describe the development of an integrated testing strategy comprised of two or three tier systems for evaluating sensitization potential of chemicals.

QSAR modelling of integrin antagonists using enhanced Bayesian regularised genetic neural networks

Bayesian regularised genetic neural network (BRGNN) has been used for modelling the inhibition activity of 141 biphenylalanine derivatives as integrin antagonists. Three local pattern search (PS) methods, simulated annealing and threshold acceptance were combined with BRGNN in the form of a hybrid genetic algorithm (HGA). The results obtained revealed that PS is a suitable method for improving the ability of BRGNN to break out from the local minima. The proposed HGA technique is able to retrieve important variables from complex systems and nonlinear search spaces for optimisation. Two models with 8-3-1 artificial neural network (ANN) architectures were developed for describing α₄β₇ and α₄β₁ modulatory activities of integrin antagonists. Monte Carlo cross-validation was performed to validate the models and Q₂ values of 0.75 and 0.74 were obtained for α₄β₇ and α₄β₁ inhibitory activities, respectively. The scrambling technique was used for sensitivity analysis of descriptors appearing in ANN models. Frequencies of repetition and sensitivity analysis of molecular descriptors revealed that 3D-Morse descriptors are influential factors for describing α₄β₇ inhibitory activity, while in the case of α₄β₁ inhibitory activity, the Randic shape index, the lowest eigenvalue of the Burden matrix and the number of rotatable bonds are important parameters.