# Mode of Action of Neonicotinoid Insecticides Imidacloprid and Thiacloprid to the Cockroach Pamea7 Nicotinic Acetylcholine Receptor

Alison Cartereau, Emiliane Taillebois, Jean-Yves Le Questel, Steeve H. Thany

► To cite this version:

Alison Cartereau, Emiliane Taillebois, Jean-Yves Le Questel, Steeve H. Thany. Mode of Action of Neonicotinoid Insecticides Imidacloprid and Thiacloprid to the Cockroach Pamea7 Nicotinic Acetylcholine Receptor. International Journal of Molecular Sciences, 2021, 22 (18), pp.9880. 10.3390/ijms22189880. [hal-05137474]

HAL Id: hal-05137474

https://hal.science/hal-05137474v1

Submitted on 1 Jul 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

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

Alison Cartereau<fnref n="1" />, Emiliane Taillebois<fnref n="1" /> ID, Jean-Yves Le Questel<fnref n="2" />,* ID and Steeve H. Thany<fnref n="1" />,∗ ID

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Cartereau, A.; Taillebois, E.; Le Questel, J.-Y.; Thany, S.H. Mode of Action of Neonicotinoid Insecticides Imidacloprid and Thiacloprid to the Cockroach Pameα7 Nicotinic Acetylcholine Receptor. *Int. J. Mol. Sci.* **2021**, *22*, 9880. https://doi.org/10.3390/ijms22189880

Academic Editor: T. G. Emyr Davies

Received: 4 August 2021

Accepted: 8 September 2021

Published: 13 September 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** 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.

**Keywords:** acetylcholine; nicotinic receptors; neonicotinoids; imidacloprid; thiacloprid; molecular modeling; binding studies

## 1. Introduction

Nicotinic acetylcholine receptors (nAChRs) are representative members of the Cys-loop ligand-gated ion channels (LGICs) superfamily, alongside GABA, glycine, and 5-HT3 receptors [1–3]. These neuroreceptors are essential for synaptic transmission processes in both vertebrate and insect nervous systems, receiving the chemical signals between neurons in the nervous system and converting them into electrical outputs. nAChRs are pentameric LGICs, with the five subunits being symmetrically or pseudosymmetrically arranged around a central ion-conducting pore, from homo- or heteropentamers of the related subunits [1–3]. These properties have made nAChRs important targets for human and veterinary drugs, as well as insecticides. However, the functional organization of nAChRs, as well as their diversity in terms of subunit composition and stoichiometry is much better understood in vertebrates than insects [4–6]. Pharmacological studies show that there at least two pharmacologically distinct classes of nAChRs α-bungarotoxin (α-Bgt)-sensitive and -insensitive, and as in vertebrates, insect α7-like subunits are potential candidates to form α-Bgt-sensitive receptors. Insect neuronal nAChRs are of particular interest because they are the main target of neonicotinoid insecticides, which are important in agriculture and veterinary medicine for controlling insect pests and preventing the transmission of insect borne diseases [7–9]. In general, the pharmacological properties of insect nAChRs are studied using electrophysiological approaches with isolated neurons expressing nAChR subtypes. In particular, cockroach neurons from thoracic ganglia and dorsal unpaired

<footnote n="1">1 Université d'Orléans, LBLGC USC INRAE 1328, 45067 Orléans, France; alison.cartereau@univ-orleans.fr (A.C.); emiliane.taillebois@univ-orleans.fr (E.T.)</footnote>

<footnote n="2">2 Université de Nantes, CNRS, CEISAM UMR 6230, 44000 Nantes, France</footnote>

<footnote>* Correspondence: jean-yves.le-questel@univ-nantes.fr (J.-Y.L.Q.); steeve.thany@univ-orleans.fr (S.H.T.); Tel.: +33-(0)2-3841-7058 (J.-Y.L.Q.); +33-(0)2-5112-5563 (S.H.T.)</footnote>
median (DUM) neurons have been used as models to characterize the pharmacological properties of insect native nAChR subtypes [10–13]. Using cockroach thoracic ganglia, two α-Bgt-sensitive nAChR subtypes were characterized as desensitizing (nAChD), which are selectively inhibited by imidacloprid (IMI), and non-desensitizing (nAChN), which are selectively inhibited by methyllycaconitine [10,14]. Moreover, nAChD receptors have been shown to be significantly inhibited by neonicotinoid insecticides, while in contrast, nAChN are activated by neonicotinoids [10]. α-Bgt-sensitive and -insensitive nAChR subtypes were also found in the DUM neurons. Two α-Bgt-insensitive receptors were identified as nAChR1, which are sensitive to IMI and selectively blocked by d-tubocurarine (d-TC), while nAChR2 is inhibited by mecamylamine [13,15–17]. Together, these studies have been used to monitor and analyze the mode of action of neonicotinoid insecticides on insect nAChRs, in order to understand the mechanisms of resistance.

In the aim of characterizing the nAChR subunits involved in the pharmacological properties of cockroach native nAChR subtypes, eight nAChR subunits (seven α and one β subunits) were cloned in the cockroach *Periplaneta americana* [18], while around ten subunits were found in other insect species [19–23]. Indeed, no *Periplaneta americana* α5 (Pameα5) subunit was identified in the cockroach nervous system, suggesting that the cockroach lacked this subunit. Pameα6 and Pameα7 form a cluster with drosophila Dα7 while Pameα8 is close to drosophila Dβ2, suggesting that it corresponds more with a β rather than to an α subunit [18]. The inhibition of cockroach Pameα3, Pameα8 and Pameβ1 expression by antisense oligonucleotides reduced IMI and nicotine current amplitudes, which suggests that they could be included in nAChR1 subtype. Whereas, the inhibition of Pameα1, Pameα2 and Pameβ1 decreased nicotine-activated currents, which correspond to the current-voltage curve of nAChR2 subtype [24]. Moreover, the specific inhibition of the cockroach Pameβ1 subunit resulted in a decrease in nicotine-induced currents, confirming that this subunit accounts for functional nAChR subtypes [18]. These data suggest that more than three neuronal nAChR subtypes were expressed in the cockroach DUM neurons [24]. Based on the number of subunits, and given that the pharmacological properties of these receptors were recorded using whole-cell patch-clamp recording, we hypothesize that nAChD, nAChN, nAChR1 and nAChR2 may be considered as receptor populations with different neuronal nAChR subtypes. While the nature of the subunits has been shown to be crucial in the assembly and function of heteromeric receptors [25], this study aimed to examine the relevant electrophysiological and structural properties of a homomeric cockroach Pameα7 nAChR.

We examined the voltage-clamp recordings and molecular docking studies of the interactions of two relevant representatives of neonicotinoid insecticides (Figure 1), IMI and thiacloprid (THI), with cockroach homomeric Pameα7 nAChR [26]. We demonstrated that the orientation and interactions of the two insecticides in the Pameα7 nAChR binding pocket are very similar. The relative contributions of the various amino acid residues composing the binding site are presented and discussed, according to the nature of the neonicotinoid.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Chemical structures of the neonicotinoids, imidacloprid (left) and thiacloprid (right).</figcaption></figure>
## 2. Results and Discussion

### 2.1. Effect of Neonicotinoid Insecticides on Nicotine-Induced Currents

We previously demonstrated that cockroach nervous system expressed homomeric Pameα7 nAChR [26]. To gain a better understanding of the functional properties of cockroach Pameα7 nAChRs, we have studied the effect of neonicotinoid insecticides, since they have been understood to be agonists of cockroach nAChRs [12,27]. Here, we found that IMI did not induce inward currents on Pameα7 nAChRs. No macroscopic currents were detected within the recorded limits of detection (Figure 2). The lack of agonist action with IMI is not surprising as it is considered a partial agonist of cockroach [12,15] and honeybee nAChRs [28]. Therefore, it may not activate the cockroach Pameα7 receptor, but as suggested, could activate a cockroach receptor composed of Pameα3, Pameα8 and Pameβ1 subunits [24]. Indeed, the blocking expression of these subunits caused a decrease in both nicotine and IMI-induced current amplitudes [24]. Interestingly, THI was able to induce inward currents which were 17% of the nicotine-evoked current amplitudes (Figure 2).

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Currents induced by nicotine (NIC), imidacloprid (IMI) and thiacloprid (THI) on Pameα7 homomeric receptors expressed on *Xenopus laevis* oocytes.</figcaption></figure>

We proposed that IMI and THI showed distinct effects on the cockroach Pameα7 receptor. THI was a partial agonist of the Pameα7 receptor, while IMI showed no agonist effect. Moreover, analysis of the modulatory effect induced by both IMI and THI revealed that co-application or 5 min pretreatment of 10 µM IMI with 10 mM nicotine had a significant effect on nicotine current amplitudes, IMI increased nicotine-evoked current amplitudes to 30% and 40%, respectively (*p* < 0.05, *n* = 12, Figure 3). Moreover, no statistical difference was found between 5 min pretreatment and co-application of IMI (One-way ANOVA, *p* > 0.05, *n* = 10).

The co-application of 10 µM THI with 10 mM nicotine significantly decreased nicotine-induced current amplitudes to 54% of the control currents when they were co-applied (*p* < 0.05, *n* = 6, Figure 4A,B), or 28% after 5 min pretreatment (*p* < 0.05, *n* = 6, Figure 4C,D) with 10 µM THI. These results confirm that the two neonicotinoids could be divided into two subgroups having distinct effects on the Pameα7 receptor.

### 2.2. Sequence Alignments and Homology Modeling

The pairwise sequence alignments of the extracellular domain of the cockroach Pameα7 nAChR subunits against the template crystal structure (PDB ID: 3C79) [29] show 28.4% identity and 36.7% similarity, respectively (Figure 5). The sequence identity value appears slightly below the 30% “threshold” for high accuracy template-based 3D modeling. However, despite this relatively low level of sequence identity, accurate alignment can be obtained in this case as reflected by the sequence similarity value. Based on the sequence alignments, the 3D models of cockroach Pameα7 nAChR homo-pentamers were built using the Prime module [30] of the Schrodinger suite 2014-1 [31] (Figures 6 and 7).
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Effect of co-application and pretreatment of IMI on nicotine-evoked current amplitudes. (A) Currents induced by 10 mM nicotine (NIC) and the co-application of both 10 µM imidacloprid and 10 mM nicotine. (B) Histograms illustrating co-application of 10 µM imidacloprid (IMI) with 10 mM nicotine. (C,D) currents and histograms under 10 mM nicotine and 5 min pretreatment with 10 µM imidacloprid (IMI). Data are mean ± SEM. \* *p* < 0.05. Each histogram represents *n* = 12.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Effects of co-application or pretreatment of THI on nicotine-evoked currents. (A,B) show the effect of co-application of 10 µM THI on currents induced by bath application of 10 mM nicotine. (C,D) illustrate currents-evoked by 10 mM nicotine after pretreatment with 10 µM THI. In all histograms, data are mean ± SEM, *n* = 6 and \* *p* < 0.05. Note that no significant difference was found when in the THI effects with or without 5 min pretreatment. While, (E) represents the inhibition-response curve under co-application of THI with 10 mM nicotine. Each point represents *n* = 6.</figcaption></figure>
<table><caption>**Figure 5.** Sequence alignment of Ac-AChBP and cockroach Pameα7 receptor used to generate the homology model. (*) indicate conserved residues. (:) and (.) indicate conserved and semi-conserved substitutions.</caption><tbody><tr><td>Ac-AChBP</td><td>1</td><td>AHSQANLMRLKSDLFNRSP--MYEPGPTKDDPLITVLGFTLQDIVKADSSTNEVDLVYYEQ</td><td>58</td></tr><tr><td>Pameα7</td><td>1</td><td>SLQGPHEKRILLHTLLDHYNVLERPVANESDPIQLSFGLTLIMQIIDVDEKNOLLITNIWLK</td><td>60</td></tr><tr><td>Ac-AChBP</td><td>59</td><td>QRWKLNSLMWDPNEYGNIIDFRTSAADIWTPDITAYSSTRP-VQVLSPQIAVVTHDGSVM</td><td>117</td></tr><tr><td>Pameα7</td><td>61</td><td>LEWNDVNLRWSTEYGGVKDLRIPPHRIWKPDVLMYNSAEGFDGTYPTNVVVRDSGSCL</td><td>120</td></tr><tr><td>Ac-AChBP</td><td>118</td><td>FIPAQRLSFMCDDPTG-VDSEE GATCAVKFGSWVYSGEIDLKTDTD-QVDLSSYYASSKY</td><td>175</td></tr><tr><td>Pameα7</td><td>121</td><td>YVPPGIFKSTCKIDITWFPFDDQRC EMKFGSWTYDGFQLDLQDDTGGDISSFITNGEW</td><td>180</td></tr><tr><td>Ac-AChBP</td><td>176</td><td>EILSATQTRQVQHYS CCPEPYIDVNLLVKFRERRA</td><td>210</td></tr><tr><td>Pameα7</td><td>181</td><td>DLLGVPGKRNEIYNNCCPEPYIDITFIIIIRRRTL</td><td>215</td></tr></tbody></table>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** 2D interaction diagrams of the ligand-receptor interactions established within the CD interface of the cockroach Pameα7 receptor for (A) IMI and (B) THI.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>Figure 7. *Cont.*</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Figure 7.** Predicted binding modes of IMI and THI with the cockroach Pameα7 receptor viewed according to four orientations: from the top of the homomeric pentameric nAChR (**top** of the **A,B**); from the side of the C and D subunits (**C,D**); with the pyridinic ring (**left**: **E,G**) and the push-pull fragment (**right**: **F,H**), in the plane of the figure, respectively. The docking scores and glide energies computed were −7.2 and −75.1 kcal/mol for IMI and −4.6 and −71.7 kcal/mol for THI, respectively. For clarity, the labelling of the key amino acids involved in intermolecular interactions with the two ligands have been added to (**E–G**), with these amino acids being represented as ball-and-stick models. For the same reasons, only the polar hydrogen atoms are represented. Tyr: tyrosine, Trp: tryptophan, Val: valine, Asn: asparagine, Lys: lysine.</figcaption></figure>

The generated models were energy minimized using the OPLS-2005 force field [32]. The quality of the 3D model, in terms of geometry was assessed using the MolProbity program [33]. Both the final models showed 99.6% residues (except proline and glycine residues), with their φ,ψ torsional angles in the most energetically favored regions of the Ramachandran plot [34]. The protein backbone root mean square deviations (RMSD) of the final cockroach Pameα7 nAChR model was of 0.38Å, using the template structure (3C79) as the reference, and no deviation was observed in the important loops within the binding site.

### 2.3. Investigating the Binding of Imidacloprid and Thiacloprid to the Cockroach Pameα7nAChR through Molecular Docking Studies

In order to determine the binding interactions of neonicotinoids to the homology models of the cockroach Pameα7 nAChR at the atomic level, we carried out docking studies with relevant neonicotinoids. We selected IMI and THI for these investigations, given that; (i) they are forerunners of two chemical classes of neonicotinoid, namely nitro- and cyano-substituted neonicotinoids; (ii) their behavior was notably different from an electrophysiological point of view. The same trends were observed in all five interfaces, and only the results obtained for one of the subunits (CD) will be presented and discussed here.

The geometric and docking parameters are reported in Table S1 of the Supplementary Materials Section. Figure 6A,B show the 2D interaction diagrams of IMI and THI respectively bound to the cockroach Pameα7 nAChR. Whereas, Figure 7 shows complementary 3D views of the predicted binding modes of the two neonicotinoids.

The docking scores and glide energies of the top ranked poses show very similar trends for both IMI and THI insecticides (Table S1, Supplementary Materials and Figure 7).
Indeed, the slight preference for IMI evident for the two parameters (of about 2.5 and 3 kcal mol⁻¹ for docking score and glide energy, respectively) must not be used to draw conclusions and should be considered cautiously. Nevertheless, such studies allow the key interactions established between the IMI and THI insecticides and their target to be delineated. Figure 7A–D show top and side views of the homomeric cockroach Pameα7 nAChR, allowing the location of the ligand binding site to be visualized at the interface between two subunits. Figure 7E–H show the key interactions established between each neonicotinoid and the binding site amino acids. Therefore, the pyridine ring of the two ligands forms π-π interactions with Trp189 of the cockroach Pameα7 nAChR. In addition, and in agreement with the trends observed in the crystallographic structures of neonicotinoid-AChBP complexes [29,35–42], the pyridine ring is consistently involved in water-mediated HB interactions with the main chain groups of amino acid residues, and in the case of the cockroach Pameα7 nAChR, corresponding to Asn147 (CO) and Val159 (NH) (Figure 7). Comparable trends were obtained for the end groups of the push-pull fragment of the two ligands, given that no specific contact is predicted between this moiety of the two insecticides and the receptor residues. In fact, in concordance with experimental data, the stability of the two ligands inside the binding pocket appears to be realized through CH ... π and π ... π interactions, which involve positively polarized CH groups of the five membered saturated ring of the ligands and aromatic amino acid residues of the cockroach Pameα7 nAChR, such as Trp189 and Tyr231 (see Table S1 in Supplementary Materials). It is also worth noting that short contacts are established between the ligand and one sulfur atom of Cys233. These trends are consistent with the key role played by Trp and dicysteine residues in the binding of neonicotinoids insecticides to insect nAChR models.

### 3. Conclusions

Using electrophysiological measurements based on voltage-clamp recordings, we found that, as opposed to THI, IMI could not activate cockroach Pameα7 nAChRs, as no macroscopic currents were detected in our recording conditions. Our results are compatible with the neonicotinoids' modulatory effect on nicotine-induced currents. Indeed, we found that THI decreases nicotine evoked current amplitudes, while IMI induces an increase. Nevertheless, we cannot exclude that IMI has a weak partial agonist action, which was not detected in our experimental conditions. In the same way, the decrease in nicotine-evoked currents under bath application of THI could be associated to a competitive effect on the binding site. Indeed, we found a strong modulatory effect under co-application, which was not present with pre-treatment. Complementary molecular docking investigations led to consistent results with the currently available structural data for the interactions of IMI and THI with 3D nAChR models. Therefore, it provides a starting point for analyzing the relative behavior of neonicotinoid insecticides to cockroach Pameα7 homomeric receptor at the atomic level. The trends that have been highlighted by the electrophysiological results, which suggest that IMI and THI are considered to be positive and negative modulators of the cockroach Pameα7 receptor, respectively. Further electrophysiological, molecular modeling and binding studies on these compounds and other representatives of neonicotinoid insecticides (acetamiprid, clothianidin, thiametoxam) are ongoing in our groups. An exciting prospect for Pameα7 nAChR is to confirm the importance of amino acid residues in the ACh binding domain, particularly in the loop B-C interval in the neonicotinoid selectivity. Finally, the present data reinforce our findings that neonicotinoid insecticides, including nitro- and cyano-substituted neonicotinoids act as agonists and/or modulators of insect homomeric nAChR subtypes. Moreover, insect α-7-like receptors could be used to analyze and compare the binding properties of neonicotinoid insecticides, similar to studies conducted with hybrid receptors.
# 4. Materials and Methods

## 4.1. Compounds

ACh, nicotine, IMI, THI and α-Bgt were purchased from Sigma Chemical Co, (St Quentin, France).

## 4.2. Oocytes Injection

*Xenopus laevis* oocytes were purchased from the University of Rennes, France (CRB Xenope, Rennes, France) and kept in standard oocyte saline (SOS) of the following composition: in mM, 100 NaCl, 2 KCl, 1 MgCl₂, 1.8 CaCl₂ and 5 HEPES, pH 7.5. Stages V and VI oocytes were harvested and defolliculated after treatment with 2 mg·mL⁻¹ collagenase IA (Sigma, Lezennes, France) in Ca²⁺-free SOS solution, supplemented with 0.8 mg·mL⁻¹ trypsin inhibitor. Defolliculated oocytes were injected into nuclei with 2 ng of α7 cDNA cloned in pSGEM [43–45]. Injected oocytes were maintained at 18 °C in SOS solution supplemented with penicillin (100 U·mL⁻¹), streptomycin (100 mg·mL⁻¹), gentamycin (50 mg·mL⁻¹) and sodium pyruvate (2.5 mM).

## 4.3. Voltage-Clamp Recordings

Currents were recorded 4 days after injection, using two microelectrodes filled with 3 M KCl. The oocyte membrane potential was held at −80 mV [46], and perfused continuously with recording buffer at room temperature (20–22 °C). To suppress any possible endogenous muscarinic responses, saline containing 0.5 µM atropine (Sigma, Saint Quentin, France) was employed [47,48]. Oocytes were challenged with a test compound for 5 min intervals to minimize receptor desensitization [44]. To evaluate neonicotinoids modulation of nicotine responses, they were tested with and without 5 min pretreatment. Experimental data were digitized with a Dfidata–1322A A/D converter and later analyzed with pCLAMP (Molecular Devices, Union City, CA, USA). All compound solutions were prepared with the SOS.

## 4.4. Homology Modeling

The amino acid sequence of the cockroach *Periplaneta americana* Pameα7 subunit were extracted from the Uniprot server (www expasy org, access date: 15 July 2021) [49]. The nearest homology of the target sequences were characterized using the BLAST program [50]. The crystal structure of the acetylcholine binding protein (Ac-AChBP) extracted from *Aplysia californica* (Ac) (PDB ID: 3C79) [35] was selected as a template to build the three-dimensional (3D) model of cockroach nAChRs. Indeed, Ac-AChBP is recognized as the surrogate of the ligand binding domain of the extracellular domain of insect nAChRs [35]. The crystal structures of the Ac-AChBP-neonicotinoid complexes (PDB ID: 3C79 and PDB ID: 3C84 for IMI and THI ligands, respectively) were downloaded from the Protein Data Bank (www.pdb.org) [35]. The pairwise sequence alignments were performed to align the target and the template sequence. The 3D homology models were built using the Prime v3.6 (Schrödinger, LLC, New York, NY, USA) module of the Schrodinger suite 2014-1 (Schrödinger, LLC, New York, NY, USA). The rotamers of the conserved amino acid residues are preserved in the homology model in a way that the final 3D model does not significantly deviate from the template structure. The stereochemical quality of the model was verified using MolProbity [33]. Further detail on the methodology based on Molecular Modeling has been described in a previous study [51].

## 4.5. Docking

The chemical structures of the neonicotinoids IMI and THI are shown in Figure 1. The structures were converted to 3D at pH 7.0 ± 0.2 using the LigPrep v3.0 (Schrödinger, LLC, New York, NY, USA) module of the Schrodinger suite 2014-1 (Schrödinger, LLC, New York, NY, USA). The confgen program (Schrödinger, LLC, New York, NY, USA) was then used to retrieve the lowest energy conformer of each insecticide for docking. The docking was carried out using the Glide v6.3 [52] program of the Schrodinger suite 2014-1 (Schrödinger,
LLC, New York, NY, USA). The active site was defined by all residues located in a sphere of 6 Å centered on the ligand. The same residues were selected for receptor grid generation. The extra-precision (XP) [53] mode of the docking algorithm was employed to dock IMI and THI ligands. It is worth remembering that nAChRs are organized as pentamers, that is, there are five identical ligand binding sites and the ligand binds between the cleft of the two subunits. The ligands were therefore docked in all the five interfaces of the cockroach Pameα7 nAChR model. The docking results were validated by comparing the predicted ligand binding modes to the crystallographic Ac-AChBP-neonicotinoid structures.

## 4.6. Statistical Analysis

Data were shown as mean ± SEM and analyzed using Prism 7 (GraphPad Software, La Jolla, CA, USA). Responses to experimental neonicotinoid applications were determined relative to the preceding nicotine control responses in the same experimental conditions. The dose response curves were derived from the fitted curve following the equation: $Y = Imin + (Imax – Imin)/(1 + 10(log(EC50 X)H))$, where Y is the normalized response, Imax and Imin are the maximum and minimum responses, H is the Hill coefficient, EC50 is the concentration giving half the maximum response and X is the logarithm of the compound concentration.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/ijms22189880/s1.

**Author Contributions:** Conceptualization, S.H.T. and J.-Y.L.Q.; methodology, S.H.T., J.-Y.L.Q.; investigation, A.C., J.-Y.L.Q.; writing—original draft preparation, S.H.T., J.-Y.L.Q., E.T., A.C.; writing—review and editing, S.H.T., J.-Y.L.Q.; funding acquisition, S.H.T. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by ANSES, grant number EVALSUB N° 2020/01/214.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** All data from the present study can be found from the corresponding authors.

**Acknowledgments:** Authors acknowledge all students who participated to the study.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

## References

1. Corringer, P.J.; Le Novere, N.; Changeux, J.P. Nicotinic receptors at the amino acid level. *Annu. Rev. Pharmacol. Toxicol.* **2000**, *40*, 431–458. [CrossRef]

2. Changeux, J.P. The nicotinic acetylcholine receptor: A typical ‘allosteric machine’. *Philos. Trans. R. Soc. Lond. B Biol. Sci.* **2018**, *373*, 20170174. [CrossRef] [PubMed]

3. Bertrand, D.; Terry, A.V., Jr. The wonderland of neuronal nicotinic acetylcholine receptors. *Biochem. Pharmacol.* **2018**, *151*, 214–225. [CrossRef] [PubMed]

4. Thany, S.H.; Lenaers, G.; Raymond-Delpech, V.; Sattelle, D.B.; Lapied, B. Exploring the pharmacological properties of insect nicotinic acetylcholine receptors. *Trends Pharmacol. Sci.* **2007**, *28*, 14–22. [CrossRef] [PubMed]

5. Ihara, M.; Buckingham, S.D.; Matsuda, K.; Sattelle, D.B. Modes of Action, Resistance and Toxicity of Insecticides Targeting Nicotinic Acetylcholine Receptors. *Curr. Med. Chem.* **2017**, *24*, 2925–2934. [CrossRef] [PubMed]

6. Taillebois, E.; Cartereau, A.; Jones, A.K.; Thany, S.H. Neonicotinoid insecticides mode of action on insect nicotinic acetylcholine receptors using binding studies. *Pestic. Biochem. Physiol.* **2018**, *151*, 56–66. [CrossRef] [PubMed]

7. Casida, J.E. Pest toxicology: The primary mechanisms of pesticide action. *Chem. Res. Toxicol.* **2009**, *22*, 609–619. [CrossRef] [PubMed]

8. Jeschke, P.; Nauen, R.; Schindler, M.; Elbert, A. Overview of the status and global strategy for neonicotinoids. *J. Agric. Food Chem.* **2011**, *59*, 2897–2908. [CrossRef] [PubMed]

9. Jeschke, P.; Nauen, R.; Beck, M.E. Nicotinic acetylcholine receptor agonists: A milestone for modern crop protection. *Angew. Chem. Int. Ed.* **2013**, *52*, 9464–9485. [CrossRef] [PubMed]
10. Salgado, V.L.; Saar, R. Desensitizing and non-desensitizing subtypes of alpha-bungarotoxin-sensitive nicotinic acetylcholine receptors in cockroach neurons. *J. Insect Physiol.* **2004**, *50*, 867–879. [CrossRef]

11. Tan, J.; Galligan, J.J.; Hollingworth, R.M. Agonist actions of neonicotinoids on nicotinic acetylcholine receptors expressed by cockroach neurons. *Neurotoxicology* **2007**, *28*, 829–842. [CrossRef] [PubMed]

12. Houchat, J.N.; Dissanamossi, B.M.; Landagaray, E.; Mathe-Allainmat, M.; Cartereau, A.; Graton, J.; Lebreton, J.; Le Questel, J.Y.; Thany, S.H. Mode of action of sulfoxaflor on alpha-bungarotoxin-insensitive nAChR1 and nAChR2 subtypes: Inhibitory effect of imidacloprid. *Neurotoxicology* **2019**, *74*, 132–138. [CrossRef] [PubMed]

13. Courjaret, R.; Grolleau, F.; Lapied, B. Two distinct calcium-sensitive and -insensitive PKC up- and down-regulate an alpha-bungarotoxin-resistant nAChR1 in insect neurosecretory cells (DUM neurons). *Eur. J. Neurosci.* **2003**, *17*, 2023–2034. [CrossRef] [PubMed]

14. Salgado, V.L. Antagonist pharmacology of desensitizing and non-desensitizing nicotinic acetylcholine receptors in cockroach neurons. *Neurotoxicology* **2016**, *56*, 188–195. [CrossRef] [PubMed]

15. Courjaret, R.; Lapied, B. Complex intracellular messenger pathways regulate one type of neuronal alpha-bungarotoxin-resistant nicotinic acetylcholine receptors expressed in insect neurosecretory cells (dorsal unpaired median neurons). *Mol. Pharmacol.* **2001**, *60*, 80–91. [CrossRef] [PubMed]

16. Bodereau-Dubois, B.; List, O.; Calas-List, D.; Marques, O.; Communal, P.Y.; Thany, S.H.; Lapied, B. Transmembrane potential polarization, calcium influx, and receptor conformational state modulate the sensitivity of the imidacloprid-insensitive neuronal insect nicotinic acetylcholine receptor to neonicotinoid insecticides. *J. Pharmacol. Exp. Ther.* **2012**, *341*, 326–339. [CrossRef] [PubMed]

17. Thany, S.H.; Courjaret, R.; Lapied, B. Effect of calcium on nicotine-induced current expressed by an atypical alpha-bungarotoxin-insensitive nAChR2. *Neurosci. Lett.* **2008**, *438*, 317–321. [CrossRef]

18. Taillebois, E.; Thany, S.H. Characterization of nicotine acetylcholine receptor subunits in the cockroach *Periplaneta americana* mushroom bodies reveals a strong expression of beta1 subunit: Involvement in nicotine-induced currents. *Arch. Insect Biochem. Physiol.* **2016**, *93*, 40–54. [CrossRef]

19. Sattelle, D.B.; Jones, A.K.; Sattelle, B.M.; Matsuda, K.; Reenan, R.; Biggin, P.C. Edit, cut and paste in the nicotinic acetylcholine receptor gene family of *Drosophila melanogaster*. *Bioessays* **2005**, *27*, 366–376. [CrossRef]

20. Jones, A.K.; Raymond-Delpech, V.; Thany, S.H.; Gauthier, M.; Sattelle, D.B. The nicotinic acetylcholine receptor gene family of the honey bee, *Apis mellifera*. *Genome Res.* **2006**, *16*, 1422–1430. [CrossRef]

21. Jones, A.K.; Grauso, M.; Sattelle, D.B. The nicotinic acetylcholine receptor gene family of the malaria mosquito, *Anopheles gambiae*. *Genomics* **2005**, *85*, 176–187. [CrossRef]

22. Dale, R.P.; Jones, A.K.; Tamborindeguy, C.; Davies, T.G.; Amey, J.S.; Williamson, S.; Wolstenholme, A.; Field, L.M.; Williamson, M.S.; Walsh, T.K.; et al. Identification of ion channel genes in the *Acyrthosiphon pisum* genome. *Insect Mol. Biol.* **2010**, *19* (Suppl. 2), 141–153. [CrossRef]

23. Taillebois, E.; Beloula, A.; Quinchard, S.; Jaubert-Possamai, S.; Daguin, A.; Servent, D.; Tagu, D.; Thany, S.H.; Tricoire-Leignel, H. Neonicotinoid Binding, Toxicity and Expression of Nicotinic Acetylcholine Receptor Subunits in the Aphid *Acyrthosiphon pisum*. *PLoS ONE* **2014**, *9*, e96669. [CrossRef]

24. Sun, H.; Liu, Y.; Li, J.; Cang, X.; Bao, H.; Liu, Z. The potential subunits involved in two subtypes of alpha-Bgt-resistant nAChRs in cockroach dorsal unpaired median (DUM) neurons. *Insect Biochem. Mol. Biol.* **2017**, *81*, 32–40. [CrossRef]

25. Phulera, S.; Zhu, H.; Yu, J.; Claxton, D.P.; Yoder, N.; Yoshioka, C.; Gouaux, E. Cryo-EM structure of the benzodiazepine-sensitive alpha1beta1gamma2S tri-heteromeric GABAA receptor in complex with GABA. *Elife* **2018**, *7*, e39383. [CrossRef]

26. Cartereau, A.; Taillebois, E.; Selvam, B.; Martin, C.; Graton, J.; Le Questel, J.Y.; Thany, S.H. Cloning and Expression of Cockroach alpha7 Nicotinic Acetylcholine Receptor Subunit. *Front. Physiol.* **2020**, *11*, 418. [CrossRef] [PubMed]

27. Thany, S.H. Agonist actions of clothianidin on synaptic and extrasynaptic nicotinic acetylcholine receptors expressed on cockroach sixth abdominal ganglion. *Neurotoxicology* **2009**, *30*, 1045–1052. [CrossRef] [PubMed]

28. Deglise, P.; Grunewald, B.; Gauthier, M. The insecticide imidacloprid is a partial agonist of the nicotinic receptor of honeybee Kenyon cells. *Neurosci. Lett.* **2002**, *321*, 13–16. [CrossRef]

29. Talley, T.T.; Harel, M.; Hibbs, R.E.; Radic, Z.; Tomizawa, M.; Casida, J.E.; Taylor, P. Atomic interactions of neonicotinoid agonists with AChBP: Molecular recognition of the distinctive electronegative pharmacophore. *Proc. Natl. Acad. Sci. USA* **2008**, *105*, 7606–7611. [CrossRef] [PubMed]

30. Schrödinger Release 2014-1: Prime, version 3.5, Schrödinger, LLC, New York, NY, 2014. [CrossRef] [PubMed]

31. Schrödinger Release 2014-1, Schrödinger, LLC, New York, NY, 2014. [CrossRef]

32. Jorgensen, W.L.; Maxwell, D.S.; Tirado-Rives, J. Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids. *J. Am. Chem. Soc.* **1996**, *118*, 11225–11236. [CrossRef]

33. Chen, V.B.; Arendall, W.B., 3rd; Headd, J.J.; Keedy, D.A.; Immormino, R.M.; Kapral, G.J.; Murray, L.W.; Richardson, J.S.; Richardson, D.C. MolProbity: All-atom structure validation for macromolecular crystallography. *Acta Crystallogr. D Biol. Crystallogr.* **2010**, *66*, 12–21. [CrossRef] [PubMed]

34. Ramachandran, G.N.; Ramakrishnan, C.; Sasisekharan, V. Stereochemistry of polypeptide chain configurations. *J. Mol. Biol.* **1963**, *7*, 95–99. [CrossRef]
35. Tomizawa, M.; Casida, J.E. Molecular recognition of neonicotinoid insecticides: The determinants of life or death. *Acc. Chem. Res.* **2009**, *42*, 260–269. [CrossRef] [PubMed]

36. Sixma, T.K.; Smit, A.B. Acetylcholine binding protein (AChBP): A secreted glial protein that provides a high-resolution model for the extracellular domain of pentameric ligand-gated ion channels. *Annu. Rev. Biophys. Biomol. Struct.* **2003**, *32*, 311–334. [CrossRef] [PubMed]

37. Celie, P.H.; van Rossum-Fikkert, S.E.; van Dijk, W.J.; Brejc, K.; Smit, A.B.; Sixma, T.K. Nicotine and carbamylcholine binding to nicotinic acetylcholine receptors as studied in AChBP crystal structures. *Neuron* **2004**, *41*, 907–914. [CrossRef]

38. Celie, P.H.; Klaassen, R.V.; van Rossum-Fikkert, S.E.; van Elk, R.; van Nierop, P.; Smit, A.B.; Sixma, T.K. Crystal structure of acetylcholine-binding protein from *Bulinus truncatus* reveals the conserved structural scaffold and sites of variation in nicotinic acetylcholine receptors. *J. Biol. Chem.* **2005**, *280*, 26457–26466. [CrossRef]

39. Hansen, S.B.; Sulzenbacher, G.; Huxford, T.; Marchot, P.; Taylor, P.; Bourne, Y. Structures of *Aplysia* AChBP complexes with nicotinic agonists and antagonists reveal distinctive binding interfaces and conformations. *EMBO J.* **2005**, *24*, 3635–3646. [CrossRef] [PubMed]

40. Ihara, M.; Okajima, T.; Yamashita, A.; Oda, T.; Hirata, K.; Nishiwaki, H.; Morimoto, T.; Akamatsu, M.; Ashikawa, Y.; Kuroda, S.; et al. Crystal structures of *Lymnaea stagnalis* AChBP in complex with neonicotinoid insecticides imidacloprid and clothianidin. *Invert. Neurosci.* **2008**, *8*, 71–81. [CrossRef] [PubMed]

41. Rucktooa, P.; Smit, A.B.; Sixma, T.K. Insight in nAChR subtype selectivity from AChBP crystal structures. *Biochem. Pharmacol.* **2009**, *78*, 777–787. [CrossRef] [PubMed]

42. Ulens, C.; Akdemir, A.; Jongejan, A.; van Elk, R.; Bertrand, S.; Perrakis, A.; Leurs, R.; Smit, A.B.; Sixma, T.K.; Bertrand, D.; et al. Use of acetylcholine binding protein in the search for novel alpha7 nicotinic receptor ligands. In silico docking, pharmacological screening, and X-ray analysis. *J. Med. Chem.* **2009**, *52*, 2372–2383. [CrossRef] [PubMed]

43. Bertrand, D.; Ballivet, M.; Gomez, M.; Bertrand, S.; Phannavong, B.; Gundelfinger, E.D. Physiological properties of neuronal nicotinic receptors reconstituted from the vertebrate beta 2 subunit and Drosophila alpha subunits. *Eur. J. Neurosci.* **1994**, *6*, 869–875. [CrossRef] [PubMed]

44. Ihara, M.; Matsuda, K.; Otake, M.; Kuwamura, M.; Shimomura, M.; Komai, K.; Akamatsu, M.; Raymond, V.; Sattelle, D.B. Diverse actions of neonicotinoids on chicken alpha7, alpha4beta2 and Drosophila-chicken SADbeta2 and ALSbeta2 hybrid nicotinic acetylcholine receptors expressed in *Xenopus laevis* oocytes. *Neuropharmacology* **2003**, *45*, 133–144. [CrossRef]

45. Couturier, S.; Bertrand, D.; Matter, J.M.; Hernandez, M.C.; Bertrand, S.; Millar, N.; Valera, S.; Barkas, T.; Ballivet, M. A neuronal nicotinic acetylcholine receptor subunit (alpha 7) is developmentally regulated and forms a homo-oligomeric channel blocked by alpha-BTX. *Neuron* **1990**, *5*, 847–856. [CrossRef]

46. Taylor-Wells, J.; Hawkins, J.; Colombo, C.; Bermudez, I.; Jones, A.K. Cloning and functional expression of intracellular loop variants of the honey bee (*Apis mellifera*) RDL GABA receptor. *Neurotoxicology* **2017**, *60*, 207–213. [CrossRef]

47. Matsuda, K.; Buckingham, S.D.; Freeman, J.C.; Squire, M.D.; Baylis, H.A.; Sattelle, D.B. Effects of the alpha subunit on imidacloprid sensitivity of recombinant nicotinic acetylcholine receptors. *Br. J. Pharmacol.* **1998**, *123*, 518–524. [CrossRef]

48. Matsuda, K.; Shimomura, M.; Kondo, Y.; Ihara, M.; Hashigami, K.; Yoshida, N.; Raymond, V.; Mongan, N.P.; Freeman, J.C.; Komai, K.; et al. Role of loop D of the alpha7 nicotinic acetylcholine receptor in its interaction with the insecticide imidacloprid and related neonicotinoids. *Br. J. Pharmacol.* **2000**, *130*, 981–986. [CrossRef] [PubMed]

49. Gasteiger, E.; Gattiker, A.; Hoogland, C.; Ivanyi, I.; Appel, R.D.; Bairoch, A. ExPASy: The proteomics server for in-depth protein knowledge and analysis. *Nucleic Acids Res.* **2003**, *31*, 3784–3788. [CrossRef]

50. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. *J. Mol. Biol.* **1990**, *215*, 403–410. [CrossRef]

51. Selvam, B.; Graton, J.; Laurent, A.D.; Alamidine, Z.; Mathe-Allainmat, M.; Lebreton, J.; Coqueret, O.; Olivier, C.; Thany, S.H.; Le Questel, J.Y. Imidacloprid and thiacloprid neonicotinoids bind more favourably to cockroach than to honeybee alpha6 nicotinic acetylcholine receptor: Insights from computational studies. *J. Mol. Graph. Model.* **2015**, *55*, 1–12. [CrossRef] [PubMed]

52. Friesner, R.A.; Banks, J.L.; Murphy, R.B.; Halgren, T.A.; Klicic, J.J.; Mainz, D.T.; Repasky, M.P.; Knoll, E.H.; Shelley, M.; Perry, J.K.; et al. Glide: A new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. *J. Med. Chem.* **2004**, *47*, 1739–1749. [CrossRef] [PubMed]

53. Friesner, R.A.; Murphy, R.B.; Repasky, M.P.; Frye, L.L.; Greenwood, J.R.; Halgren, T.A.; Sanschagrin, P.C.; Mainz, D.T. Extra precision glide: Docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes. *J. Med. Chem.* **2006**, *49*, 6177–6196. [CrossRef] [PubMed]