Structure Of Bovine Rhodopsin Research Articles

Binding mode analysis and enrichment studies on homology models of the human histamine H4 receptor

Ligand-supported homology models of the human histamine H4 receptor (hH4R) were developed based on the crystal structure of bovine rhodopsin and different known H4 ligands (histamine, OUP-16, JNJ7777120). Enrichment tests were performed to analyze whether our hH4R models can select known actives from random decoys. The impact of receptor conformation and the effect of different sets of random decoys, docking methods (FlexX, FlexX-Pharm) and scoring functions (FlexX-Score, D-Score, PMF-Score, G-Score, ChemScore) were investigated. We found that two agonists (histamine and OUP-16) form complementary interactions with Asp94 (3.32), Glu182 (5.46) and Thr323 (6.55), whereas JNJ7777120 interacts with Asp94 (3.32) and Glu182 (5.46) only. These results suggest a role of Thr323 (6.55) in ligand binding and presumably also in receptor activation. The models optimized in the presence of an agonist (histamine) and an antagonist (JNJ7777120) were compared in more detail. We conclude that the ligand used in the model building process can significantly influence the efficacy of virtual screening.

Generation of a homology model of the human histamine H3 receptor for ligand docking and pharmacophore-based screening

The human histamine H(3) receptor (hH(3)R) is a G-protein coupled receptor (GPCR), which modulates the release of various neurotransmitters in the central and peripheral nervous system and therefore is a potential target in the therapy of numerous diseases. Although ligands addressing this receptor are already known, the discovery of alternative lead structures represents an important goal in drug design. The goal of this work was to study the hH(3)R and its antagonists by means of molecular modelling tools. For this purpose, a strategy was pursued in which a homology model of the hH(3)R based on the crystal structure of bovine rhodopsin was generated and refined by molecular dynamics simulations in a dipalmitoylphosphatidylcholine (DPPC)/water membrane mimic before the resulting binding pocket was used for high-throughput docking using the program GOLD. Alternatively, a pharmacophore-based procedure was carried out where the alleged bioactive conformations of three different potent hH(3)R antagonists were used as templates for the generation of pharmacophore models. A pharmacophore-based screening was then carried out using the program Catalyst. Based upon a database of 418 validated hH(3)R antagonists both strategies could be validated in respect of their performance. Seven hits obtained during this screening procedure were commercially purchased, and experimentally tested in a [(3)H]N(alpha)-methylhistamine binding assay. The compounds tested showed affinities at hH(3)R with K ( i ) values ranging from 0.079 to 6.3 muM.

Critical Role for the Second Extracellular Loop in the Binding of Both Orthosteric and Allosteric G Protein-coupled Receptor Ligands

The second extracellular (E2) loop of G protein-coupled receptors (GPCRs) plays an essential but poorly understood role in the binding of non-peptidic small molecules. We have utilized both orthosteric ligands and allosteric modulators of the M2 muscarinic acetylcholine receptor, a prototypical Family A GPCR, to probe possible E2 loop binding dynamics. We developed a homology model based on the crystal structure of bovine rhodopsin and predicted novel cysteine substitutions that should dramatically reduce E2 loop flexibility via disulfide bond formation and significantly inhibit the binding of both types of ligands. This prediction was validated experimentally using radioligand binding, dissociation kinetics, and cell-based functional assays. The results argue for a flexible "gatekeeper" role of the E2 loop in the binding of both allosteric and orthosteric GPCR ligands.

Systematic Analysis of the Entire Second Extracellular Loop of the V1a Vasopressin Receptor: KEY RESIDUES, CONSERVED THROUGHOUT A G-PROTEIN-COUPLED RECEPTOR FAMILY, IDENTIFIED

The roles of extracellular residues of G-protein-coupled receptors (GPCRs) are not well defined compared with residues in transmembrane helices. Nevertheless, it has been established that extracellular domains of both peptide-GPCRs and amine-GPCRs incorporate functionally important residues. Extracellular loop 2 (ECL2) has attracted particular interest, because the x-ray structure of bovine rhodopsin revealed that ECL2 projects into the binding crevice within the transmembrane bundle. Our study provides the first comprehensive investigation into the role of the individual residues comprising the entire ECL2 domain of a small peptide-GPCR. Using the V(1a) vasopressin receptor, systematic substitution of all of the ECL2 residues by Ala generated 30 mutant receptors that were characterized pharmacologically. The majority of these mutant receptor constructs (24 in total) had essentially wild-type ligand binding and intracellular signaling characteristics, indicating that these residues are not critical for normal receptor function. However, four aromatic residues Phe(189), Trp(206), Phe(209), and Tyr(218) are important for agonist binding and receptor activation and are highly conserved throughout the neurohypophysial hormone subfamily of peptide-GPCRs. Located in the middle of ECL2, juxtaposed to the highly conserved disulfide bond, Trp(206) and Phe(209) project into the binding crevice. Indeed, Phe(209) is part of the Cys-X-X-X-Ar (where Ar is an aromatic residue) motif, which is well conserved in both peptide-GPCRs and amine-GPCRs. In contrast, Phe(189) and Tyr(218), located at the extreme ends of ECL2, may be important for determining the position of the ECL2 cap over the binding crevice. This study provides mechanistic insight into the roles of highly conserved ECL2 residues.

Molecular modeling of histamine H3 receptor and QSAR studies on arylbenzofuran derived H3 antagonists

Histamine H3 receptors are presynaptic autoreceptors found in both central and peripheral nervous systems of many species. The central effects of these receptors suggest a potential therapeutic role for their antagonists in treatment of several neurological disorders such as epilepsy, schizophrenia, Alzheimer's and Parkinson's diseases. The purpose of this study was to identify the structural requirements for H3 antagonistic activity via quantitative structure–activity relationship (QSAR) studies and receptor modeling/docking techniques. A combination of partial least squares (PLS) and genetic algorithm (GA) was used in the QSAR approach to select the structural descriptors relevant to the receptor binding affinity of a series of 58 H3 antagonists. The descriptors were selected out of a pool of >1000 descriptors calculated by DRAGON, Hyperchem and ACD labs suite of programs. The resulting QSAR models for rat and human H3 binding affinities were validated using different strategies. QSAR models generated in the current work suggested the role of charge transfer interactions in the ligand–receptor interaction verified using the molecular modeling of the receptor and docking two antagonists to the binding site. The 3D model of human H3 receptor was built based on bovine rhodopsin structure and evaluated by molecular dynamics (MD) simulation in a mixed water–vacuum–water environment. The results were indicative of the stability of the model relating the observed structural changes during the MD simulation to the suggested ligand–receptor interactions. The results of this investigation are expected to be useful in the process of design and development of new potent H3 receptor antagonists.

Understanding interactions of gastric inhibitory polypeptide (GIP) with its G‐protein coupled receptor through NMR and molecular modeling

Gastric inhibitory polypeptide (GIP, or glucose-dependent insulinotropic polypeptide) is a 42-amino acid incretin hormone moderating glucose-induced insulin secretion. Antidiabetic therapy based on GIP holds great promise because of the fact that its insulinotropic action is highly dependent on the level of glucose, overcoming the sideeffects of hypoglycemia associated with the current therapy of Type 2 diabetes. The truncated peptide, GIP(1-30)NH2, has the same activity as the full length native peptide. We have studied the structure of GIP(1-30)NH2 and built a model of its G-protein coupled receptor (GPCR). The structure of GIP(1-30)NH2 in DMSO-d6 and H2O has been studied using 2D NMR (total correlation spectroscopy (TOCSY), nuclear overhauser effect spectroscopy (NOESY), double quantum filtered-COSY (DQF-COSY), 13C-heteronuclear single quantum correlation (HSQC) experiments, and its conformation built by MD simulations with the NMR data as constraints. The peptide in DMSO-d6 exhibits an alpha-helix between residues Ile12 and Lys30 with a discontinuity at residues Gln19 and Gln20. In H2O, the alpha-helix starts at Ile7, breaks off at Gln19, and then continues right through to Lys30. GIP(1-30)NH2 has all the structural features of peptides belonging to family B1 GPCRs, which are characterized by a coil at the N-terminal and a long C-terminal alpha-helix with or without a break. A model of the seven transmembrane (TM) helices of the GIP receptor (GIPR) has been built on the principles of comparative protein modeling, using the crystal structure of bovine rhodopsin as a template. The N-terminal domain of GIPR has been constructed from the NMR structure of the N-terminal of corticoptropin releasing factor receptor (CRFR), a family B1 GCPR. The intra and extra cellular loops and the C-terminal have been modeled from fragments retrieved from the PDB. On the basis of the experimental data available for some members of family B1 GPCRs, four pairs of constraints between GIP(1-30)NH2 and its receptor were used in the FTDOCK program, to build the complete model of the GIP(1-30)NH2:GIPR complex. The model can rationalize the various experimental observations including the potency of the truncated GIP peptide. This work is the first complete model at the atomic level of GIP(1-30)NH2 and of the complex with its GPCR.

Multiple Residues in the Second Extracellular Loop Are Critical for M3 Muscarinic Acetylcholine Receptor Activation

Recent studies suggest that the second extracellular loop (o2 loop) of bovine rhodopsin and other class I G protein-coupled receptors (GPCRs) targeted by biogenic amine ligands folds deeply into the transmembrane receptor core where the binding of cis-retinal and biogenic amine ligands is known to occur. In the past, the potential role of the o2 loop in agonist-dependent activation of biogenic amine GPCRs has not been studied systematically. To address this issue, we used the M(3) muscarinic acetylcholine receptor (M3R), a prototypic class I GPCR, as a model system. Specifically, we subjected the o2 loop of the M3R to random mutagenesis and subsequently applied a novel yeast genetic screen to identity single amino acid substitutions that interfered with M3R function. This screen led to the recovery of about 20 mutant M3Rs containing single amino acid changes in the o2 loop that were inactive in yeast. In contrast, application of the same strategy to the extracellular N-terminal domain of the M3R did not yield any single point mutations that disrupted M3R function. Pharmacological characterization of many of the recovered mutant M3Rs in mammalian cells, complemented by site-directed mutagenesis studies, indicated that the presence of several o2 loop residues is important for efficient agonist-induced M3R activation. Besides the highly conserved Cys(220) residue, Gln(207), Gly(211), Arg(213), Gly(218), Ile(222), Phe(224), Leu(225), and Pro(228) were found to be of particular functional importance. In general, mutational modification of these residues had little effect on agonist binding affinities. Our findings are therefore consistent with a model in which multiple o2 loop residues are involved in stabilizing the active state of the M3R. Given the high degree of structural homology found among all biogenic amine GPCRs, our findings should be of considerable general relevance.

Targeting HIV‐1 Through Molecular Modeling and Docking Studies of CXCR4: Leads for Therapeutic Development

The chemokine receptor CXCR4 is the receptor for several chemokines and major co-receptor for X4 human immunodeficiency virus type-1 strains entry into cell. A three-dimensional model of human CXCR4 was developed by homology modeling using the high-resolution bovine rhodopsin structure as template. Interactions between CXCR4 and flavonoids were investigated using in silico docking studies. The results underscore the potential of these compounds that they may become important new antiviral drugs to combat AIDS. It is worth mentioning also that apart from these existing flavonoids, there are many new compounds that may also be useful as topical agents to inactivate virus, or may act as adjuvants with other antiviral drugs.

Identification of the Intracellular Region of the Leukotriene B4 Receptor Type 1 That Is Specifically Involved in Gi Activation

Many G-protein-coupled receptors can activate more than one G-protein subfamily member. Leukotriene B4 receptor type 1 (BLT1) is a high affinity G-protein-coupled receptors for leukotriene B4 functioning in host defense, inflammation, and immunity. Previous studies have shown that BLT1 utilizes different G-proteins (the Gi family and G16 G-proteins) in mediating diverse cellular events and that truncation of the cytoplasmic tail of BLT1 does not impair activation of Gi and G16 proteins. To determine responsive regions of BLT1 for G-protein coupling, we performed an extensive mutagenesis study of its intracellular loops. Three intracellular loops (i1, i2, and i3) of BLT1 were found to be important for both Gi and G16 coupling, as judged by Gi-dependent guanosine 5'-(gamma-thio) triphosphate (GTPgammaS) binding and G16-dependent inositol phosphate accumulation assays. The i3-1 mutant, with a mutation at the i3 amino terminus, exhibited greatly reduced GTPgammaS binding but intact inositol phosphate accumulation triggered by leukotriene B4 stimulation. These results suggest that the i3-1 region is required only for Gi activation. Moreover, in the i3-1 mutant, the deficiency in Gi activation was accompanied by a loss of the high affinity leukotriene B4 binding state seen with the wild type receptor. A three-dimensional model of BLT1 constructed based on the structure of bovine rhodopsin suggests that the i3-1 region may consist of the cytoplasmic end of the transmembrane helix V, which protrudes the helix into the cytoplasm. From mutational studies and three-dimensional modeling, we propose that the extended cytoplasmic helix connected to the transmembrane helix V of BLT1 might be a key region for selective activation of Gi proteins.

Exploring the binding site of the human muscarinic M3 receptor: Homology modeling and docking study

AbstractThe human muscarinic M3 receptor (hM3) and its interactions with selective agonists and antagonists were investigated by means of combined homology and docking approach. Also, two pharmacophoric models for the hM3 agonist and antagonist binding sites were proposed. The three‐dimensional (3D) structure of hM3 receptor was modeled based on the high‐resolution X‐ray structure of bovine rhodopsin from the Protein Data Bank (PDB). To validate the reliability of the model obtained, the main chain torsion angles phi (Ψ) and psi (Φ) were examined in a Ramachandran plot, and all omega angles were measured for peptidic bond planarity. The characteristics of the active site, the position, and the orientation of ligands in situ, as well as the binding modes of the representative agonists and antagonists, were analyzed by applying a molecular docking technique using the AutoDock 3.0.5 program. Specific interactions responsible for recognition of the hM3 receptor, like ionic bond formed between protonated amine of the ligands and the Asp3.6 side chain were identified. Structure–reactivity relationships have been explained by analyzing the 3D structure of the hM3 model and the ligand conformations resulted from molecular docking simulation. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007

Identification of Leu276 of the S1P1Receptor and Phe263 of the S1P3Receptor in Interaction with Receptor Specific Agonists by Molecular Modeling, Site-Directed Mutagenesis, and Affinity Studies

Sphingosine-1-phosphate (S1P) receptor agonists are novel immunosuppressive agents. The selectivity of S1P1 against S1P3 is strongly correlated with lymphocyte sequestration and minimum acute toxicity and bradycardia. This study describes molecular modeling, site-directed mutagenesis, and affinity studies exploring the molecular basis for selectivity between S1P1 and S1P3 receptors. Computational models of human S1P1 and S1P3 receptors bound with two nonselective agonists or two S1P1-selective agonists were developed based on the X-ray crystal structure of bovine rhodopsin. The models predict that S1P1 Leu276 and S1P3 Phe263 contribute to the S1P1/S1P3 selectivity of the two S1P1-selective agonists. These residues were subjected to site-directed mutagenesis. The wild-type and mutant S1P receptors were expressed in Chinese hamster ovary cells and examined for their abilities to bind to and be activated by agonists in vitro. The results indicate that the mutations have minimal effects on the activities of the two nonselective agonists, although they have dramatic effects on the S1P1-selective agonists. These studies provide a fundamental understanding of how these two receptor-selective agonists bind to the S1P1 and S1P3 receptors, which should aid development of more selective S1P1 receptor agonists with immunosuppressive properties and improved safety profiles.

Versatility and Differential Roles of Cysteine Residues in Human Prostacyclin Receptor Structure and Function

Prostacyclin plays important roles in vascular homeostasis, promoting vasodilatation and inhibiting platelet thrombus formation. Previous studies have shown that three of six cytoplasmic cysteines, particularly those within the C-terminal tail, serve as important lipidation sites and are differentially conjugated to palmitoyl and isoprenyl groups (Miggin, S. M., Lawler, O. A., and Kinsella, B. T. (2003) J. Biol. Chem. 278, 6947-6958). Here we report distinctive roles for extracellular- and transmembrane-located cysteine residues in human prostacyclin receptor structure-function. Within the extracellular domain, all cysteines (4 of 4) appear to be involved in disulfide bonding interactions (i.e. a highly conserved Cys-92-Cys-170 bond and a putative non-conserved Cys-5-Cys-165 bond), and within the transmembrane (TM) region there are several cysteines (3 of 8) that maintain critical hydrogen bonding interactions (Cys-118 (TMIII), Cys-251 (TMVI), and Cys-202 (TMV)). This study highlights the necessity of sulfhydryl (SH) groups in maintaining the structural integrity of the human prostacyclin receptor, as 7 of 12 extracellular and transmembrane cysteines studied were found to be differentially indispensable for receptor binding, activation, and/or trafficking. Moreover, these results also demonstrate the versatility and reactivity of these cysteine residues within different receptor environments, that is, extracellular (disulfide bonds), transmembrane (H-bonds), and cytoplasmic (lipid conjugation).

Molecular dynamics simulation of the human adenosine A3 receptor: agonist induced conformational changes of Trp243

The adenosine A(3) receptor together with rhodopsin belongs to Class A of the G-protein coupled receptors. As the crystal structure of bovine rhodopsin represents the dark (inactive) state of the receptor, the details of GPCR activation are still unknown. In this molecular dynamics study we investigate how the homology model of the human adenosine A(3) receptor responds to ligand exposure. To this end we placed the homology model in a POPC membrane model. After equilibrating for 13 ns an agonist (Cl-IB-MECA) and an inverse agonist (PSB-10) were placed inside the putative binding pocket. In the following 10 ns molecular dynamics simulation we observed a different behaviour of the side-chain torsions of Trp243(6.48), depending on the presence or absence of the agonist or inverse agonist. This conformational change of Trp243 correlates with the assumed influence of ligands on receptor activation. Other predicted conformational changes of the receptor could not be observed yet. So Trp243 may represent the first switch in receptor activation.

Mapping the Binding Site of Arginine Vasopressin to V1aand V1bVasopressin Receptors

Starting from the 2.8-A resolution x-ray structure of bovine rhodopsin, three-dimensional molecular models of the complexes between arginine vasopressin and two receptor subtypes (V1a, V1b) have been built. Amino acid sequence alignment and docking studies suggest that four key residues (1.35, 2.65, 4.61, and 5.35) fine tune the binding of vasopressin and related peptide agonists to both receptor subtypes. To validate these predictions, a series of single or double mutants were engineered at V1a and V1b receptor subtypes and tested for their binding and functional properties. Two negatively charged amino acids at positions 1.35 and 2.65 are key anchoring residues to the Arg8 residue of arginine vasopressin. Moreover, two amino acids (V(4.61) and P(5.35)) delineating a hydrophobic subsite at the human V1b receptor are responsible for the recognition of V1b selective peptide agonists. Last, one of the latter positions (5.35) is hypothesized to explain the pharmacological species differences between rat and human vasopressin receptors for a V1b peptide agonist. Altogether these refined three-dimensional models of V1a and V1b human receptors should enable the identification of further new selective V1a and V1b agonists as pharmacological but also therapeutic tools.

Differences in Binding Sites of Two Melatonin Receptors Help to Explain Their Selectivity to Some Melatonin Analogs: A Molecular Modeling Study

Numerous diseases have been linked to the malfunction of G-protein coupled receptors (GP-CRs). Their adequate treatment requires rational design of new high-affinity and high-selectivity drugs targeting these receptors. In this work, we report three-dimensional models of the human MT1 and MT2 melatonin receptors, members of the GPCR family. The models are based on the X-ray structure of bovine rhodopsin. The computational approach employs an original procedure for optimization of receptor-ligand structures. It includes rotation of one of the transmembrane α-helices around its axis with simultaneous assessment of quality of the resulting complexes according to a number of criteria we have developed for this purpose. The optimal geometry of the receptor-ligand binding is selected based on the analysis of complementarity of hydrophobic/hydrophilic properties between the ligand and its protein environment in the binding site. The elaborated “optimized” models are employed to explore the details of protein-ligand interactions for melatonin and a number of its analogs with known affinity to MT1 and MT2 receptors. The models permit rationalization of experimental data, including those that were not used in model building. The perspectives opened by the constructed models and by the optimization procedure in the design of new drugs are discussed.

Exploring the Molecular Mechanism for Color Distinction in Humans

We examine here the role of the red, green, and blue human opsin structures in modulating the absorption properties of 11-cis-retinal bonded to the protein via a protonated Schiff base (PSB). We built the three-dimensional structures of the human red, green, and blue opsins using homology modeling techniques with the crystal structure of bovine rhodopsin as the template. We then used quantum mechanics (QM) combined with molecular mechanics (MM) (denoted as QM/MM) techniques in conjunction with molecular dynamics to determine how the room temperature molecular structures of the three human color opsin proteins modulate the absorption frequency of the same bound 11-cis-retinal chromophore to account for the differences in the observed absorption spectra. We find that the conformational twisting of the 11-cis-retinal PSB plays an important role in the green to blue opsin shift, whereas the dipolar side chains in the binding pocket play a surprising role of red-shifting the blue opsin with respect to the green opsin, as a fine adjustment to the opsin shift. The dipolar side chains play a large role in the opsin shift from red to green.

Proposal of a New Binding Orientation for Non-Peptide AT1 Antagonists: Homology Modeling, Docking and Three-Dimensional Quantitative Structure−Activity Relationship Analysis

A three-dimensional model of the AT1 receptor was constructed by means of a homology modeling procedure, using the X-ray structure of bovine rhodopsin as the initial template and taking into account the available site-directed mutagenesis data. The docking of losartan and its active metabolite EXP3174, followed by 1 ns of molecular dynamics (MD) simulation inserted into the phospholipid bilayer, suggested a different binding orientation for these antagonists from those previously proposed. Furthermore, the docking of several non-peptide antagonists was used as an alignment tool for the development of a three-dimensional quantitative structure-activity relationship (3D-QSAR) model, and the good results confirmed our binding hypothesis and the reliability of the model.

Stereoselectivity of 8-OH-DPAT toward the serotonin 5-HT 1A receptor: Biochemical and molecular modeling study

The great majority of pharmacological investigations of 5-HT 1A receptors’ reactivity has been performed using racemic 8-OH-DPAT, therefore the biochemical as well as behavioral profiles of both 8-OH-DPAT enantiomers are not circumstantiated. In the biochemical study capability of racemic 8-OH-DPAT (0.05, 0.1 mg/kg s.c.) and its counterparts R-8-OH-DPAT (0.05, 0.1 mg/kg s.c.) and S-8-OH-DPAT (0.05, 0.1 mg/kg s.c.) to influence 5-HT synthesis rate in rats’ prefrontal cortex, hypothalamus, hippocampus and brainstem was evaluated by HPLC/ED technique. Biochemical results are supported by the exhaustive computational study of possible differences between R- and S-enantiomer toward the 5-HT 1A receptor. A reliable 3D model of the rat 5-HT 1A receptor was constructed from the amino acid sequence using the crystal structure of bovine rhodopsin as a structural template. The structure of the receptor model was validated through docking studies and molecular dynamics simulations that gave results consistent with experimental data. Docking studies and the dynamics of ligand–receptor complexes emphasized different profiles of both enantiomers at the molecular level. The results of both biochemical and computational studies confirmed that R-enantiomer in contrast to S-8-OH-DPAT acts as full and potent agonist, whilst racemic form may display similar pharmacological profile to R-8-OH-DPAT.

Characterization of the Residues in Helix 8 of the Human β1-Adrenergic Receptor That Are Involved in Coupling the Receptor to G Proteins

Several key amino acids within amphipathic helix 8 of the human beta1-adrenergic receptor (beta1-AR) were mutagenized to characterize their role in signaling by G protein-coupled receptors. Mutagenesis of phenylalanine at position 383 in the hydrophobic interface to histidine (F383H) prevented the biosynthesis of the receptor, indicating that the orientation of helix 8 is important for receptor biosynthesis. Mutagenesis of aspartic acid at position 382 in the hydrophilic interface to leucine (D382L) reduced the binding and uncoupled the receptor from G protein activation. Mutagenesis of the basic arginine residue at position 384 to glutamine (R384Q) or to glutamic acid (R384E) increased basal and agonist-stimulated adenylyl cyclase activities. R384Q and R384E displayed features associated with constitutively active receptors because inverse agonists markedly reduced their elevated basal adenylyl cyclase activities. Isoproterenol increased the phosphorylation and promoted the desensitization of the Gly389 or Arg389 allelic variants of the wild type beta1-AR but failed to produce these effects in R384Q and R384E, because these receptors were maximally phosphorylated and desensitized under basal conditions. In contrast to the membranous distribution of the wild type beta1-AR, R384Q and R384E were localized mostly within intracellular punctate structures. Inverse agonists restored the membranous distribution of R384Q and R384E, indicating that they recycled normally when their constitutive internalization was blocked by inverse agonists. These data combined with computer modeling of the putative three-dimensional organization of helix 8 indicated that the amphipathic character of helix 8 and side chain projections of Asp382 and Arg384 within the hydrophilic interface might serve as a tethering site for the G protein.

Muscarinic Receptors: A Comparative Analysis of Structural Features and Binding Modes through Homology Modelling and Molecular Docking

Three-dimensional models of the five human muscarinic receptors were obtained from their known sequences. Homology modelling based on the crystallographic structure of bovine rhodopsin yielded models compatible with known results from site-directed mutagenesis studies. The only exceptions were the cytoplasmic loop 3 (CL3) in the five receptors, and the large C-terminal domain in M(1). Here, homology modelling with other closely related proteins allowed to solve these gaps. A detailed comparative discussion of the five models is given. The second part of the work involved docking experiments with the physiological ligand acetylcholine, again yielding results entirely compatible with results from mutagenesis experiments. The study revealed analogies and differences between the five receptors in the residues, and interactions leading to the recognition and binding of acetylcholine.