Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2013 Jun;27(6):525-38.
doi: 10.1007/s10822-013-9660-8. Epub 2013 Jun 28.

On homology modeling of the M₂ muscarinic acetylcholine receptor subtype

Jan Jakubík  1 , Alena Randáková, Vladimír Doležal
Affiliations

On homology modeling of the M₂ muscarinic acetylcholine receptor subtype

Jan Jakubík et al. J Comput Aided Mol Des. 2013 Jun.

Abstract

Twelve homology models of the human M2 muscarinic receptor using different sets of templates have been designed using the Prime program or the modeller program and compared to crystallographic structure (PDB:3UON). The best models were obtained using single template of the closest published structure, the M3 muscarinic receptor (PDB:4DAJ). Adding more (structurally distant) templates led to worse models. Data document a key role of the template in homology modeling. The models differ substantially. The quality checks built into the programs do not correlate with the RMSDs to the crystallographic structure and cannot be used to select the best model. Re-docking of the antagonists present in crystallographic structure and relative binding energy estimation by calculating MM/GBSA in Prime and the binding energy function in YASARA suggested it could be possible to evaluate the quality of the orthosteric binding site based on the prediction of relative binding energies. Although estimation of relative binding energies distinguishes between relatively good and bad models it does not indicate the best one. On the other hand, visual inspection of the models for known features and knowledge-based analysis of the intramolecular interactions allows an experimenter to select overall best models manually.

PubMed Disclaimer

Figures

Fig. 1
Fig. 1
Alignments of templates to target structure. Alignment of templates for homology modeling labeled by their PDB entry code to the target sequence of the human M2 muscarinic acetylcholine receptor. Stars denote conserved and dots consensual residues. Colors denote secondary structure: red—helix; white—coil; yellow—strand; green—turn. Secondary structure of the target was predicted by PsiPred (http://bioinf.cs.ucl.ac.uk/psipred/) and secondary structures of templates were taken from respective crystal structures. For orientation transmembrane (TM) helices, inner (i) and outer (o) loops are indicated
Fig. 2
Fig. 2
Twelve homology models of the M2 receptor superposed on the 3UON crystal structure. Homology models (color) of the M2 receptor based on the templates listed in Table 1 are superposed using MUSTANG [27] implemented in YASARA on the crystal structure 3UON (gray). Orientation: extracellular site up, TM VI and TM VII front. Colors: purple—α-helix; yellow—β-sheet; cyan—turn; white—coil. RMSDs of the models to the target structure (3UON) are listed in Table 4
Fig. 3
Fig. 3
Correlation of calculated binding energies and RMSD of redocked QNB to 3UON. Calculated YASARA binding energies (top), AutoDock energies (middle) and Prime MMGB/SA energies (bottom) of poses from QNB re-docking to 3UON are plotted against their RMSDs. Detail with poses of RMSD of ligand smaller than 3 Å is on left and all poses are on right. Correlation coefficients are indicated in the graphs
Fig. 4
Fig. 4
Calculated binding energies of muscarinic antagonists docked to homology models. Calculated binding energies of the antagonists QNB (squares), NMQNB (downward triangles), NMS (circles) and atropine (upward triangles) in YASARA (y-axis) are plotted against the binding energies calculated in Prime MMGB/SA (x-axis). Good poses are black, and bad poses are red. Closed symbols denote the best poses
Fig. 4
Fig. 4
Calculated binding energies of muscarinic antagonists docked to homology models. Calculated binding energies of the antagonists QNB (squares), NMQNB (downward triangles), NMS (circles) and atropine (upward triangles) in YASARA (y-axis) are plotted against the binding energies calculated in Prime MMGB/SA (x-axis). Good poses are black, and bad poses are red. Closed symbols denote the best poses
Fig. 5
Fig. 5
The best QNB docking poses of 12 homology models superposed on the 3UON crystal structure. View from the extracellular site, TM II down, TM VI and TM VII up, of the best docking poses, according to the binding energy estimates (Fig. 4), of QNB (green carbons) and residues of the orthosteric binding site (color) superposed on the crystal structure of 3UON (gray). Colors: Cyan—carbon; red—oxygen; blue—nitrogen; white—hydrogen; yellow—hydrogen bonds. The residue labels correspond to the M2 sequence. The calculated RMSDs of QNB and residues of the orthosteric binding site of the models to the target structure (3UON) are shown in Table 5
See this image and copyright information in PMC

Similar articles

See all similar articles

Cited by

See all "Cited by" articles

References

    1. Jakubík J, Doležal V, El-Fakahany EE, Janíčková H, Randáková A, Šantrůčková E. Perspective for design of selective muscarinic agonists. In: Babušíková E, Dobrota D, Lehotský J, editors. New frontiers in molecular mechanisms in psychiatric and neurologic disorders. Martin: Jesseniova lekárska fakulta UK; 2011.
    1. Jakubík J, El-Fakahany EE. Allosteric modulation of muscarinic acetylcholine receptors. Pharmaceuticals. 2010;9:2838–2860. doi: 10.3390/ph3092838. - DOI - PMC - PubMed
    1. Christopoulos A, El-Fakahany EE. Novel persistent activation of muscarinic M1 receptors by xanomeline. Eur J Pharmacol. 1997;334:R3–R4. doi: 10.1016/S0014-2999(97)01162-X. - DOI - PubMed
    1. Caffrey M. Membrane protein crystallization. J Struct Biol. 2003;142:108–132. doi: 10.1016/S1047-8477(03)00043-1. - DOI - PubMed
    1. Haga K, Kruse AC, Asada H, Yurugi-Kobayashi T, Shiroishi M, et al. Structure of the human M2 muscarinic acetylcholine receptor bound to an antagonist. Nature. 2012;482:547–551. doi: 10.1038/nature10753. - DOI - PMC - PubMed

Publication types

Substances

LinkOut - more resources