Jean-Yves Le Questel (CEISAM, Nantes)

Insights into the binding mode of neonicotinoids to the insect nicotinic acetylcholine receptor from Molecular Modeling

Molecular modeling methodologies are efficient and powerful tools to complement experimental observations and predict the properties of original compounds. In this framework, quantum chemistry methods are well-suited for an accurate estimation of key physicochemical properties. Indeed, force field based procedures are not appropriate for a detailed description of ligand-receptor complexes. The combination of various theoretical approaches is therefore necessary in order to get a comprehensive and precise picture of the properties investigated.

In this presentation, through a selection of recent examples studied in our group, the strategy that we use, based on an approach combining a wide range of molecular modeling methods: (QM (DFT), homology modeling, docking, Molecular Dynamics, QM/QM’) and experimental observations, will be presented. This methodology will be applied to the elucidation of conformational and binding properties of neonicotinoids, insecticides targeting nicotinic acetylcholine Receptors (nAChR). First, conformational analyses of the ligands are carried out using Density Functional Theory (DFT) calculations in the isolated state and implicit (SMD) solvent (water, dichloromethane) models. The various conformers are then docked using an induced fit protocol into the binding site of Aplysia californica Acetylcholine Binding Protein (Ac-AChBP), the recognized surrogate of the ligand binding domain of insect nAChR or homology models. Relevant complexes are thereafter submitted to molecular dynamic (MD) simulations. In a last step, QM/QM' calculations are carried out to refine the description (geometries, energies) of the neonicotinoid- nAChR models interactions (Figure 1).

Figure 1.  - Schematic view of the 3D model considered for the QM/QM’ calculations in the case of thiacloprid, one cyanoimine neonicotinoid.

The applied protocol allows pointing out the key components, in terms of amino acid residues, of the binding interactions, according to the chemical structure of the various neonicotinoids. The contribution of weak intermolecular interactions such as halogen bonding, as well as the role of water mediated interactions, is highlighted and quantified. On the whole, the presentation will illustrate the complementarity and the limits of the various methods to delineate and quantify ligand-protein interactions.


Laboratoire CEISAM UMR-CNRS 6230, Université de Nantes
2 rue de la Houssinière, BP 92208