Structure and dynamics of the full-length M1 muscarinic acetylcholine receptor studied by molecular dynamics simulations

The Repellent DEET Potentiates Carbamate Effects via Insect Muscarinic Receptor Interactions: An Alternative Strategy to Control Insect Vector-Borne Diseases

Insect vector-borne diseases remain one of the principal causes of human mortality. In addition to conventional measures of insect control, repellents continue to be the mainstay for personal protection. Because of the increasing pyrethroid-resistant mosquito populations, alternative strategies to reconstitute pyrethroid repellency and knock-down effects have been proposed by mixing the repellent DEET (N,N-Diethyl-3-methylbenzamide) with non-pyrethroid insecticide to better control resistant insect vector-borne diseases. By using electrophysiological, biochemichal, in vivo toxicological techniques together with calcium imaging, binding studies and in silico docking, we have shown that DEET, at low concentrations, interacts with high affinity with insect M1/M3 mAChR allosteric site potentiating agonist effects on mAChRs coupled to phospholipase C second messenger pathway. This increases the anticholinesterase activity of the carbamate propoxur through calcium-dependent regulation of acetylcholinesterase. At high concentrations, DEET interacts with low affinity on distinct M1/M3 mAChR site, counteracting the potentiation. Similar dose-dependent dual effects of DEET have also been observed at synaptic mAChR level. Additionally, binding and in silico docking studies performed on human M1 and M3 mAChR subtypes indicate that DEET only displays a low affinity antagonist profile on these M1/M3 mAChRs. These results reveal a selective high affinity positive allosteric site for DEET in insect mAChRs. Finally, bioassays conducted on Aedes aegypti confirm the synergistic interaction between DEET and propoxur observed in vitro, resulting in a higher mortality of mosquitoes. Our findings reveal an unusual allosterically potentiating action of the repellent DEET, which involves a selective site in insect. These results open exciting research areas in public health particularly in the control of the pyrethroid-resistant insect-vector borne diseases. Mixing low doses of DEET and a non-pyrethroid insecticide will lead to improvement in the efficiency treatments thus reducing both the concentration of active ingredients and side effects for non-target organisms. The discovery of this insect specific site may pave the way for the development of new strategies essential in the management of chemical use against resistant mosquitoes.

Toward activated homology models of the human M1 muscarinic acetylcholine receptor

Structure-based virtual screening offers a good opportunity for the discovery of selective M1 muscarinic acetylcholine receptor (mAChR) agonists for the treatment of Alzheimer's disease. However, no 3-D structure of an M1 mAChR is yet available and the homology models that have been previously reported are only able to identify antagonists in virtual screening experiments. In this study, we generated a homology model of the human M1 mAChR, based on the crystal structure of an M3 mAChR as the template. This initial model was modified, using the agonist-bound crystal structure of a β2-adrenergic receptor as a guide, to give two possible activated structures. The T192 side chain was adjusted in both structures and one of the structures also had the whole of transmembrane (TM) 5 rotated and tilted toward the inner channel of the transmembrane region. The binding sites of all three structures were then refined by induced-fit docking (IFD) with acetylcholine. Virtual screening experiments showed that all three refined models could efficiently differentiate agonists from decoy molecules, with the TM5-modified models also giving good agonist/antagonist selectivity. The whole range of agonists and antagonists was observed to bind within the orthosteric site of the structure obtained by IFD refinement alone, implying that it has inactive state character. In contrast, the two TM5-modified structures were unable to accommodate the antagonists, supporting the proposition that they possess activated state character.

Homology modeling of human muscarinic acetylcholine receptors

We have developed homology models of the acetylcholine muscarinic receptors M₁R-M₅R, based on the β₂-adrenergic receptor crystal as the template. This is the first report of homology modeling of all five subtypes of acetylcholine muscarinic receptors with binding sites optimized for ligand binding. The models were evaluated for their ability to discriminate between muscarinic antagonists and decoy compounds using virtual screening using enrichment factors, area under the ROC curve (AUC), and an early enrichment measure, LogAUC. The models produce rational binding modes of docked ligands as well as good enrichment capacity when tested against property-matched decoy libraries, which demonstrates their unbiased predictive ability. To test the relative effects of homology model template selection and the binding site optimization procedure, we generated and evaluated a naïve M₂R model, using the M₃R crystal structure as a template. Our results confirm previous findings that binding site optimization using ligand(s) active at a particular receptor, i.e. including functional knowledge into the model building process, has a more pronounced effect on model quality than target-template sequence similarity. The optimized M₁R-M₅R homology models are made available as part of the Supporting Information to allow researchers to use these structures, compare them to their own results, and thus advance the development of better modeling approaches.

Activation of the ghrelin receptor is described by a privileged collective motion: a model for constitutive and agonist-induced activation of a sub-class A G-protein coupled receptor (GPCR)

Three homology models of the human ghrelin receptor (GHS-R1a) have been generated from the available X-ray structures of rhodopsin (RHO model), opsin (OPS model) and beta-2 adrenergic receptor (B2 model). The latter was used as a starting point for combined molecular dynamics simulation (MDS) and full atom normal modes analysis (NMA). A low-frequency normal mode (mode 16) perfectly reproduced the intracellular motions observed between B2 and RHO models; in the opposite direction along the same mode, the generated structures are closer to the OPS model, suggesting a direct link with GHS-R1a activation. This was in agreement with motions of the seven transmembranous segments, increase of the solvent accessibility of the 140-ERY-142 sequence, and flip of the Trp276 (C WLP) residue, some features related to GPCRs activation. According to our model, His280 was proposed to stabilize Trp276 in the active state; this was verified by site-directed mutagenesis and biochemical characterization of the resulting H280A and H280S mutants, which were fully functional but sharing an important decrease of their basal activities. Docking performed with short ghrelin derivatives Gly-Ser-Ser ([octa])-Phe-NH (2) and Gly-Ser-Ser ([octa])-Phe-Leu-NH (2) allowed the identification of a robust position of these peptides in the active site of the receptor. This model was refined by MDS and validated by docking experiments performed on a set of 55 ghrelin receptor ligands based on the 1,2,4- triazole scaffold. Finally, NMA performed on the obtained peptide-receptor complex suggested stabilization of the Trp276 residue and of the whole receptor in the active state, preventing the motion observed along mode 16 computed for the unbound receptor. Our results show that NMA offers a powerful approach to study the conformational diversity and the activation mechanism of GPCRs.