Molecular Dynamics Simulations Study on the Resistant Mechanism of Insects to Imidacloprid due to Y151-S and R81T Mutations in nAChRs

In-silico selection of peptides for the recognition of imidacloprid

The sensitive detection of pesticides using low-cost receptors designed from peptides can widen their uses in the environmental surveillance for emerging pollutants. In-silico selection of peptides can help accelerate the design of receptor sequence banks for a given target of interest. In this work, we started from Lymnaea stagnalis acetylcholine-binding protein Q55R mutant receptor-imidacloprid complex, available in the PDB databank, to select three primary short peptides (YSP09, DMR12, WQW13 respectively having 9, 12 and 13 amino acids (AA) in length) from the pesticide interacting zones with the A, B and C chains of the nicotinic receptor. Using molecular docking and molecular dynamics (MD) simulations, we showed that the three peptides can form complexes with the target imidacloprid, having energies close to that obtained from a reference RNR12 peptide. Combination of these peptides allowed preparing a new set of longer peptides (YSM21, PSM22, PSW31 and WQA34) that have higher stability and affinity as shown by the MM-PBSA calculations. In particular, the WQA34 peptide displayed an average binding free energy of -6.44±0.27 kcal/mol, which is three times higher than that of the reference RNR12 peptide (-2.29±0.25 kcal/mol) and formed a stable complex with imidacloprid. Furthermore, the dissociation constants (Kd), calculated from the binding free energy, showed that WQA32 (40 μM) has three orders of magnitude lower Kd than the reference RNR12 peptide (3.4 × 104 μM). Docking and RMSD scores showed that the WQA34 peptide is potentially selective to the target imidacloprid with respect to acetamiprid and clothianidin. Therefore, this peptide can be used in wet-lab experiments to prepare a biosensor to selectively detect imidacloprid.

V101I and R81T mutations in the nicotinic acetylcholine receptor β1 subunit are associated with neonicotinoid resistance in Myzus persicae

BACKGROUND

The peach-potato aphid, Myzus persicae (Sulzer) (Hemiptera: Aphididae), is a major pest worldwide. The intensive use of insecticides has led to the development of resistance against neonicotinoid insecticides. The R81T mutation in the nicotinic acetylcholine receptor (nAChR) beta1 subunit is considered a crucial mechanism adaptation to neonicotinoid resistance in M. persicae and Aphis gossypii.

RESULTS

Resistance-related mutations (R81T and V101I) were detected in the imidacloprid-resistant M. persicae AH19 population. The V101I mutation is reported for the first time. The V101I and R81T mutations existed separately, indicating that the two mutations evolved independently. Imidacloprid resistance in the AH19 population was stable without insecticide exposure. Four mutant strains were selected from the population with stable resistance. The resistance of the AH19-T, AH19-I, and AH19-TI strains to imidacloprid, thiamethoxam, and dinotefuran was significantly increased compared with the AH19-W strain. Synergism bioassays showed that the inhibition of three detoxification enzymes did not affect imidacloprid resistance in the AH19-T and AH19-I strains. Expression of nAChR β1 subunits in the AH19-W, AH19-T, and AH19-I strains remained unchanged.

CONCLUSION

The V101I mutation is associated with neonicotinoid resistance in M. persicae. The resistance of the AH19-T and AH19-I strains to neonicotinoids appears to be mainly due to the R81T and V101I mutations, whereas these mutations, together with changes in the cytochrome P450 monooxygenases and nAChR expression may be responsible for the development of neonicotinoid resistance in the AH19-TI strain. © 2021 Society of Chemical Industry.

The relaxation behavior of supercooled and glassy imidacloprid

Employing dielectric spectroscopy, oscillatory shear rheology, and calorimetry, the present work explores the molecular dynamics of the widely used insecticide imidacloprid above and below its glass transition temperature. In its supercooled liquid regime, the applied techniques yield good agreement regarding the characteristic structural (alpha) relaxation times of this material. In addition, the generalized Gemant-DiMarzio-Bishop model provides a good conversion between the frequency-dependent dielectric and shear mechanical responses in its viscous state, allowing for an assessment of imidacloprid's molecular hydrodynamic radius. In order to characterize the molecular dynamics in its glassy regime, we employ several approaches. These include the application of frequency-temperature superposition (FTS) to its isostructural dielectric and rheological responses as well as use of dielectric and calorimetric physical aging and the Adam-Gibbs-Vogel model. While the latter approach and dielectric FTS provide relaxation times that are close to each other, the other methods predict notably longer times that are closer to those reflecting a complete recovery of ergodicity. This seemingly conflicting dissimilarity demonstrates that the molecular dynamics of glassy imidacloprid strongly depends on its thermal history, with high relevance for the use of this insecticide as an active ingredient in technological applications.