Constitutive activation of the beta2 adrenergic receptor alters the orientation of its sixth membrane-spanning segment

Does the Lipid Bilayer Orchestrate Access and Binding of Ligands to Transmembrane Orthosteric/Allosteric Sites of G Protein-Coupled Receptors?

The ligand-binding sites of many G protein-coupled receptors (GPCRs) are situated around and deeply embedded within the central pocket formed by their seven transmembrane-spanning α-helical domains. Generally, these binding sites are assumed accessible to endogenous ligands from the aqueous phase. Recent advances in the structural biology of GPCRs, along with biophysical and computational studies, suggest that amphiphilic and lipophilic molecules may gain access to these receptors by first partitioning into the membrane and then reaching the binding site via lateral diffusion through the lipid bilayer. In addition, several crystal structures of class A and class B GPCRs bound to their ligands offer unprecedented details on the existence of lipid-facing allosteric binding sites outside the transmembrane helices that can only be reached via lipid pathways. The highly organized structure of the lipid bilayer may direct lipophilic or amphiphilic drugs to a specific depth within the bilayer, changing local concentration of the drug near the binding site and affecting its binding kinetics. Additionally, the constraints of the lipid bilayer, including its composition and biophysical properties, may play a critical role in "pre-organizing" ligand molecules in an optimal orientation and conformation to facilitate receptor binding. Despite its clear involvement in molecular recognition processes, the critical role of the membrane in binding ligands to lipid-exposed transmembrane binding sites remains poorly understood and warrants comprehensive investigation. Understanding the mechanistic basis of the structure-membrane interaction relationship of drugs will not only provide useful insights about receptor binding kinetics but will also enhance our ability to take advantage of the apparent membrane contributions when designing drugs that target transmembrane proteins with improved efficacy and safety. In this minireview, we summarize recent structural and computational studies on membrane contributions to binding processes, elucidating both lipid pathways of ligand access and binding mechanisms for several orthosteric and allosteric ligands of class A and class B GPCRs.

Methods for the Development of In Silico GPCR Models

The Reggio group has constructed computer models of the inactive and G-protein-activated states of the cannabinoid CB1 and CB2 receptors, as well as, several orphan receptors that recognize a subset of cannabinoid compounds, including GPR55 and GPR18. These models have been used to design ligands, mutations, and covalent labeling studies. The resultant second-generation models have been used to design ligands with improved affinity, efficacy, and subtype selectivity. Herein, we provide a guide for the development of GPCR models using the most recent orphan receptor studied in our lab, GPR3. GPR3 is an orphan receptor that belongs to the Class A family of G-protein-coupled receptors. It shares high sequence similarity with GPR6, GPR12, the lysophospholipid receptors, and the cannabinoid receptors. GPR3 is predominantly expressed in mammalian brain and oocytes and it is known as a Gαs-coupled receptor activated constitutively in cells. GPR3 represents a possible target for the treatment of different pathological conditions such as Alzheimer's disease, oocyte maturation, or neuropathic pain. However, the lack of potent and selective GPR3 ligands is delaying the exploitation of this promising therapeutic target. In this context, we aim to develop a homology model that helps us to elucidate the structural determinants governing ligand-receptor interactions at GPR3. In this chapter, we detail the methods and rationale behind the construction of the GPR3 active-and inactive-state models. These homology models will enable the rational design of novel ligands, which may serve as research tools for further understanding of the biological role of GPR3.

Investigation of the influence of external factors on the conformational dynamics of rhodopsin-like receptors by means of molecular dynamics simulation

The study reports about the influence of binding of orthosteric ligands on the conformational dynamics of β-2-adrenoreceptor. Using molecular dynamics (MD) simulation, we found that there was a little fraction of active states of the receptor in its apo (ligand-free) ensemble. Analysis of MD trajectories indicated that such spontaneous activation of the receptor is accompanied by the motion in intracellular part of its alpha-helices. Thus, receptor's constitutive activity directly results from its conformational dynamics. On the other hand, the binding of a full agonist resulted in a significant shift of the initial equilibrium towards its active state. Finally, the binding of the inverse agonist stabilized the receptor in its inactive state. It is likely that the binding of inverse agonists might be a universal way of constitutive activity inhibition in vivo. Our results indicate that ligand binding redistribute pre-existing conformational degrees of freedom (in accordance to the Monod-Wyman-Changeux Model) of the receptor rather than cause induced fit in it. Therefore, the ensemble of biologically relevant receptor conformations is encoded in its spatial structure, and individual conformations from that ensemble might be used by the cell in conformity with the physiological behavior.

Crucial positively charged residues for ligand activation of the GPR35 receptor

GPR35 is a G protein-coupled receptor expressed in the immune, gastrointestinal, and nervous systems in gastric carcinomas and is implicated in heart failure and pain perception. We investigated residues in GPR35 responsible for ligand activation and the receptor structure in the active state. GPR35 contains numerous positively charged amino acids that face into the binding pocket that cluster in two distinct receptor regions, TMH3-4-5-6 and TMH1-2-7. Computer modeling implicated TMH3-4-5-6 for activation by the GPR35 agonists zaprinast and pamoic acid. Mutation results for the TMH1-2-7 region of GPR35 showed no change in ligand efficacies at the K1.32A, R2.65A, R7.33A, and K7.40A mutants. However, mutation of arginine residues in the TMH3-4-5-6 region (R4.60, R6.58, R3.36, R(164), and R(167) in the EC2 loop) had effects on signaling for one or both agonists tested. R4.60A resulted in a total ablation of agonist-induced activation in both the β-arrestin trafficking and ERK1/2 activation assays. R6.58A increased the potency of zaprinast 30-fold in the pERK assay. The R(167)A mutant decreased the potency of pamoic acid in the β-arrestin trafficking assay. The R(164)A and R(164)L mutants decreased potencies of both agonists. Similar trends for R6.58A and R(167)A were observed in calcium responses. Computer modeling showed that the R6.58A mutant has additional interactions with zaprinast. R3.36A did not express on the cell surface but was trapped in the cytoplasm. The lack of surface expression of R3.36A was rescued by a GPR35 antagonist, CID2745687. These results clearly show that R4.60, R(164), R(167), and R6.58 play crucial roles in the agonist initiated activation of GPR35.

Identification of the GPR55 antagonist binding site using a novel set of high-potency GPR55 selective ligands

GPR55 is a class A G protein-coupled receptor (GPCR) that has been implicated in inflammatory pain, neuropathic pain, metabolic disorder, bone development, and cancer. Initially deorphanized as a cannabinoid receptor, GPR55 has been shown to be activated by non-cannabinoid ligands such as l-α-lysophosphatidylinositol (LPI). While there is a growing body of evidence of physiological and pathophysiological roles for GPR55, the paucity of specific antagonists has limited its study. In collaboration with the Molecular Libraries Probe Production Centers Network initiative, we identified a series of GPR55 antagonists using a β-arrestin, high-throughput, high-content screen of ~300000 compounds. This screen yielded novel, GPR55 antagonist chemotypes with IC50 values in the range of 0.16-2.72 μM [Heynen-Genel, S., et al. (2010) Screening for Selective Ligands for GPR55: Antagonists (ML191, ML192, ML193) (Bookshelf ID NBK66153; PMID entry 22091481)]. Importantly, many of the GPR55 antagonists were completely selective, with no agonism or antagonism against GPR35, CB1, or CB2 up to 20 μM. Using a model of the GPR55 inactive state, we studied the binding of an antagonist series that emerged from this screen. These studies suggest that GPR55 antagonists possess a head region that occupies a horizontal binding pocket extending into the extracellular loop region, a central ligand portion that fits vertically in the receptor binding pocket and terminates with a pendant aromatic or heterocyclic ring that juts out. Both the region that extends extracellularly and the pendant ring are features associated with antagonism. Taken together, our results provide a set of design rules for the development of second-generation GPR55 selective antagonists.

Long range effect of mutations on specific conformational changes in the extracellular loop 2 of angiotensin II type 1 receptor

The topology of the second extracellular loop (ECL2) and its interaction with ligands is unique in each G protein-coupled receptor. When the orthosteric ligand pocket located in the transmembrane (TM) domain is occupied, ligand-specific conformational changes occur in the ECL2. In more than 90% of G protein-coupled receptors, ECL2 is tethered to the third TM helix via a disulfide bond. Therefore, understanding the extent to which the TM domain and ECL2 conformations are coupled is useful. To investigate this, we examined conformational changes in ECL2 of the angiotensin II type 1 receptor (AT1R) by introducing mutations in distant sites that alter the activation state equilibrium of the AT1R. Differential accessibility of reporter cysteines introduced at four conformation-sensitive sites in ECL2 of these mutants was measured. Binding of the agonist angiotensin II (AngII) and inverse agonist losartan in wild-type AT1R changed the accessibility of reporter cysteines, and the pattern was consistent with ligand-specific "lid" conformations of ECL2. Without agonist stimulation, the ECL2 in the gain of function mutant N111G assumed a lid conformation similar to AngII-bound wild-type AT1R. In the presence of inverse agonists, the conformation of ECL2 in the N111G mutant was similar to the inactive state of wild-type AT1R. In contrast, AngII did not induce a lid conformation in ECL2 in the loss of function D281A mutant, which is consistent with the reduced AngII binding affinity in this mutant. However, a lid conformation was induced by [Sar(1),Gln(2),Ile(8)] AngII, a specific analog that binds to the D281A mutant with better affinity than AngII. These results provide evidence for the emerging paradigm of domain coupling facilitated by long range interactions at distant sites on the same receptor.

Emerging opportunities for allosteric modulation of G-protein coupled receptors

Their ubiquitous nature, wide cellular distribution and versatile molecular recognition and signalling help make G-protein binding receptors (GPCRs) the most important class of membrane proteins in clinical medicine, accounting for ∼40% of all current therapeutics. A large percentage of current drugs target the endogenous ligand binding (orthosteric) site, which are structurally and evolutionarily conserved, particularly among members of the same GPCR subfamily. With the recent advances in GPCR X-ray crystallography, new opportunities for developing novel subtype selective drugs have emerged. Given the increasing recognition that the extracellular surface conformation changes in response to ligand binding, it is likely that all GPCRs possess an allosteric site(s) capable of regulating GPCR signalling. Allosteric sites are less structurally conserved than their corresponding orthosteric site and thus provide new opportunities for the development of more selective drugs. Constitutive oligomerisation (dimerisation) identified in many of the GPCRs investigated, adds another dimension to the structural and functional complexity of GPCRs. In this review, we compare 60 crystal structures of nine GPCR subtypes (rhodopsin, ß₂-AR, ß₁-AR, A(2a)-AR, CXCR4, D₃R, H₁R, M₂R, M₃R) across four subfamilies of Class A GPCRs, and discuss mechanisms involved in receptor activation and potential allosteric binding sites across the highly variable extracellular surface of these GPCRs. This analysis has identified a new extracellular salt bridge (ESB-2) that might be exploited in the design of allosteric modulators.

Investigation of D₂ receptor-agonist interactions using a combination of pharmacophore and receptor homology modeling

A combined modeling approach was used to identify structural factors that underlie the structure–activity relationships (SARs) of full dopamine D₂ receptor agonists and structurally similar inactive compounds. A 3D structural model of the dopamine D₂ receptor was constructed, with the agonist (−)-(R)-2-OH-NPA present in the binding site during the modeling procedure. The 3D model was evaluated and compared with our previously published D₂ agonist pharmacophore model. The comparison revealed an inconsistency between the projected hydrogen bonding feature (Ser-TM5) in the pharmacophore model and the TM5 region in the structure model. A new refined pharmacophore model was developed, guided by the shape of the binding site in the receptor model and with less emphasis on TM5 interactions. The combination of receptor and pharmacophore modeling also identified the importance of His393^6.55 for agonist binding. This convergent 3D pharmacophore and protein structure modeling strategy is considered to be general and can be highly useful in less well-characterized systems to explore ligand–receptor interactions. The strategy has the potential to identify weaknesses in the individual models and thereby provides an opportunity to improve the discriminating predictivity of both pharmacophore searches and structure-based virtual screens.

Identification of essential cannabinoid-binding domains: structural insights into early dynamic events in receptor activation

The classical cannabinoid agonist HU210, a structural analog of (-)-Δ(9)-tetrahydrocannabinol, binds to brain cannabinoid (CB1) receptors and activates signal transduction pathways. To date, an exact molecular description of the CB1 receptor is not yet available. Utilizing the minor binding pocket of the CB1 receptor as the primary ligand interaction site, we explored HU210 binding using lipid bilayer molecular dynamics (MD) simulations. Among the potential ligand contact residues, we identified residues Phe-174(2.61), Phe-177(2.64), Leu-193(3.29), and Met-363(6.55) as being critical for HU210 binding by mutational analysis. Using these residues to guide the simulations, we determined essential cannabinoid-binding domains in the CB1 receptor, including the highly sought after hydrophobic pocket important for the binding of the C3 alkyl chain of classical and nonclassical cannabinoids. Analyzing the simulations of the HU210-CB1 receptor complex, the CP55940-CB1 receptor complex, and the (-)-Δ(9)-tetrahydrocannabinol-CB1 receptor complex, we found that the positioning of the C3 alkyl chain and the aromatic stacking between Trp-356(6.48) and Trp-279(5.43) is crucial for the Trp-356(6.48) rotamer change toward receptor activation through the rigid-body movement of H6. The functional data for the mutant receptors demonstrated reductions in potency for G protein activation similar to the reductions seen in ligand binding affinity for HU210.

Identification of the GPR55 agonist binding site using a novel set of high-potency GPR55 selective ligands

Marijuana is the most widely abused illegal drug, and its spectrum of effects suggests that several receptors are responsible for the activity. Two cannabinoid receptor subtypes, CB1 and CB2, have been identified, but the complex pharmacological properties of exogenous cannabinoids and endocannabinoids are not fully explained by their signaling. The orphan receptor GPR55 binds a subset of CB1 and CB2 ligands and has been proposed as a cannabinoid receptor. This designation, however, is controversial as a result of recent studies in which lysophosphatidylinositol (LPI) was identified as a GPR55 agonist. Defining a biological role for GPR55 requires GPR55 selective ligands that have been unavailable. From a β-arrestin, high-throughput, high-content screen of 300000 compounds run in collaboration with the Molecular Libraries Probe Production Centers Network initiative (PubChem AID1965), we identified potent GPR55 selective agonists. By modeling of the GPR55 activated state, we compared the GPR55 binding conformations of three of the novel agonists obtained from the screen, CID1792197, CID1172084, and CID2440433 (PubChem Compound IDs), with that of LPI. Our modeling indicates the molecular shapes and electrostatic potential distributions of these agonists mimic those of LPI; the GPR55 binding site accommodates ligands that have inverted-L or T shapes with long, thin profiles that can fit vertically deep in the receptor binding pocket while their broad head regions occupy a horizontal binding pocket near the GPR55 extracellular loops. Our results will allow the optimization and design of second-generation GPR55 ligands and provide a means for distinguishing GPR55 selective ligands from those interacting with cannabinoid receptors.

Understanding functional residues of the cannabinoid CB1

The brain cannabinoid (CB(1)) receptor that mediates numerous physiological processes in response to marijuana and other psychoactive compounds is a G protein coupled receptor (GPCR) and shares common structural features with many rhodopsin class GPCRs. For the rational development of therapeutic agents targeting the CB(1) receptor, understanding of the ligand-specific CB(1) receptor interactions responsible for unique G protein signals is crucial. For a more than a decade, a combination of mutagenesis and computational modeling approaches has been successfully employed to study the ligand-specific CB(1) receptor interactions. In this review, after a brief discussion about recent advances in understanding of some structural and functional features of GPCRs commonly applicable to the CB(1) receptor, the CB(1) receptor functional residues reported from mutational studies are divided into three different types, ligand binding (B), receptor stabilization (S) and receptor activation (A) residues, to delineate the nature of the binding pockets of anandamide, CP55940, WIN55212-2 and SR141716A and to describe the molecular events of the ligand-specific CB(1) receptor activation from ligand binding to G protein signaling. Taken these CB(1) receptor functional residues, some of which are unique to the CB(1) receptor, together with the biophysical knowledge accumulated for the GPCR active state, it is possible to propose the early stages of the CB(1) receptor activation process that not only provide some insights into understanding molecular mechanisms of receptor activation but also are applicable for identifying new therapeutic agents by applying the validated structure-based approaches, such as virtual high throughput screening (HTS) and fragment-based approach (FBA).

Structural and functional roles of small group-conserved amino acids present on helix-H7 in the β(2)-adrenergic receptor

Sequence analysis of the class A G protein-coupled receptors (GPCRs) reveals that most of the highly conserved sites are located in the transmembrane helices. A second level of conservation exists involving those residues that are conserved as a group characterized by small and/or weakly polar side chains (Ala, Gly, Ser, Cys, Thr). These positions can have group conservation levels of up to 99% across the class A GPCRs and have been implicated in mediating helix-helix interactions in membrane proteins. We have previously shown that mutation of group-conserved residues present on transmembrane helices H2-H4 in the β(2)-adrenergic receptor (β(2)-AR) can influence both receptor expression and function. We now target the group-conserved sites, Gly315(7.42) and Ser319(7.46), on H7 for structure-function analysis. Replacing Ser319(7.46) with smaller amino acids (Ala or Gly) did not influence the ability of the mutant receptors to bind to the antagonist dihydroalprenolol (DHA) but resulted in ~15-20% agonist-independent activity. Replacement of Ser319(7.46) with the larger amino acid leucine lowered the expression of the S319L mutant and its ability to bind DHA. Both the G315A and G315S mutants also exhibited agonist-independent signaling, while the G315L mutant did not show specific binding to DHA. These data indicate that Gly315(7.42) and Ser319(7.46) are stabilizing β(2)-AR in an inactive conformation. We discuss our results in the context of van der Waals interactions of Gly315(7.42) with Trp286(6.48) and hydrogen bonding interactions of Ser319(7.46) with amino acids on H1-H2-H7 and with structural water.

GPCR Conformations: Implications for Rational Drug Design

G protein-coupled receptors (GPCRs) comprise a large class of transmembrane proteins that play critical roles in both normal physiology and pathophysiology. These critical roles offer targets for therapeutic intervention, as exemplified by the substantial fraction of current pharmaceutical agents that target members of this family. Tremendous contributions to our understanding of GPCR structure and dynamics have come from both indirect and direct structural characterization techniques. Key features of GPCR conformations derived from both types of characterization techniques are reviewed.

What site-directed labeling studies tell us about the mechanism of rhodopsin activation and G-protein binding

Rhodopsin is the photoreceptor protein responsible for dim-light vision in mammals. Due to extensive biophysical, structural and computational analysis of this membrane protein, it is presently the best understood G-protein coupled receptor (GPCR). Here I briefly review one approach that has been extensively used to identify dynamic and structural changes in rhodopsin--the use of site-directed labeling methods (SDL) coupled with electron paramagnetic resonance (EPR) and fluorescence spectroscopy. These SDL studies involve introducing individual cysteine residues into the receptor, then labeling them with cysteine-reactive probes for subsequent analysis by the appropriate spectroscopy. I will give a brief overview of how SDL methods are carried out and how the data is analyzed. Then, I will discuss how SDL studies were carried out on rhodopsin, and how they were used to identify a key structural change that occurs in rhodopsin upon activation--movement of transmembrane helix 6 (TM6). I will also briefly discuss how the SDL studies of rhodopsin compare with SDL studies of other GPCRs, and compare the SDL data with early and recent crystal structures of rhodopsin. Finally, I will discuss why the TM6 movement is required for rhodopsin to couple with the G-protein transducin, and speculate how this mechanism might be a universal method used by all GPCRs to bind G-proteins and perhaps other proteins involved in visual signal transduction, such as arrestin and rhodopsin kinase.

Interactions of the melanocortin-4 receptor with the peptide agonist NDP-MSH

Melanocortin-4 receptor (MC4R) has an important regulatory role in energy homeostasis and food intake. Peptide agonists of the MC4R are characterized by the conserved sequence His₆-Phe₇-Arg₈-Trp₉, which is crucial for their interaction with the receptor. This investigation utilized the covalent attachment approach to identify receptor residues in close proximity to the bound ligand [Nle₄,d-Phe₇]melanocyte-stimulating hormone (NDP-MSH), thereby differentiating between residues directly involved in ligand binding and those mutations that compromise ligand binding by inducing conformational changes in the receptor. Also, recent X-ray structures of G-protein-coupled receptors were utilized to refine a model of human MC4R in the active state (R^⁎), which was used to generate a better understanding of the binding mode of the ligand NDP-MSH at the atomic level.

The mutation of residues in the human MC4R—such as Leu106 of extracellular loop 1, and Asp122, Ile125, and Asp126 of transmembrane (TM) helix 3, His264 (TM6), and Met292 (TM7)—to Cys residues produced definitive indications of proximity to the side chains of residues in the core region of the peptide ligand. Of particular interest was the contact between d-Phe₇ on the ligand and Ile125 of TM3 on the MC4R. Additionally, Met292 (TM7) equivalent to Lys(7.45) (Ballesteros numbering scheme) involved in covalently attaching retinal in rhodopsin is shown to be in close proximity to Trp₉.

For the first time, the interactions between the terminal regions of NDP-MSH and the receptor are described. The amino-terminus appears to be adjacent to a series of hydrophilic residues with novel interactions at Cys196 (TM5) and Asp189 (extracellular loop 2). These interactions are reminiscent of sequential ligand binding exhibited by the β₂-adrenergic receptor, with the former interaction being equivalent to the known interaction involving Ser204 of the β₂-adrenergic receptor.

A lipid pathway for ligand binding is necessary for a cannabinoid G protein-coupled receptor

Recent isothiocyanate covalent labeling studies have suggested that a classical cannabinoid, (-)-7'-isothiocyanato-11-hydroxy-1',1'dimethylheptyl-hexahydrocannabinol (AM841), enters the cannabinoid CB2 receptor via the lipid bilayer (Pei, Y., Mercier, R. W., Anday, J. K., Thakur, G. A., Zvonok, A. M., Hurst, D., Reggio, P. H., Janero, D. R., and Makriyannis, A. (2008) Chem. Biol. 15, 1207-1219). However, the sequence of steps involved in such a lipid pathway entry has not yet been elucidated. Here, we test the hypothesis that the endogenous cannabinoid sn-2-arachidonoylglycerol (2-AG) attains access to the CB2 receptor via the lipid bilayer. To this end, we have employed microsecond time scale all-atom molecular dynamics (MD) simulations of the interaction of 2-AG with CB2 via a palmitoyl-oleoyl-phosphatidylcholine lipid bilayer. Results suggest the following: 1) 2-AG first partitions out of bulk lipid at the transmembrane alpha-helix (TMH) 6/7 interface; 2) 2-AG then enters the CB2 receptor binding pocket by passing between TMH6 and TMH7; 3) the entrance of the 2-AG headgroup into the CB2 binding pocket is sufficient to trigger breaking of the intracellular TMH3/6 ionic lock and the movement of the TMH6 intracellular end away from TMH3; and 4) subsequent to protonation at D3.49/D6.30, further 2-AG entry into the ligand binding pocket results in both a W6.48 toggle switch change and a large influx of water. To our knowledge, this is the first demonstration via unbiased molecular dynamics that a ligand can access the binding pocket of a class A G protein-coupled receptor via the lipid bilayer and the first demonstration via molecular dynamics of G protein-coupled receptor activation triggered by a ligand binding event.

Identification of orthosteric and allosteric site mutations in M2 muscarinic acetylcholine receptors that contribute to ligand-selective signaling bias

Muscarinic acetylcholine receptors contain at least one allosteric site that is topographically distinct from the acetylcholine, orthosteric binding site. Although studies have investigated the basis of allosteric modulation at these receptors, less is known about putative allosteric ligands that activate the receptor in their own right. We generated M(2) muscarinic acetylcholine receptor mutations in either the orthosteric site in transmembrane helices 3 and 6 (TM3 and -6) or part of an allosteric site involving the top of TM2, the second extracellular (E2) loop, and the top of TM7 and investigated their effects on the binding and function of the novel selective (putative allosteric) agonists (AC-42 (4-n-butyl-1-(4-(2-methylphenyl)-4-oxo-1-butyl)piperidine HCl), 77-LH-28-1 (1-(3-(4-butyl-1-piperidinyl)propyl)-3,3-dihydro-2(1H)-quinolinone), and N-desmethylclozapine) as well as the bitopic orthosteric/allosteric ligand, McN-A-343 (4-(m-chlorophenyl-carbamoyloxy)-2-butynyltrimethylammonium). Four classes of agonists were identified, depending on their response to the mutations, suggesting multiple, distinct modes of agonist-receptor interaction. Interestingly, with the exception of 77-LH-28-1, allosteric site mutations had no effect on the affinity of any of the agonists tested, but some mutations in the E2 loop influenced the efficacy of both orthosteric and novel selective agonists, highlighting a role for this region of the receptor in modulating activation status. Two point mutations (Y104(3.33)A (Ballesteros and Weinstein numbers in superscript) in the orthosteric and Y177A in the allosteric site) unmasked ligand-selective and signaling pathway-selective effects, providing evidence for the existence of pathway-specific receptor conformations. Molecular modeling of 77-LH-28-1 and N-desmethylclozapine yielded novel binding poses consistent with the possibility that the functional selectivity of such agents may arise from a bitopic mechanism.

Modulation of the constitutive activity of the ghrelin receptor by use of pharmacological tools and mutagenesis

Ghrelin and its receptor are important regulators of metabolic functions, including appetite, energy expenditure, fat accumulation, and growth hormone (GH) secretion. The ghrelin receptor is characterized by an ability to signal even without any ligand present with approximately 50% of the maximally ghrelin-induced efficacy-a feature that may have important physiological implications. The high basal signaling can be modulated either by administration of specific ligands or by engineering of mutations in the receptor structure. [D-Arg(1), D-Phe(5), D-Trp(7,9), Leu(11)]-substance P was the first inverse agonist to be identified for the ghrelin receptor, and this peptide has been used as a starting point for identification of the structural requirements for inverse agonist properties in the ligand. The receptor binding core motif was identified as D-Trp-Phe-D-Trp-Leu-Leu, and elongation of this peptide in the amino-terminal end determined the efficacy. Attachment of a positively charged amino acid was responsible for full inverse agonism, whereas an alanin converted the peptide into a partial agonist. Importantly, by use of mutational mapping of the residues critical for the modified D-Trp-Phe-D-Trp-Leu-Leu peptides, it was found that space-generating mutations in the deeper part of the receptor improved inverse agonism, whereas similar mutations located in the more extracellular part improved agonism. Modulation of the basal signaling by mutations in the receptor structure is primarily obtained by substitutions in an aromatic cluster that keep TMs VI and VII in close proximity to TM III and thus stabilize the active conformation. Also, substitution of a Phe in TM V is crucial for the high basal activity of the receptor as this residue serves as a partner for Trp VI:13 in the active conformation. It is suggested that inverse agonist and antagonist against the ghrelin receptor provide an interesting possibility in the development of drugs for treatment of obesity and diabetes and that improved structural understanding of the receptor function facilitates the drug development.

Conserved waters mediate structural and functional activation of family A (rhodopsin-like) G protein-coupled receptors

G protein-coupled receptors with seven transmembrane alpha-helices (GPCRs) comprise the largest receptor superfamily and are involved in detecting a wide variety of extracellular stimuli. The availability of high-resolution crystal structures of five prototypical GPCRs, bovine and squid rhodopsin, engineered A(2A)-adenosine, beta(1)- and beta(2)-adrenergic receptors, permits comparative analysis of features common to these and likely all GPCRs. We provide an analysis of the distribution of water molecules in the transmembrane region of these GPCR structures and find conserved contacts with microdomains demonstrated to be involved in receptor activation. Colocalization of water molecules associating with highly conserved and functionally important residues in several of these GPCR crystal structures supports the notion that these waters are likely to be as important to proper receptor function as the conserved residues. Moreover, in the absence of large conformational changes in rhodopsin after photoactivation, we propose that ordered waters contribute to the functional plasticity needed to transmit activation signals from the retinal-binding pocket to the cytoplasmic face of rhodopsin and that fundamental features of the mechanism of activation, involving these conserved waters, are shared by many if not all family A receptors.

Activation switch in the transmembrane domain of metabotropic glutamate receptor

Metabotropic glutamate receptors (mGluRs), members of family 3 G protein-coupled receptors, play pivotal roles in the modulation of synaptic transmission and are important drug targets for various neurological diseases. The structures of the extracellular ligand-binding domain (ECD) of mGluRs and its changes upon ligand binding have been well studied by various techniques, including X-ray crystallography. In contrast, little is known about the structure and structural changes of the seven-transmembrane domain (TMD). Here we searched for constitutively active mutation (CAM) sites in the TMD of mGluR8 to get insight into the epicenter of the structural changes in TMD, a potential target for allosteric ligands. Mutational analyses based on the knowledge of activating mutations of calcium-sensing receptor showed the presence of several CAM sites in the TMD of mGluR8. Among them, the site at position 764 in helix V, where threonine is present, is unique in that any substitution resulted in elevation of the basal activity, and some substitutions caused a loss of responsiveness to agonist. Further comprehensive mutational analyses indicated that the additional mutation of the CAM site at position 663 in helix III, where a tyrosine residue is present, generated a revertant phenotype. Mutations at these sites also increased the agonist binding affinity, although these sites are far from the ECD. These results indicated that the specific pair of amino acids at these CAM sites forms an activation switch that stabilizes the inactive state of mGluR8 and mediates the signal flow from the ECD to the cytoplasmic G protein-interacting site.

Ligand-binding architecture of human CB2 cannabinoid receptor: evidence for receptor subtype-specific binding motif and modeling GPCR activation

The extensive physiological influence of transmission through the CB2 cannabinoid receptor makes this G protein-coupled receptor (GPCR) a promising therapeutic target for treating neuropathic pain, inflammation, and immune disorders. However, there is little direct structural information pertaining to either GPCR or CB2-receptor ligand recognition and activation. The present work helps characterize experimentally the ligand-binding interactions of the human CB2 (hCB2) receptor. This study illustrates how our overall experimental approach, "ligand-assisted protein structure" (LAPS), affords direct determination of the requirements for ligand binding to the hCB2 receptor and discrimination among the binding motifs for ligands that activate therapeutically relevant GPCRs.

Residues accessible in the binding-site crevice of transmembrane helix 6 of the CB2 cannabinoid receptor

We have used the substituted-cysteine accessibility method (SCAM) to map the residues in the sixth membrane-spanning segment of the CB2 cannabinoid receptor that contribute to the surface of the water-accessible binding-site crevice. Using a background of the mutant C2.59S which is relatively insensitive to the methanethiosulfonate (MTS) reagents, we mutated to cysteine, one at a time, 34 consecutive residues in TMH6 of the CB2 receptor. These mutant receptors were then expressed in HEK293 cells. By incubating HEK293 cells stably transfected with CB2 receptors with the small, charged, hydrophilic, thiol-specific reagent methanethiosulfonate ethylammonium (MTSEA), [(3)H]CP55940 binding was significantly inhibited for six mutant receptors. All six of the mutants that reacted with MTSEA were protected from the reaction when pretreated with the cannabinoid agonist WIN55212-2, suggesting that MTSEA modification occurred within the binding crevice. Therefore, the side chains of the residues at these reactive loci (V6.51, L6.52, L6.54, M6.55, L6.59, and T6.62) are on the water-accessible surface of the binding-site crevice. These residues are extracellular to the TMH6 CWXP hinge motif. The pattern of accessibility is consistent with a alpha-helical conformation for this segment of TMH6. Molecular modeling studies performed in the context of the CB2 model show that V6.51, L6.52, L6.54, M6.55, L6.59, and T6.62 face into the CB2 binding pocket, further confirming our SCAM results. These results are similar to the accessibility patterns determined by SCAM studies of TMH6 in the opioid and dopamine D2 receptors.

Molecular modeling of adenosine receptors: new results and trends

Adenosine is a ubiquitous neuromodulator, which carries out its biological task by stimulating four cell surface receptors (A(1), A(2A), A(2B), and A(3)). Adenosine receptors (ARs) are members of the superfamily of G protein-coupled receptors (GPCRs). Their discovery opened up new avenues for potential drug treatment of a variety of conditions such as asthma, neurodegenerative disorders, chronic inflammatory diseases, and many other physiopathological states that are believed to be associated with changes in adenosine levels. Knowledge of the 3D structure of ARs could be of great help in the task of understanding their function and in the rational design of specific ligands. However, since GPCRs are membrane-bound proteins, high-resolution structural characterization is still an extremely difficult task. For this reason, great importance has been placed on molecular modeling studies and, particularly in the last few years, on homology modeling (HM) techniques. The publication of the first high-resolution crystal structure for bovine rhodopsin (bRh), a GPCR superfamily member, provides the option of utilizing HM to generate 3D models based on detailed structural information. In this review we report, analyze, and compare the main experimental data, computational HM procedures and validation methods used for ARs, describing in detail the most successful results.

Conformational change of the transmembrane helices II and IV of metabotropic glutamate receptor involved in G protein activation

G protein-coupled receptors are classified into several families on the basis of their amino acid sequences and the members of the same family exhibit sequence similarity but those of different families do not. In family 1 GPCRs such as rhodopsin and adrenergic receptor, extensive studies have revealed the stimulus-dependent conformational change of the receptor: the rearrangement of transmembrane helices III and VI is essential for G protein activation. In contrast, in family 3 GPCRs such as metabotropic glutamate receptor (mGluR), the inter-protomer relocation upon ligand binding has been observed but there is much less information about the structural changes of the transmsmbrane helices and the cytoplasmic domains. Here we identified constitutively active mutation sites at the cytoplasmic borders of helices II and IV of mGluR8 and successfully inhibited the G protein activation ability by engineering disulfide cross-linking between these cytoplasmic regions. The analysis of all possible single substitution mutants of these residues revealed that some steric interactions around these sites would be important to keep the receptor protein inactive. These results provided the model that the conformational changes at the cytoplasmic ends of helices II and IV of mGluR are involved in the efficient G protein coupling.

Ligand-stabilized conformational states of human beta(2) adrenergic receptor: insight into G-protein-coupled receptor activation

G-protein-coupled receptors (GPCRs) are known to exist in dynamic equilibrium between inactive- and several active-state conformations, even in the absence of a ligand. Recent experimental studies on the beta(2) adrenergic receptor (beta(2)AR) indicate that structurally different ligands with varying efficacies trigger distinct conformational changes and stabilize different receptor conformations. We have developed a computational method to study the ligand-induced rotational orientation changes in the transmembrane helices of GPCRs. This method involves a systematic spanning of the rotational orientation of the transmembrane helices (TMs) that are in the vicinity of the ligand for predicting the helical rotations that occur on ligand binding. The predicted ligand-stabilized receptor conformations are characterized by a simultaneous lowering of the ligand binding energy and a significant gain in interhelical and receptor-ligand hydrogen bonds. Using the beta(2)AR as a model, we show that the receptor conformational state depends on the structure and efficacy of the ligand for a given signaling pathway. We have studied the ligand-stabilized receptor conformations of five different ligands, a full agonist, norepinephrine; a partial agonist, salbutamol; a weak partial agonist, dopamine; a very weak agonist, catechol; and an inverse agonist, ICI-115881. The predicted ligand-stabilized receptor models correlate well with the experimentally observed conformational switches in beta(2)AR, namely, the breaking of the ionic lock between R131(3.50) at the intracellular end of TM3 (part of the DRY motif) and E268(6.30) on TM6, and the rotamer toggle switch on W286(6.48) on TM6. In agreement with trp-bimane quenching experiments, we found that norepinephrine and dopamine break the ionic lock and engage the rotamer toggle switch, whereas salbutamol, a noncatechol partial agonist only breaks the ionic lock, and the weak agonist catechol only engages the rotamer toggle switch. Norepinephrine and dopamine occupy the same binding region, between TM3, TM5, and TM6, whereas the binding site of salbutamol is shifted toward TM4. Catechol binds deeper into the protein cavity compared to the other ligands, making contact with TM5 and TM6. A part of the catechol binding site overlaps with those of dopamine and norepinephrine but not with that of salbutamol. Virtual ligand screening on 10,060 ligands on the norepinephrine-stabilized receptor conformation shows an enrichment of 38% compared to ligand unbound receptor conformation. These results show that ligand-induced conformational changes are important for developing functionally specific drugs that will stabilize a particular receptor conformation. These studies represent the first step toward a more universally applicable computational method for studying ligand efficacy and GPCR activation.

Computational modeling study of functional microdomains in cannabinoid receptor type 1

The seven transmembrane helices (TMH) G-protein-coupled receptors (GPCRs) constitute one of the largest superfamily of signaling proteins found in mammals. Some of its members, in which the cannabinoid (CB) receptors are included, stand out because their functional states can be modulated by a broad spectrum of effector molecules. The relative ligand promiscuity exhibited by these receptors could be related with particular attributes conferred by their molecular architecture and represents a motivating issue to be explored. In this regard, this study represents an effort to investigate the cannabinoid receptor type 1 (CB1) ligand recognition plasticity, using comparative modeling, molecular dynamics (MD) simulations and docking. Our results suggest that a cooperative set of subtle structural rearrangements within the TMHs provide to the CB1 protein the plasticity to reach alternate configurations. These changes include the relaxation of intramolecular constraints, the rotations, translations and kinks of the majority of TMHs and the reorganization of the ligand binding cavities.

Galpha-subunits differentially alter the conformation and agonist affinity of kappa-opioid receptors

Although ligand-induced conformational changes in G protein-coupled receptors (GPCRs) are well-documented, there is little direct evidence for G protein-induced changes in GPCR conformation. To investigate this possibility, the effects of overexpressing Galpha-subunits (Galpha16 or Galphai2) with the kappa-opioid receptor (KOR) were examined. The changes in KOR conformation were subequently examined via the substituted cysteine accessibility method (SCAM) in transmembrane domains 6 (TM6) and 7 (TM7) and extracellular loop 2 (EL2). Significant conformational changes were observed on TM7, the extracellular portion of TM6, and EL2. Seven SCAM-sensitive residues (S3107.33, F3147.37, and I3167.39 to Y3207.43) on TM7 presented a cluster pattern when the KOR was exposed to baseline amounts of G protein, and additional residues became sensitive upon overexpression of various G proteins. In TM7, S3117.34 and N3267.49 were found to be sensitive in Galpha16-overexpressed cells and Y3137.36, N3227.45, S3237.46, and L3297.52 in Galphai2-overexpressed cells. In addition, the degree of sensitivity for various TM7 residues was augmented, especially in Galphai2-overexpressed cells. A similar phenomenon was also observed for residues in TM6 and EL2. In addition to an enhanced sensitivity of certain residues, our findings also indicated that a slight rotation was predicted to occur in the upper part of TM7 upon G protein overexpression. These relatively modest conformational changes engendered by G protein overexpression had both profound and differential effects on the abilities of agonists to bind to KOR. These data are significant because they demonstrate that Galpha-subunits differentially modulate the conformation and agonist affinity of a prototypical GPCR.

Computational studies of Family A and Family B GPCRs

A full picture of the similarities between Family A and Family B GPCRs (G-protein coupled receptors) has been frustrated by the lack of clear homology between the respective sequences. Here, we review previous computational studies on GPCR dimerization in which the putative dimerization interfaces have been analysed using entropy, the ET (evolutionary trace) method and related methods. The results derived from multiple sequence alignments of Family A subfamilies have been mapped on to the rhodopsin crystal structure using standard alignments. Similarly, the results for the Family B alignments have been mapped on to the rhodopsin crystal structure using the 'cold-spot' alignment. For both Family A and Family B GPCRs, the sequence analysis indicates that there are functional sites on essentially all transmembrane helices, consistent with the parallel daisy chain model of GPCR oligomerization in which each GPCR makes interactions with a number of neighbouring GPCRs. The results are not too sensitive to the quality of the alignment. Molecular Dynamics simulations of the activation process within a single transmembrane bundle of the rhodopsin and the beta(2)-adrenergic receptor have been reviewed; the key observation, which is consistent with other computational studies, is that there is a translation and bending of helix 6, which contributes to a significant opening out of the intracellular face of the receptor, as shown in the accompanying movies. The simulations required the application of specific experiment-derived harmonic and half-harmonic distance restraints and so the application of such simulations to Family B GPCRs requires considerable care because of the alignment problem. Thus, in order to address the alignment problem, we have exploited the observation that GCR1, a plant GPCR, has homology with Family A, Family B and Family E GPCRs. The resulting alignment for transmembrane helix 3 is presented.

Distinct Structural Changes in a G Protein-coupled Receptor Caused by Different Classes of Agonist Ligands

The activity of G protein-coupled receptors can be modulated by different classes of ligands, including agonists that promote receptor signaling and inverse agonists that reduce basal receptor activity. The conformational changes in receptor structure induced by different agonist ligands are not well understood at present. In this study, we employed an in situ disulfide cross-linking strategy to monitor ligand-induced conformational changes in a series of cysteine-substituted mutant M(3) muscarinic acetylcholine receptors. The observed disulfide cross-linking patterns indicated that muscarinic agonists trigger a separation of the N-terminal segment of the cytoplasmic tail (helix 8) from the cytoplasmic end of transmembrane domain I. In contrast, inverse muscarinic agonists were found to increase the proximity between these two receptor regions. These findings provide a structural basis for the opposing biological effects of muscarinic agonists and inverse agonists. This study also provides the first piece of direct structural information as to how the conformations induced by these two functionally different classes of ligands differ at the molecular level. Given the high degree of structural homology found among most G protein-coupled receptors, our findings should be of broad general relevance.

Substituted NDP-MSH peptides paired with mutant melanocortin-4 receptors demonstrate the role of transmembrane 6 in receptor activation

The melanocortin-4 receptor (MC4R) is involved in regulating energy homeostasis and is a potential therapeutic target for obesity and cachexia. Molecular interactions between peptide ligands and MC4R have been studied in detail. Less is known regarding the role of these interactions in the mechanism of MC4R activation. The aim of this study was to investigate the molecular mechanism of human MC4R activation by [Nle4, d-Phe7]alpha-melanocyte-stimulating hormone (NDP-MSH), by first defining the role of the His6-d-Phe7-Arg8-Trp9 residues in receptor activation (Emax for stimulation of cAMP accumulation) using modified peptides, then understanding how their interaction with the receptor modulates activation using site-directed mutagenesis and a molecular model of NDP-MSH bound to the active state of the receptor. Alanine substitution indicated that the d-Phe7, Arg8, and Trp9 side chains contribute binding energy but are not essential for the receptor activation event. Conversely, His6 to Ala6 substitution reduced receptor activation but did not affect affinity. Chlorine substitutions on the d-Phe7 side chain also inhibited receptor activation. F261(6.51)A and F284(7.35)A receptor mutations acted as gain-of-function mutations, restoring efficacy to the His6 and d-Phe7 substituted peptides that had lost efficacy at the wild-type receptor. Based on a model of NDP-MSH and MC4R interaction, the antagonist behavior of these peptides is consistent with the prevention of transmembrane 6 (TM6) rotation. This data supports the hypothesis that increasing the size of d-Phe7 directly interferes with TM6 rotation, preventing receptor activation. We further propose that removing the interaction with the His6 side chain reorients the peptide within the binding pocket, indirectly impeding TM6 rotation by strengthening peptide interaction with F261(6.51) and F284(7.35). These findings refine the molecular basis for the mechanism of ligand-stimulated hMC4R activation and will be useful for the development of hMC4R agonists and antagonists.

The Angiotensin II AT1 Receptor Structure-Activity Correlations in the Light of Rhodopsin Structure

The most prevalent physiological effects of ANG II, the main product of the renin-angiotensin system, are mediated by the AT1 receptor, a rhodopsin-like AGPCR. Numerous studies of the cardiovascular effects of synthetic peptide analogs allowed a detailed mapping of ANG II's structural requirements for receptor binding and activation, which were complemented by site-directed mutagenesis studies on the AT1 receptor to investigate the role of its structure in ligand binding, signal transduction, phosphorylation, binding to arrestins, internalization, desensitization, tachyphylaxis, and other properties. The knowledge of the high-resolution structure of rhodopsin allowed homology modeling of the AT1 receptor. The models thus built and mutagenesis data indicate that physiological (agonist binding) or constitutive (mutated receptor) activation may involve different degrees of expansion of the receptor's central cavity. Residues in ANG II structure seem to control these conformational changes and to dictate the type of cytosolic event elicited during the activation. 1) Agonist aromatic residues (Phe8 and Tyr4) favor the coupling to G protein, and 2) absence of these residues can favor a mechanism leading directly to receptor internalization via phosphorylation by specific kinases of the receptor's COOH-terminal Ser and Thr residues, arrestin binding, and clathrin-dependent coated-pit vesicles. On the other hand, the NH2-terminal residues of the agonists ANG II and [Sar1]-ANG II were found to bind by two distinct modes to the AT1 receptor extracellular site flanked by the COOH-terminal segments of the EC-3 loop and the NH2-terminal domain. Since the [Sar1]-ligand is the most potent molecule to trigger tachyphylaxis in AT1 receptors, it was suggested that its corresponding binding mode might be associated with this special condition of receptors.

Mutation studies of Ser7.39 and Ser2.60 in the human CB1 cannabinoid receptor: evidence for a serine-induced bend in CB1 transmembrane helix 7

Ligands of structurally diverse natures are able to bind at the CB(1) cannabinoid receptor, suggesting the existence of multiple binding sites on the receptor. Modeling studies have implicated Ser2.60(173) and Ser7.39(383) as possible interaction site(s) for CB(1) agonists. To test the importance of these residues for receptor recognition, recombinant human CB(1) receptors, stably expressed in human embryonic kidney 293 cells, were used to investigate the consequences of mutating Ser2.60 (to S2.60A) or Ser7.39 (to S7.39A) in radioligand binding and guanosine 5'-3-O-(thio)triphosphate functional assays. The S7.39A mutant resulted in a total ablation of [(3)H](-)-3-[2-hydroxyl-4-(1,1-dimethylheptyl)phenyl]-4-[3-hydroxylpropyl] cyclohexan-1-ol (CP55,940) high-affinity binding. However, [(3)H](R)-(+)-[2,3-dihydro-5-methyl-3-[(4-morpholinyl)methyl]-pyrrolo[1,2,3-de]-1,4-benzoxazin-6-yl](1-naphthalenyl)methanone (WIN55,212-2) binding properties at S7.39A were comparable with those of the wild-type (WT) receptor. The binding affinity of (-)-11beta-hydroxy-3-(1',1'-dimethylheptyl)hexahydrocannabinol (AM4056) and (-)-11-hydroxydimethylheptyl-Delta(8)-tetrahydrocannabinol (HU210) were drastically reduced (50- to 100-fold) at the S7.39A mutant. Likewise, the EC(50) for HU210 and AM4056-mediated activation of the S7.39A receptor was increased by >200-fold. In contrast, the binding affinity and potency of WIN55,212-2, CP55,940, HU210, and AM4056 were unaltered at the S2.60A mutant compared with WT human CB(1) receptors. These results clearly suggest that Ser7.39, but not Ser2.60, plays a crucial role in mediating ligand specific interactions for CP55,940, HU210, and AM4056 at the human CB(1) receptor. Our modeling studies predict that Ser7.39 in a g-chi1 conformation may induce a helix bend in TMH7 that provides docking space for CP55,940 binding; the S7.39A mutation may alter this binding space, precluding CP55,940 binding.

Activation of G protein-coupled receptors

G protein-coupled receptors (GPCRs) mediate responses to hormones and neurotransmitters, as well as the senses of sight, smell, and taste. These remarkably versatile signaling molecules respond to structurally diverse ligands. Many GPCRs couple to multiple G protein subtypes, and several have been shown to activate G protein-independent signaling pathways. Drugs acting on GPCRs exhibit efficacy profiles that may differ for different signaling cascades. The functional plasticity exhibited by GPCRs can be attributed to structural flexibility and the existence of multiple ligand-specific conformational states. This chapter will review our current understanding of the mechanism by which agonists bind and activate GPCRs.

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.

G protein coupled receptor structure and activation

G protein coupled receptors (GPCRs) are remarkably versatile signaling molecules. The members of this large family of membrane proteins are activated by a spectrum of structurally diverse ligands, and have been shown to modulate the activity of different signaling pathways in a ligand specific manner. In this manuscript I will review what is known about the structure and mechanism of activation of GPCRs focusing primarily on two model systems, rhodopsin and the beta(2) adrenoceptor.

Cannabinoid CB1 receptor recognition of endocannabinoids via the lipid bilayer: molecular dynamics simulations of CB1 transmembrane helix 6 and anandamide in a phospholipid bilayer

The phospholipid bilayer plays a central role in the lifecycle of the endogenous cannabinoid, N-arachidonoylethanolamine (anandamide, AEA). Therefore, the orientation and location of AEA in the phospholipid bilayer with respect to key membrane associated proteins, is a central issue in understanding the mechanism of endocannabinoid signaling. In this paper, we report a test of the hypothesis that a betaXXbeta motif (formed by beta branching amino acids, V6.43 and I6.46) on the lipid face of the cannabinoid CB1 receptor in its inactive state may serve as an initial CB1 interaction site for AEA. Eight 6 ns NAMD2 molecular dynamics simulations of AEA were conducted in a model system composed of CB1 transmembrane helix 6 (TMH6) in a 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) bilayer. In addition, eight 6 ns NAMD2 molecular dynamics simulations of a low CB1 affinity (20:2, n-6) analog of AEA were conducted in the same model system. AEA was found to exhibit a higher incidence of V6.43/I6.46 groove insertion than did the (20:2, n-6) analog. In certain cases, AEA established a high energy of interaction with TMH6 by first associating with the V6.43/I6.46 groove and then molding itself to the lipid face of TMH6 to establish a hydrogen bonding interaction with the exposed backbone carbonyl of P6.50. Based upon these results, we propose that the formation of this hydrogen bonded AEA/TMH6 complex may be the initial step in CB1 recognition of AEA in the lipid bilayer.

Biarylpyrazole inverse agonists at the cannabinoid CB1 receptor: importance of the C-3 carboxamide oxygen/lysine3.28(192) interaction

The biarylpyrazole, N-(piperidin-1-yl)-5-(4-chlorophenyl)-1-(2,4-dichlorophenyl)-4-methyl-1H-pyrazole-3-carboxamide (SR141716; 1) has been shown to act as an inverse agonist/antagonist at the cannabinoid CB1 receptor. Our previous mutant cycle study suggested that K3.28(192) is involved in a direct interaction with the C-3 substituent of 1 in wild-type (WT) CB1.(1) However, these results did not establish what part of the C-3 substituent of 1 is involved in the K3.28(192) hydrogen bond, the carboxamide oxygen or the piperidine nitrogen. Furthermore, our previous calcium channel assay results for 5-(4- chlorophenyl)-3-[(E)-2-cyclohexylethenyl]-1-(2,4-dichlorophenyl)-4- methyl-1H-pyrazole (VCHSR; 2) (an analogue of 1 that lacks hydrogen-bonding capability in its C-3 substituent) showed that this compound acts as a neutral antagonist, a result that is in contrast to 1, which acts as an inverse agonist in this same assay.(1) These results suggested a relationship between biarylpyrazole interaction with K3.28(192) at CB1 and inverse agonism, but these results were for a single pair of compounds (1 and 2). The work presented here was designed to test two hypotheses derived from our modeling and mutant cycle results. The hypotheses are as follows: (1) it is the carboxamide oxygen of the C-3 substituent of 1 that interacts directly with K3.28(192) and (2) the interaction with K3.28(192) is crucial for the production of inverse agonism for biarylpyrazoles such as 1. To determine whether the carboxamide oxygen or the piperidine nitrogen of the C-3 substituent may be the interaction site for K3.28(192), we designed, synthesized, and evaluated a new set of analogues of 1 (3-6, Chart 1) in which modifications only to the C-3 substituent of 1 have been made. In each case, the modifications that were made preserved the geometry of this substituent in 1. The absence of the piperidine nitrogen was not found to affect affinity, whereas the absence of the carboxamide oxygen resulted in a reduction in affinity. CB1 docking studies in an inactive state model of CB1 resulted in the trend, 3,1<5,4<2<6 for ligand/CB1 interaction energies. This trend was consistent with the trend in WT CB1 Ki values versus [3H]CP55,940 reported here. In calcium channel assays, all analogues with carboxamide oxygens (1, 3, and 4) were found to be inverse agonists, whereas those that lacked this group (2, 5, and 6) were found to be neutral antagonists. Taken together, these results support the hypothesis that it is the carboxamide oxygen of the C-3 substituent of 1 that engages in a hydrogen bond with K3.28(192) in WT CB1. Furthermore, functional results for 1-6 support the hypothesis that the interaction of 1 with K3.28(192) may be key to its inverse agonism.

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

The great majority of pharmacological investigations of 5-HT1A 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-HT1A receptor. A reliable 3D model of the rat 5-HT1A 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.

Regulation of dopamine transporter trafficking by intracellular amphetamine

The dopamine (DA) transporter (DAT) mediates the removal of released DA. DAT is the major molecular target responsible for the rewarding properties and abuse potential of the psychostimulant amphetamine (AMPH). AMPH has been shown to reduce the number of DATs at the cell surface, and this AMPH-induced cell surface DAT redistribution may result in long-lasting changes in DA homeostasis. The molecular mechanism by which AMPH induces trafficking is not clear. Because AMPH is a substrate, we do not know whether extracellular AMPH stimulates trafficking through its interaction with DAT and subsequent alteration in DAT function, thereby triggering intracellular signaling or whether AMPH must be transported and then act intracellularly. In agreement with our previous studies, extracellular AMPH caused cytosolic redistribution of the wild-type human DAT (WT-hDAT). However, AMPH did not induce cytosolic redistribution in an uptake-impaired hDAT (Y335A-hDAT) that still binds AMPH. The divalent cation zinc (Zn(2+)) inhibits WT-hDAT activity, but it restores Y335A-hDAT uptake. Coadministration of Zn(2+) and AMPH consistently reduced WT-hDAT trafficking but stimulated cytosolic redistribution of Y335A-hDAT. Furthermore, direct intracellular application of AMPH, via a whole-cell patch pipette, stimulated the trafficking of Y335A-hDAT. Taken together, these data suggest that the DAT transport cycle is not required for AMPH-induced down-regulation and that an increase of intracellular AMPH is an essential component of DAT redistribution.

Mutations of CB1 T210 produce active and inactive receptor forms: correlations with ligand affinity, receptor stability, and cellular localization

Human cannabinoid receptor 1 (CB(1)) has attracted substantial interest as a potential therapeutic target for treating obesity and other obsessive disorders. An understanding of the mechanism governing the transition of the CB(1) receptor between its inactive and active states is critical for understanding how therapeutics can selectively regulate receptor activity. We have examined the importance of the Thr at position 210 in CB(1) in this transition, a residue predicted to be on the same face of the helix as the Arg of the DRY motif highly conserved in the G protein-coupled receptor superfamily. This Thr was substituted with Ile and Ala via mutagenesis, and the receptors, T210I and T210A, were expressed in HEK 293 cells. The T210I receptor exhibited enhanced agonist and diminished inverse agonist affinity relative to the wild type, consistent with a shift toward the active form. However, treatment with GTPgammaS to inhibit G protein coupling diminished the affinity change for the inverse agonist SR141716A. The decreased thermal stability of the T210I receptor and increased level of internalization of a T210I receptor-GFP chimera were also observed, consistent with constitutive activity. In contrast, the T210A receptor exhibited the opposite profile: diminished agonist and enhanced inverse agonist affinity. The T210A receptor was found to be more thermally stable than the wild type, and high levels of a T210A receptor-GFP chimera were localized to the cell surface as predicted for an inactive receptor form. These results suggest that T210 plays a key role in governing the transition between inactive and active CB(1) receptor states.

Metal Ion Site Engineering Indicates a Global Toggle Switch Model for Seven-transmembrane Receptor Activation

Much evidence indicates that, during activation of seven-transmembrane (7TM) receptors, the intracellular segments of the transmembrane helices (TMs) move apart with large amplitude, rigid body movements of especially TM-VI and TM-VII. In this study, AspIII:08 (Asp113), the anchor point for monoamine binding in TM-III, was used as the starting point to engineer activating metal ion sites between the extracellular segments of the beta2-adrenergic receptor. Cu(II) and Zn(II) alone and in complex with aromatic chelators acted as potent (EC50 decreased to 0.5 microm) and efficacious agonists in sites constructed between positions III:08 (Asp or His), VI:16 (preferentially Cys), and/or VII:06 (preferentially Cys). In molecular models built over the backbone conformation of the inactive rhodopsin structure, the heavy atoms that coordinate the metal ion were located too far away from each other to form high affinity metal ion sites in both the bidentate and potential tridentate settings. This indicates that the residues involved in the main ligand-binding pocket will have to move closer to each other during receptor activation. On the basis of the distance constraints from these activating metal ion sites, we propose a global toggle switch mechanism for 7TM receptor activation in which inward movement of the extracellular segments of especially TM-VI and, to some extent, TM-VII is coupled to the well established outward movement of the intracellular segments of these helices. We suggest that the pivots for these vertical seesaw movements are the highly conserved proline bends of the involved helices.

Molecular mechanism of 7TM receptor activation--a global toggle switch model

The multitude of chemically highly different agonists for 7TM receptors apparently do not share a common binding mode or active site but nevertheless act through induction of a common molecular activation mechanism. A global toggle switch model is proposed for this activation mechanism to reconcile the accumulated biophysical data supporting an outward rigid-body movement of the intracellular segments, as well as the recent data derived from activating metal ion sites and tethered ligands, which suggests an opposite, inward movement of the extracellular segments of the transmembrane helices. According to this model, a vertical see-saw movement of TM-VI-and to some degree TM-VII-around a pivot corresponding to the highly conserved prolines will occur during receptor activation, which may involve the outer segment of TM-V in an as yet unclear fashion. Small-molecule agonists can stabilize such a proposed active conformation, where the extracellular segments of TM-VI and -VII are bent inward toward TM-III, by acting as molecular glue deep in the main ligand-binding pocket between the helices, whereas larger agonists, peptides, and proteins can stabilize a similar active conformation by acting as Velcro at the extracellular ends of the helices and the connecting loops.