Homology Modeling and Validation of the Human M1 Muscarinic Acetylcholine Receptor

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.