Imidacloprid and thiacloprid neonicotinoids bind more favourably to cockroach than to honeybee α6 nicotinic acetylcholine receptor: Insights from computational studies

Thiamethoxam at environmentally relevant concentrations induces neurotoxicity in zebrafish larvae through binding with multiple receptors

Thiamethoxam (THM) is one of the most widely used insecticides globally, which was designed to selectively target nicotinic acetylcholine receptors (nAChRs) in the insect nervous system and is generally considered safe for non-targeted organisms. However, increasing evidence has demonstrated its neurotoxicity in aquatic organisms, though the underlying mechanisms, especially at environmentally relevant concentrations, remain largely unclear. In this study, the swimming distance of zebrafish was significantly shortened by 14.06%-21.64% after exposure to THM at 10-1000 ng/L. This behavioral impairment may result from the damage to nervous and visual systems, as confirmed by notable apoptosis, histological analysis of the eyes, and differential expression of numerous genes. Molecular docking and biomarkers assays found that THM can bind with nAChR and multiple hormone receptors, with binding energies varying from -3.75 to -6.74 kcal/mol. Consequently, the concentrations of a neurotransmitter (acetylcholine) and related hormones (cortisol, triiodothyronine, thyroxine, and thyroid-stimulating hormone) were significantly affected. Further investigations using a weighted gene correlation network and metabolomics suggest that THM may enter the cell via endocytosis and bind with multiple hormone receptors, potentially activating the MAPK signaling pathway. This activation may disrupt purine and pyrimidine metabolism in the cell nucleus, ultimately leading to cell apoptosis and neurotoxicity. This study reveals that THM, even at environmentally relevant concentrations, poses neurological risks to zebrafish and underscore the need for urgent attention to the ecological impacts of THM in aquatic environments.

Discovery of piperonyl-tethered sulfoximines as novel low bee-toxicity aphicides targeting Amelα1/ratβ2 complex

Nicotinic acetylcholine receptor (nAChR) is recognized as a significant insecticide target for neonicotinoids and some agonists. In this study, the nAChR α1 subunit from Apis mellifera was first found to be narrowly tuned to different bee toxicity insecticides, namely, sulfoxaflor (SFX) and flupyradifurone (FPF). Hence, novel sulfoximine derivatives 7a-h were rationally designed and synthesized by introducing a benzo[d][1,3]dioxole moiety into a unique sulfoximine skeleton based on the binding cavity characteristics of Amelα1/ratβ2. The two electrode voltage clamp responses of 7a-h were obviously lower than that of SFX, indicating their potentially low bee toxicity. Besides, representative compounds 7b and 7g exhibited low bee toxicity (LD50 > 11.0 μg/bee at 48 h) revealed by acute contact toxicity bioassays. Molecular modelling results indicated that Ile152, Ala151, and Val160 from honeybee subunit Amelα1 and Lys144 and Trp80 from aphid subunit Mpα1 may be crucial for bee toxicity and aphicidal activity, respectively. These results clarify the toxic mechanism of agonist insecticides on nontargeted pollinators and reveal novel scaffold sulfoximine aphicidal candidates with low bee toxicity. These results will provide a new perspective on the rational design and highly effective development of novel eco-friendly insecticides based on the structure of the nAChR subunit.

Two-Compartmental Toxicokinetic Model Predicts Interspecies Sensitivity Variation of Imidacloprid to Aquatic Invertebrates

Interspecies sensitivity to the same chemical can be several orders of magnitude different. Quantifying toxicologically internal levels and toxicokinetic (TK) parameters is critical in elucidating the interspecies sensitivity. Herein, a two-compartmental TK model was constructed to characterize the uptake, distribution, and elimination kinetics toward interspecies sensitivity to an insecticide, imidacloprid. Imidacloprid exhibited the highest lethality to the insect Chironomus dilutus, followed by Lumbriculus variegatus, Hyalella azteca, and Daphnia magna. Interspecies sensitivity of imidacloprid to these invertebrates varied by ∼1000 folds based on water concentrations (LC50). Remarkably, the sensitivity variation decreased to ∼50 folds based on the internal residues (LR50), highlighting the critical role of TK in interspecies sensitivity. A one-compartmental TK model failed to simulate the bioaccumulation of imidacloprid in these invertebrates except for D. magna. Instead, a two-compartmental model successfully simulated the slow elimination of imidacloprid in the remaining three species by internally distinguishing the highly dynamic (C1) and toxicologically available (C2) fractions. We further showed that the species sensitivity of the invertebrates to imidacloprid was significantly related to C2, demonstrating that C2 was toxicologically available and responsible for the toxicity of imidacloprid. This mechanistic-based model bridged the internal distribution of organic contaminants in small invertebrates and the associated toxic potency.

Identification of sulfonamide compounds active on the insect nervous system: Molecular modeling, synthesis and biological evaluation

Insect nicotinic acetylcholine receptors (nAChRs) are a recognized target for insecticide design. In this work, we have identified, from a structure-based approach using molecular modeling tools, ligands with potential selective activity for pests versus pollinators. A high-throughput virtual screening with the Openeye software was performed using a library from the ZINC database, thiacloprid being used as the target structure. The top sixteen molecules were then docked in α6 cockroach and honeybee homomeric nAChRs to check from a theoretical point of view relevant descriptors in favor of pest selectivity. Among the selected molecules, one original sulfonamide compound has afterward been synthesized, together with various analogs. Two compounds of this family have been shown to behave as activators of the cockroach cholinergic synaptic transmission.

Imidacloprid degrading efficiency of Pseudomonas plecoglossicida MBSB-12 isolated from pesticide contaminated tea garden soil of Assam

Long-term use of toxic pesticides in agricultural grounds has led to adverse effects on the environment and human health. Microbe-mediated biodegradation of pollutants is considered an effective strategy for the removal of contaminants in agricultural and environmental sustainability. Imidacloprid, a neonicotinoid class of pesticides, was widely applied insecticide in the control of pests in agricultural fields including the tea gardens of Assam. Here, native bacteria from imidacloprid contaminating tea garden soils were isolated and screened for imidacloprid degradation efficiency under laboratory conditions. Out of the 30 bacterial isolates, 4 were found to tolerate high concentrations of imidacloprid (25,000 ppm), one of which isolate MBSB-12 showed the highest efficiency for imidacloprid tolerance and utilization as the sole carbon source. Morphological, biochemical, and 16 S ribosomal RNA gene sequencing-based characterization revealed the isolate as Pseudomonas plecoglossicida MBSB-12. The isolate reduced 87% of extractable imidacloprid from the treated soil in 90 days compared to the control soil (without bacterial treatment). High-Resolution Mass Spectrometry (HRMS) analysis indicated imidacloprid breakdown to comparatively less harmful products viz., imidacloprid guanidine olefin [m/z = 209.0510 (M + H)⁺], imidacloprid urea [m/z = 212.0502 (M + H)⁺] and a dechlorinated degraded product of imidacloprid with m/z value 175.0900 (M + H)⁺. Further investigation on the molecular machinery of P. plecoglossicida MBSB-12 involved in the degradation of imidacloprid is expected to provide a better understanding of the degradation pathway.

Mode of Action of Neonicotinoid Insecticides Imidacloprid and Thiacloprid to the Cockroach Pameα7 Nicotinic Acetylcholine Receptor

The functional expression of the cockroach Pameα7 nicotinic acetylcholine receptor subunit has been previously studied, and was found to be able to form a homomeric receptor when expressed in Xenopus laevis oocytes. In this study, we found that the neonicotinoid insecticide imidacloprid is unable to activate the cockroach Pameα7 receptor, although thiacloprid induces low inward currents, suggesting that it is a partial agonist. In addition, the co-application or 5 min pretreatment with 10 µM imidacloprid increased nicotine current amplitudes, while the co-application or 5 min pretreatment with 10 µM thiacloprid decreased nicotine-evoked current amplitudes by 54% and 28%, respectively. This suggesting that these two representatives of neonicotinoid insecticides bind differently to the cockroach Pameα7 receptor. Interestingly, the docking models demonstrate that the orientation and interactions of the two insecticides in the cockroach Pameα7 nAChR binding pocket are very similar. Electrophysiological results have provided evidence to suggest that imidacloprid and thiacloprid could act as modulators of the cockroach Pameα7 receptors.

Fragment-Based Local Coupled Cluster Embedding Approach for the Quantification and Analysis of Noncovalent Interactions: Exploring the Many-Body Expansion of the Local Coupled Cluster Energy

Herein, we introduce a fragment-based local coupled cluster embedding approach for the accurate quantification and analysis of noncovalent interactions in molecular aggregates. Our scheme combines two different expansions of the domain-based local pair natural orbital coupled cluster (DLPNO-CCSD(T)) energy: the many-body expansion (MBE) and the local energy decomposition (LED). The low-order terms in the MBE are initially computed in the presence of an environment that is treated at a low level of theory. Then, LED is used to decompose the energy of each term in the embedded MBE into additive fragment and fragment-pairwise contributions. This information is used to quantify the total energy of the system while providing at the same time in-depth insights into the nature and cooperativity of noncovalent interactions. Two different approaches are introduced and tested, in which the environment is treated at different levels of theory: the local coupled cluster in the Hartree-Fock (LCC-in-HF) method, in which the environment is treated at the HF level; and the electrostatically embedded local coupled cluster method (LCC-in-EE), in which the environment is replaced by point charges. Both schemes are designed to preserve as much as possible the accuracy of the parent local coupled cluster method for total energies, while being embarrassingly parallel and less memory intensive. These schemes appear to be particularly promising for the study of large and complex molecular aggregates at the coupled cluster level, such as condensed phase systems and protein-ligand interactions.

Unraveling individual host-guest interactions in molecular recognition from first principles quantum mechanics: Insights into the nature of nicotinic acetylcholine receptor agonist binding

Drug binding to a protein target is governed by a complex pattern of noncovalent interactions between the ligand and the residues in the protein's binding pocket. Here we introduce a generally applicable, parameter-free, computational method that allows for the identification, quantification, and analysis of the key ligand-residue interactions responsible for molecular recognition. Our strategy relies on Local Energy Decomposition analysis at the "gold-standard" coupled cluster DLPNO-CCSD(T) level. In the study case shown in this paper, nicotine and imidacloprid binding to the nicotinic acetylcholine receptor, our approach provides new insights into how individual amino acids in the active site determine sensitivity and selectivity of the ligands, extending and refining classical pharmacophore hypotheses. By inference, the method is applicable to any kind of host/guest interactions with potential applications in industrial biocatalysis and protein engineering.

Binding of Sulfoxaflor to Aplysia californica-AChBP: Computational Insights from Multiscale Approaches

Structural features and binding properties of sulfoxaflor (SFX) with Ac-AChBP, the surrogate of the insect nAChR ligand binding domain (LBD), are reported herein using various complementary molecular modeling approaches (QM, molecular docking, molecular dynamics, and QM/QM'). The different SFX stereoisomers show distinct behaviors in terms of binding and interactions with Ac-AChBP. Molecular docking and Molecular Dynamics (MD) simulations highlight the specific intermolecular contacts involved in the binding of the different SFX isomers and the relative contribution of the SFX functional groups. QM/QM' calculations provide further insights and a significant refinement of the geometric and energetic contributions of the various residues leading to a preference for the SS and RR stereoisomers. Notable differences in terms of binding interactions are pointed out for the four stereoisomers. The results point out the induced fit of the Ac-AChBP binding site according to the SFX stereoisomer. In this process, the water molecules-mediated contacts play a key role, their energetic contribution being among the most important for the various stereoisomers. In all cases, the interaction with Trp147 is the major binding component, through CH···π and π···π interactions. This study provides a rationale for the binding of SFX to insect nAChR, in particular with respect to the new class of sulfoximine-based insect nAChR competitive modulators, and points out the requirements of various levels of theory for an accurate description of ligand-receptor interactions.

Neonicotinoid insecticides mode of action on insect nicotinic acetylcholine receptors using binding studies

Nicotinic acetylcholine receptors (nAChRs) are the main target of neonicotinoid insecticides, which are widely used in crop protection against insect pests. Electrophysiological and molecular approaches have demonstrated the presence of several nAChR subtypes with different affinities for neonicotinoid insecticides. However, the precise mode of action of neonicotinoids on insect nAChRs remains to be elucidated. Radioligand binding studies with [³H]-α-bungarotoxin and [³H]-imidacloprid have proved instructive in understanding ligand binding interactions between insect nAChRs and neonicotinoid insecticides. The precise binding site interactions have been established using membranes from whole body and specific tissues. In this review, we discuss findings concerning the number of nAChR binding sites against neonicotinoid insecticides from radioligand binding studies on native tissues. We summarize the data available in the literature and compare the binding properties of the most commonly used neonicotinoid insecticides in several insect species. Finally, we demonstrate that neonicotinoid-nAChR binding sites are also linked to biological samples used and insect species.

Residual ground-water levels of the neonicotinoid thiacloprid perturb chemosensing of Caenorhabditis elegans

This study investigated the neurological effects of residual ground-water levels of thiacloprid on the non-target organism Caenorhabditis elegans. Nematodes treated with thiacloprid showed a dose-dependent and significantly increased twitch response at concentrations above 50 ng mL-1 that disabled their forward locomotion in liquid culture. In comparison with untreated controls, 10 ng mL-1 thiacloprid perturbed the chemosensory ability of C. elegans such that the nematodes no longer demonstrated positive chemotaxis towards a NaCl chemo-attractant, reducing their chemotaxis index from +0.48 to near to zero. Nematodes also exhibited a locomotion characteristic of those devoid of chemo-attraction, making significantly more pirouetting turns of ≥90° than the untreated controls. Compared to the untreated controls, expression of the endocytosis-associated gene, Rab-10, was also increased in C. elegans that had developed to adulthood in the presence of 10 ng mL-1 thiacloprid, suggesting their active engagement in increased recycling of affected cellular components, such as their nAChRs. Thus, even residual, low levels of this less potent neonicotinoid that may be found in field ground-water had measurable effects on a beneficial soil organism which may have environmental and ecological implications that are currently poorly understood.

Binary mixtures of neonicotinoids show different transcriptional changes than single neonicotinoids in honeybees (Apis mellifera)

Among the many factors responsible for the decline of bee populations are plant protection products such as neonicotinoids. In general, bees are exposed to not only one but mixtures of such chemicals. At environmental realistic concentrations neonicotinoids may display negative effects on the immune system, foraging activity, learning and memory formation of bees. Neonicotinoids induce alterations of gene transcripts such as nicotinic acetylcholine receptor (nAChR) subunits, vitellogenin, genes of the immune system and genes linked to memory formation. While previous studies focused on individual compounds, the effect of neonicotinoid mixtures in bees is poorly known. Here we investigated the effects of neonicotinoids acetamiprid, clothianidin, imidacloprid and thiamethoxam as single compounds, and binary mixtures thereof in honeybees. We determined transcriptional changes of nAChR subunits and vitellogenin in the brain of experimentally exposed honeybees after exposure up to 72 h. Exposure concentrations were selected on the basis of lowest effect concentrations of the single compounds. Transcriptional induction of nAChRs and vitellogenin was strongest for thiamethoxam, and weakest for acetamiprid. To a large extent, binary mixtures did not show additive transcriptional inductions but they were less than additive. Our data suggest that the joint transcriptional activity of neonicotinoids cannot be explained by concentration addition. The in vivo effects are not only governed by agonistic interaction with nAChRs alone, but are more complex as a result of interactions with other pathways as well. Further studies are needed to investigate the physiological joint effects of mixtures of neonicotinoids and other plant protection products on bees to better understand their joint effects.

Molecular recognition of thiaclopride by Aplysia californica AChBP: new insights from a computational investigation

The binding of thiaclopride (THI), a neonicotinoid insecticide, with Aplysia californica acetylcholine binding protein (Ac-AChBP), the surrogate of the extracellular domain of insects nicotinic acetylcholine receptors, has been studied with a QM/QM' hybrid methodology using the ONIOM approach (M06-2X/6-311G(d):PM6). The contributions of Ac-AChBP key residues for THI binding are accurately quantified from a structural and energetic point of view. The importance of water mediated hydrogen-bond (H-bond) interactions involving two water molecules and Tyr55 and Ser189 residues in the vicinity of the THI nitrile group, is specially highlighted. A larger stabilization energy is obtained with the THI-Ac-AChBP complex compared to imidacloprid (IMI), the forerunner of neonicotinoid insecticides. Pairwise interaction energy calculations rationalize this result with, in particular, a significantly more important contribution of the pivotal aromatic residues Trp147 and Tyr188 with THI through CH···π/CH···O and π-π stacking interactions, respectively. These trends are confirmed through a complementary non-covalent interaction (NCI) analysis of selected THI-Ac-AChBP amino acid pairs.