Mutational analysis and molecular modelling of the antagonist SR 144528 binding site on the human cannabinoid CB(2) receptor

Structural and functional insights into the G protein-coupled receptors: CB1 and CB2

The cannabinoid receptors CB1 and CB2 mediate a variety of physiological processes and continue to be explored as desirable drug targets. Both receptors are activated by the endogenous endocannabinoids and the psychoactive components of marijuana. Over the years, many efforts have been made to make selective ligands; however, the high degree of homology between cannabinoid receptor subtypes introduces challenges in studying either receptor in isolation. Recent advancements in structure biology have resulted in a surge of high-resolution structures, enriching our knowledge and understanding of receptor structure and function. In this review, of recent cannabinoid receptor structures, key features of the inactive and active state CB1 and CB2 are presented. These structures will provide additional insight into the modulation and signaling mechanism of cannabinoid receptors CB1 and CB2 and aid in the development of future therapeutics.

Design, Synthesis, and Biological Activity of New CB2 Receptor Ligands: from Orthosteric and Allosteric Modulators to Dualsteric/Bitopic Ligands

The design of dualsteric/bitopic agents as single chemical entities able to simultaneously interact with both the orthosteric and an allosteric binding site represents a novel approach in medicinal chemistry. Biased dualsteric/bitopic agents could enhance certain signaling pathways while diminishing the others that cause unwanted side effects. We have designed, synthesized, and functionally characterized the first CB2R heterobivalent bitopic ligands. In contrast to the parent orthosteric compound, our bitopic ligands selectively target CB2R versus CB1R and show a functional selectivity for the cAMP signaling pathway versus βarrestin2 recruitment. Moreover, the most promising bitopic ligand FD-22a displayed anti-inflammatory activity in a human microglial cell inflammatory model and antinociceptive activity in vivo in an experimental mouse model of neuropathic pain. Finally, computational studies clarified the binding mode of these compounds inside the CB2R, further confirming their bitopic nature.

Comprehensive review on the interaction between natural compounds and brain receptors: Benefits and toxicity

Given their therapeutic activity, natural products have been used in traditional medicines throughout the centuries. The growing interest of the scientific community in phytopharmaceuticals, and more recently in marine products, has resulted in a significant number of research efforts towards understanding their effect in the treatment of neurodegenerative diseases, such as Alzheimer's (AD), Parkinson (PD) and Huntington (HD). Several studies have shown that many of the primary and secondary metabolites of plants, marine organisms and others, have high affinities for various brain receptors and may play a crucial role in the treatment of diseases affecting the central nervous system (CNS) in mammalians. Actually, such compounds may act on the brain receptors either by agonism, antagonism, allosteric modulation or other type of activity aimed at enhancing a certain effect. The current manuscript comprehensively reviews the state of the art on the interactions between natural compounds and brain receptors. This information is of foremost importance when it is intended to investigate and develop cutting-edge drugs, more effective and with alternative mechanisms of action to the conventional drugs presently used for the treatment of neurodegenerative diseases. Thus, we reviewed the effect of 173 natural products on neurotransmitter receptors, diabetes related receptors, neurotrophic factor related receptors, immune system related receptors, oxidative stress related receptors, transcription factors regulating gene expression related receptors and blood-brain barrier receptors.

Human Cannabinoid Receptor 2 Ligand-Interaction Motif: Transmembrane Helix 2 Cysteine, C2.59(89), as Determinant of Classical Cannabinoid Agonist Activity and Binding Pose

Cannabinoid receptor 2 (CB2R)-dependent signaling is implicated in neuronal physiology and immune surveillance by brain microglia. Selective CB2R agonists hold therapeutic promise for inflammatory and other neurological disorders. Information on human CB2R (hCB2R) ligand-binding and functional domains is needed to inform the rational design and optimization of candidate druglike hCB2R agonists. Prior demonstration that hCB2R transmembrane helix 2 (TMH2) cysteine C2.59(89) reacts with small-molecule methanethiosulfonates showed that this cysteine residue is accessible to sulfhydryl derivatization reagents. We now report the design and application of two novel, pharmacologically active, high-affinity molecular probes, AM4073 and AM4099, as chemical reporters to interrogate directly the interaction of classical cannabinoid agonists with hCB2R cysteine residues. AM4073 has one electrophilic isothiocyanate (NCS) functionality at the C9 position of its cyclohexenyl C-ring, whereas AM4099 has NCS groups at that position and at the terminus of its aromatic A-ring C3 side chain. Pretreatment of wild-type hCB2R with either probe reduced subsequent [³H]CP55,940 specific binding by ∼60%. Conservative serine substitution of any hCB2R TMH cysteine residue except C2.59(89) did not affect the reduction of [³H]CP55,940 specific binding by either probe, suggesting that AM4073 and AM4099 interact irreversibly with this TMH2 cysteine. In contrast, AM841, an exceptionally potent hCB2R megagonist and direct AM4073/4099 congener bearing a single electrophilic NCS group at the terminus of its C3 side chain, had been demonstrated to bind covalently to TMH6 cysteine C6.47(257) and not C2.59(89). Molecular modeling indicates that the AM4073-hCB2R* interaction at C2.59(89) orients this classical cannabinoid away from TMH6 and toward the TMH2-TMH3 interface in the receptor's hydrophobic binding pocket, whereas the AM841-hCB2R* interaction at C6.47(257) favors agonist orientation toward TMH6/7. These data constitute initial evidence that TMH2 cysteine C2.59(89) is a component of the hCB2R binding pocket for classical cannabinoids. The results further demonstrate how interactions between classical cannabinoids and specific amino acids within the hCB2R* ligand-binding domain act as determinants of agonist pharmacological properties and the architecture of the agonist-hCB2R* conformational ensemble, allowing the receptor to adopt distinct activity states, such that interaction of classical cannabinoids with TMH6 cysteine C6.47(257) favors a binding pose more advantageous for agonist potency than does their interaction with TMH2 cysteine C2.59(89).

New pyridazinone-4-carboxamides as new cannabinoid receptor type-2 inverse agonists: Synthesis, pharmacological data and molecular docking

In the last few years, cannabinoid type-2 receptor (CB₂R) selective ligands have shown a great potential as novel therapeutic drugs in several diseases. With the aim of discovering new selective cannabinoid ligands, a series of pyridazinone-4-carboxamides was designed and synthesized, and the new derivatives tested for their affinity toward the hCB₁R and hCB₂R. The 6-(4-chloro-3-methylphenyl)-2-(4-fluorobenzyl)-N-(cis-4-methylcyclohexyl)-3-oxo-2,3-dihydropyridazine-4-carboxamide (9) displayed high CB₂-affinity (K_iCB₂ = 2.0 ± 0.81 nM) and a notable selectivity (K_iCB₁/K_iCB₂ > 2000). In addition, 9 and other active new synthesized entities have demonstrated to behave as CB₂R inverse agonists in [³⁵S]-GTPγS binding assay. ADME predictions of the newly synthesized CB₂R ligands suggest a favourable pharmacokinetic profile. Docking studies disclosed the specific pattern of interactions of these derivatives. Our results support that pyridazinone-4-carboxamides represent a new promising scaffold for the development of potent and selective CB₂R ligands.

Synthesis, molecular modeling and SAR study of novel pyrazolo[5,1-f][1,6]naphthyridines as CB2 receptor antagonists/inverse agonists

Pyrazolo[5,1-f][1,6]naphthyridine-carboxamide derivatives were synthesized and evaluated for the affinity at CB₁ and CB₂ receptors. Based on the AgOTf and proline-cocatalyzed multicomponent methodology, the ethyl 5-(p-tolyl)pyrazolo[5,1-f][1,6]naphthyridine-2-carboxylate (12) and ethyl 5-(2,4-dichlorophenyl)pyrazolo[5,1-f][1,6]naphthyridine-2-carboxylate (13) intermediates were synthesized from the appropriate o-alkynylaldehydes, p-toluenesulfonyl hydrazide and ethyl pyruvate. Most of the novel compounds feature a p-tolyl (8a-i) or a 2,4-dichlorophenyl (8j) motif at the C₅-position of the tricyclic pyrazolo[5,1-f][1,6]naphthyridine scaffold. Structural variation on the carboxamide moiety at the C₂-position includes basic monocyclic, terpenoid and adamantine-based amines. Among these derivatives, compound 8h (N-adamant-1-yl-5-(p-tolyl)pyrazolo[5,1-f][1,6]naphthyridine-2-carboxamide) exhibited the highest CB₂ receptor affinity (K_i=33nM) and a high degree of selectivity (K_iCB₁/K_iCB₂=173:1), whereas a similar trend in the near nM range was seen for the bornyl analogue (compound 8f, K_i=53nM) and the myrtanyl derivative 8j (K_i=67nM). Effects of 8h, 8f and 8j on forskolin-stimulated cAMP levels were determined, showing antagonist/inverse agonist properties for such compounds. Docking studies conducted for these derivatives and the reference antagonist/inverse agonist compound 4 (SR144528) disclosed the specific pattern of interactions probably related to the pyrazolo[5,1-f][1,6]naphthyridine scaffold as CB₂ inverse agonists.

Difference and Influence of Inactive and Active States of Cannabinoid Receptor Subtype CB2: From Conformation to Drug Discovery

Cannabinoid receptor 2 (CB2), a G protein-coupled receptor (GPCR), is a promising target for the treatment of neuropathic pain, osteoporosis, immune system, cancer, and drug abuse. The lack of an experimental three-dimensional CB2 structure has hindered not only the development of studies of conformational differences between the inactive and active CB2 but also the rational discovery of novel functional compounds targeting CB2. In this work, we constructed models of both inactive and active CB2 by homology modeling. Then we conducted two comparative 100 ns molecular dynamics (MD) simulations on the two systems-the active CB2 bound with both the agonist and G protein and the inactive CB2 bound with inverse agonist-to analyze the conformational difference of CB2 proteins and the key residues involved in molecular recognition. Our results showed that the inactive CB2 and the inverse agonist remained stable during the MD simulation. However, during the MD simulations, we observed dynamical details about the breakdown of the "ionic lock" between R131(3.50) and D240(6.30) as well as the outward/inward movements of transmembrane domains of the active CB2 that bind with G proteins and agonist (TM5, TM6, and TM7). All of these results are congruent with the experimental data and recent reports. Moreover, our results indicate that W258(6.48) in TM6 and residues in TM4 (V164(4.56)-L169(4.61)) contribute greatly to the binding of the agonist on the basis of the binding energy decomposition, while residues S180-F183 in extracellular loop 2 (ECL2) may be of importance in recognition of the inverse agonist. Furthermore, pharmacophore modeling and virtual screening were carried out for the inactive and active CB2 models in parallel. Among all 10 hits, two compounds exhibited novel scaffolds and can be used as novel chemical probes for future studies of CB2. Importantly, our studies show that the hits obtained from the inactive CB2 model mainly act as inverse agonist(s) or neutral antagonist(s) at low concentration. Moreover, the hit from the active CB2 model also behaves as a neutral antagonist at low concentration. Our studies provide new insight leading to a better understanding of the structural and conformational differences between two states of CB2 and illuminate the effects of structure on virtual screening and drug design.

Tricyclic pyrazoles. Part 8. Synthesis, biological evaluation and modelling of tricyclic pyrazole carboxamides as potential CB2 receptor ligands with antagonist/inverse agonist properties

Previous studies have investigated the relevance and structure-activity relationships (SARs) of pyrazole derivatives in relation with cannabinoid receptors, and the series of tricyclic 1,4-dihydroindeno[1,2-c]pyrazoles emerged as potent CB2 receptor ligands. In the present study, novel 1,4-dihydroindeno[1,2-c]pyrazole and 1H-benzo[g]indazole carboxamides containing a cyclopropyl or a cyclohexyl substituent were designed and synthesized to evaluate the influence of these structural modifications towards CB1 and CB2 receptor affinities. Among these derivatives, compound 15 (6-cyclopropyl-1-(2,4-dichlorophenyl)-N-(adamantan-1-yl)-1,4-dihydroindeno[1,2-c]pyrazole-3-carboxamide) showed the highest CB2 receptor affinity (Ki = 4 nM) and remarkable selectivity (KiCB1/KiCB2 = 2232), whereas a similar affinity, within the nM range, was seen for the fenchyl derivative (compound 10: Ki = 6 nM), for the bornyl analogue (compound 14: Ki = 38 nM) and, to a lesser extent, for the aminopiperidine derivative (compound 6: Ki = 69 nM). Compounds 10 and 14 were also highly selective for the CB2 receptor (KiCB1/KiCB2 > 1000), whereas compound 6 was relatively selective (KiCB1/KiCB2 = 27). The four compounds were also subjected to GTPγS binding analysis showing antagonist/inverse agonist properties (IC50 for compound 14 = 27 nM, for 15 = 51 nM, for 10 = 80 nM and for 6 = 294 nM), and this activity was confirmed for the three more active compounds in a CB2 receptor-specific in vitro bioassay consisting in the quantification of prostaglandin E2 release by LPS-stimulated BV2 cells, in the presence and absence of WIN55,212-2 and/or the investigated compounds. Modelling studies were also conducted with the four compounds, which conformed with the structural requirements stated for the binding of antagonist compounds to the human CB2 receptor.

Ligand-based homology modelling of the human CB2 receptor SR144528 antagonist binding site: a computational approach to explore the 1,5-diaryl pyrazole scaffold

SR144528 docking mode into the LBHM of the human CB₂ receptor antagonist binding site.

Modeling, molecular dynamics simulation, and mutation validation for structure of cannabinoid receptor 2 based on known crystal structures of GPCRs

The cannabinoid receptor 2 (CB2) plays an important role in the immune system. Although a few of GPCRs crystallographic structures have been reported, it is still challenging to obtain functional transmembrane proteins and high resolution X-ray crystal structures, such as for the CB2 receptor. In the present work, we used 10 reported crystal structures of GPCRs which had high sequence identities with CB2 to construct homology-based comparative CB2 models. We applied these 10 models to perform a prescreen by using a training set consisting of 20 CB2 active compounds and 980 compounds randomly selected from the National Cancer Institute (NCI) database. We then utilized the known 170 cannabinoid receptor 1 (CB1) or CB2 selective compounds for further validation. Based on the docking results, we selected one CB2 model (constructed by β1AR) that was most consistent with the known experimental data, revealing that the defined binding pocket in our CB2 model was well-correlated with the training and testing data studies. Importantly, we identified a potential allosteric binding pocket adjacent to the orthosteric ligand-binding site, which is similar to the reported allosteric pocket for sodium ion Na(+) in the A2AAR and the δ-opioid receptor. Our studies in correlation of our data with others suggested that sodium may reduce the binding affinities of endogenous agonists or its analogs to CB2. We performed a series of docking studies to compare the important residues in the binding pockets of CB2 with CB1, including antagonist, agonist, and our CB2 neutral compound (neutral antagonist) XIE35-1001. Then, we carried out 50 ns molecular dynamics (MD) simulations for the CB2 docked with SR144528 and CP55940, respectively. We found that the conformational changes of CB2 upon antagonist/agonist binding were congruent with recent reports of those for other GPCRs. Based on these results, we further examined one known residue, Val113(3.32), and predicted two new residues, Phe183 in ECL2 and Phe281(7.35), that were important for SR144528 and CP55940 binding to CB2. We then performed site-directed mutation experimental study for these residues and validated the predictions by radiometric binding affinity assay.

Examining the critical roles of human CB2 receptor residues Valine 3.32 (113) and Leucine 5.41 (192) in ligand recognition and downstream signaling activities

We performed molecular modeling and docking to predict a putative binding pocket and associated ligand-receptor interactions for human cannabinoid receptor 2 (CB2). Our data showed that two hydrophobic residues came in close contact with three structurally distinct CB2 ligands: CP-55,940, SR144528 and XIE95-26. Site-directed mutagenesis experiments and subsequent functional assays implicated the roles of Valine residue at position 3.32 (V113) and Leucine residue at position 5.41 (L192) in the ligand binding function and downstream signaling activities of the CB2 receptor. Four different point mutations were introduced to the wild type CB2 receptor: V113E, V113L, L192S and L192A. Our results showed that mutation of Val113 with a Glutamic acid and Leu192 with a Serine led to the complete loss of CB2 ligand binding as well as downstream signaling activities. Substitution of these residues with those that have similar hydrophobic side chains such as Leucine (V113L) and Alanine (L192A), however, allowed CB2 to retain both its ligand binding and signaling functions. Our modeling results validated by competition binding and site-directed mutagenesis experiments suggest that residues V113 and L192 play important roles in ligand binding and downstream signaling transduction of the CB2 receptor.

Mastering tricyclic ring systems for desirable functional cannabinoid activity

There is growing interest in using cannabinoid receptor 2 (CB2) agonists for the treatment of neuropathic pain and other indications. In continuation of our ongoing program aiming for the development of new small molecule cannabinoid ligands, we have synthesized a novel series of carbazole and γ-carboline derivatives. The affinities of the newly synthesized compounds were determined by a competitive radioligand displacement assay for human CB2 cannabinoid receptor and rat CB1 cannabinoid receptor. Functional activity and selectivity at human CB1 and CB2 receptors were characterized using receptor internalization and [(35)S]GTP-γ-S assays. The structure-activity relationship and optimization studies of the carbazole series have led to the discovery of a non-selective CB1 and CB2 agonist, compound 4. Our subsequent research efforts to increase CB2 selectivity of this lead compound have led to the discovery of CB2 selective compound 64, which robustly internalized CB2 receptors. Compound 64 had potent inhibitory effects on pain hypersensitivity in a rat model of neuropathic pain. Other potent and CB2 receptor-selective compounds, including compounds 63 and 68, and a selective CB1 agonist, compound 74 were also discovered. In addition, we identified the CB2 ligand 35 which failed to promote CB2 receptor internalization and inhibited compound CP55,940-induced CB2 internalization despite a high CB2 receptor affinity. The present study provides novel tricyclic series as a starting point for further investigations of CB2 pharmacology and pain treatment.

The importance of hydrogen bonding and aromatic stacking to the affinity and efficacy of cannabinoid receptor CB2 antagonist, 5-(4-chloro-3-methylphenyl)-1-[(4-methylphenyl)methyl]-N-[(1S,2S,4R)-1,3,3-trimethylbicyclo[2.2.1]hept-2-yl]-1H-pyrazole-3-carboxamide (SR144528)

Despite the therapeutic promise of the subnanomolar affinity cannabinoid CB2 antagonist, 5-(4-chloro-3-methylphenyl)-1-[(4-methylphenyl)methyl]-N-[(1S,2S,4R)-1,3,3-trimethylbicyclo[2.2.1]hept-2-yl]-1H-pyrazole-3-carboxamide (SR144528, 1), little is known about its binding site interactions and no primary interaction site for 1 at CB2 has been identified. We report here the results of Glide docking studies in our cannabinoid CB2 inactive state model that were then tested via compound synthesis, binding, and functional assays. Our results show that the amide functional group of 1 is critical to its CB2 affinity and efficacy and that aromatic stacking interactions in the TMH5/6 aromatic cluster of CB2 are also important. Molecular modifications that increased the positive electrostatic potential in the region between the fenchyl and aromatic rings led to more efficacious compounds. This result is consistent with the EC-3 loop negatively charged amino acid, D275 (identified via Glide docking studies) acting as the primary interaction site for 1 and its analogues.

Lead discovery, chemistry optimization, and biological evaluation studies of novel biamide derivatives as CB2 receptor inverse agonists and osteoclast inhibitors

N,N'-((4-(Dimethylamino)phenyl)methylene)bis(2-phenylacetamide) was discovered by using 3D pharmacophore database searches and was biologically confirmed as a new class of CB(2) inverse agonists. Subsequently, 52 derivatives were designed and synthesized through lead chemistry optimization by modifying the rings A-C and the core structure in further SAR studies. Five compounds were developed and also confirmed as CB(2) inverse agonists with the highest CB(2) binding affinity (CB(2)K(i) of 22-85 nM, EC(50) of 4-28 nM) and best selectivity (CB(1)/CB(2) of 235- to 909-fold). Furthermore, osteoclastogenesis bioassay indicated that PAM compounds showed great inhibition of osteoclast formation. Especially, compound 26 showed 72% inhibition activity even at the low concentration of 0.1 μM. The cytotoxicity assay suggested that the inhibition of PAM compounds on osteoclastogenesis did not result from its cytotoxicity. Therefore, these PAM derivatives could be used as potential leads for the development of a new type of antiosteoporosis agent.

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).

Mutagenesis and computer modeling studies of a GPCR conserved residue W5.43(194) in ligand recognition and signal transduction for CB2 receptor

W5.43(194), a conserved tryptophan residue among G-protein coupled receptors (GPCRs) and cannabinoid receptors (CB), was examined in the present report for its significance in CB2 receptor ligand binding and adenylyl cyclase (AC) activity. Computer modeling postulates that this site in CB2 may be involved in the affinity of WIN55212-2 and SR144528 through aromatic contacts. In the present study, we reported that a CB2 receptor mutant, W5.43(194)Y, which had a tyrosine (Y) substitution for tryptophan (W), retained the binding affinity for CB agonist CP55940, but reduced binding affinity for CB2 agonist WIN55212-2 and inverse agonist SR144528 by 8-fold and 5-fold, respectively; the CB2 W5.43(194)F and W5.43(194)A mutations significantly affect the binding activities of CP55940, WIN55212-2 and SR144528. Furthermore, we found that agonist-mediated inhibition of the forskolin-induced cAMP production was dramatically diminished in the CB2 mutant W5.43(194)Y, whereas W5.43(194)F and W5.43(194)A mutants resulted in complete elimination of downstream signaling, suggesting that W5.43(194) was essential for the full activation of CB2. These results indicate that both aromatic interaction and hydrogen bonding are involved in ligand binding for the residue W5.43(194), and the mutations of this tryptophan site may affect the conformation of the ligand binding pocket and therefore control the active conformation of the wild type CB2 receptor. W5.43(194)Y/F/A mutations also displayed noticeable enhancement of the constitutive activation probably attributed to the receptor conformational changes resulted from the mutations.

Pharmacological and molecular characterization of a dorsal root ganglion cell line expressing cannabinoid CB(1) and CB(2) receptors

The behavioral effects evoked by cannabinoids are primarily mediated by the CB(1) and CB(2) cannabinoid receptor subtypes. In vitro pharmacology of cannabinoid receptors has been elucidated using recombinant expression systems expressing either CB(1) or CB(2) receptors, with limited characterization in native cell lines endogenously expressing both CB(1) and CB(2) receptors. In the current study, we report the molecular and pharmacological characterization of the F-11 cell line, a hybridoma of rat dorsal root ganglion neurons and mouse neuroblastoma (N18TG2) cells, reported to endogenously express both cannabinoid receptors. The present study revealed that both receptors are of mouse origin in F-11 cells, and describes the relative gene expression levels between the two receptors. Pharmacological characterization of the F-11 cell line using cannabinoid agonists and antagonists indicated that the functional responses to these cannabinoid ligands are mainly mediated by CB(1) receptors. The non-selective cannabinoid ligands CP 55,940 and WIN 55212-2 are potent agonists and their efficacies in adenylate cyclase and MAPK assays are inhibited by the CB(1) selective antagonist SR141716A (SR1), but not by the CB(2) selective antagonist SR144528 (SR2). The endocannabinoid ligand 2AG, although not active in adenylate cyclase assays, was a potent activator of MAPK signaling in F-11 cells. The analysis of CB(1) and CB(2) receptor gene expression and the characterization of cannabinoid receptor pharmacology in the F-11 cell line demonstrate that it can be used as a tool for interrogating the endogenous signal transduction of cannabinoid receptor subtypes.

Modeling of ligand binding to G protein coupled receptors: cannabinoid CB1, CB2 and adrenergic β 2 AR

Cannabinoid and adrenergic receptors belong to the class A (similar to rhodopsin) G protein coupled receptors. Docking of agonists and antagonists to CB(1) and CB(2) cannabinoid receptors revealed the importance of a centrally located rotamer toggle switch and its possible participation in the mechanism of agonist/antagonist recognition. The switch is composed of two residues, F3.36 and W6.48, located on opposite transmembrane helices TM3 and TM6 in the central part of the membranous domain of cannabinoid receptors. The CB(1) and CB(2) receptor models were constructed based on the adenosine A(2A) receptor template. The two best scored conformations of each receptor were used for the docking procedure. In all poses (ligand-receptor conformations) characterized by the lowest ligand-receptor intermolecular energy and free energy of binding the ligand type matched the state of the rotamer toggle switch: antagonists maintained an inactive state of the switch, whereas agonists changed it. In case of agonists of β(2)AR, the (R,R) and (S,S) stereoisomers of fenoterol, the molecular dynamics simulations provided evidence of different binding modes while preserving the same average position of ligands in the binding site. The (S,S) isomer was much more labile in the binding site and only one stable hydrogen bond was created. Such dynamical binding modes may also be valid for ligands of cannabinoid receptors because of the hydrophobic nature of their ligand-receptor interactions. However, only very long molecular dynamics simulations could verify the validity of such binding modes and how they affect the process of activation.

hCB2 ligand-interaction landscape: cysteine residues critical to biarylpyrazole antagonist binding motif and receptor modulation

The human cannabinoid 2 GPCR (hCB2) is a prime therapeutic target. To define potential cysteine-related binding motifs critical to hCB2-ligand interaction, a library of hCB2 cysteine-substitution mutants and a novel, high-affinity biarylpyrazole hCB2 antagonist/inverse agonist (AM1336) functionalized to serve as a covalent affinity probe to target cysteine residues within (or in the microenvironment of) its hCB2 binding pocket were generated. The data provide direct experimental demonstration that both hCB2 TMH7 cysteines [i.e., C7.38(284) and C7.42(288)] are critical to optimal hCB2-AM1336 binding interaction and AM1336 pharmacological activity in a cell-based functional assay (cAMP formation). Elongating the AM1336 aliphatic side chain generated another novel hCB2 inverse agonist that binds covalently and selectively to C7.42(288) only. Identification of specific cysteine residues critical to hCB2 ligand interaction and function informs the structure-based design of hCB2-targeted medicines.

4-Oxo-1,4-dihydropyridines as selective CB2 cannabinoid receptor ligands: structural insights into the design of a novel inverse agonist series

Growing evidence shows that CB(2) receptor is an attractive therapeutic target. Starting from a series of 4-oxo-1,4-dihydroquinoline-3-carboxamide as selective CB(2) agonists, we describe here the medicinal chemistry approach leading to the development of CB(2) receptor inverse agonists with a 4-oxo-1,4-dihydropyridine scaffold. The compounds reported here show high affinity and potency at the CB(2) receptor while showing only modest affinity for the centrally expressed CB(1) cannabinoid receptor. Further, we found that the functionality of this series is controlled by its C-6 substituent because agonists bear a methyl or a tert-butyl group and inverse agonists, a phenyl or 4-chlorophenyl group, respectively. Finally, in silico studies suggest that the C-6 substituent could modulate the conformation of W6.48 known to be critical in GPCR activation.

2,3-Dihydro-1-benzofuran derivatives as a series of potent selective cannabinoid receptor 2 agonists: design, synthesis, and binding mode prediction through ligand-steered modeling

We recently discovered and reported a series of N-alkyl-isatin acylhydrazone derivatives that are potent cannabinoid receptor 2 (CB(2)) agonists. In an effort to improve the druglike properties of these compounds and to better understand and improve the treatment of neuropathic pain, we designed and synthesized a new series of 2,3-dihydro-1-benzofuran derivatives bearing an asymmetric carbon atom that behave as potent selective CB(2) agonists. We used a multidisciplinary medicinal chemistry approach with binding mode prediction through ligand-steered modeling. Enantiomer separation and configuration assignment were carried out for the racemic mixture for the most selective compound, MDA7 (compound 18). It appeared that the S enantiomer, compound MDA104 (compound 33), was the active enantiomer. Compounds MDA42 (compound 19) and MDA39 (compound 30) were the most potent at CB(2). MDA42 was tested in a model of neuropathic pain and exhibited activity in the same range as that of MDA7. Preliminary ADMET studies for MDA7 were performed and did not reveal any problems.

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.

Conformationally constrained analogues of N-(piperidinyl)-5-(4-chlorophenyl)-1-(2,4- dichlorophenyl)-4-methyl-1H-pyrazole-3-carboxamide (SR141716): design, synthesis, computational analysis, and biological evaluations

Structure-activity relationships (SARs) of 1 (SR141716) have been extensively documented, however, the conformational properties of this class have received less attention. In an attempt to better understand ligand conformations optimal for receptor recognition, we have designed and synthesized a number of derivatives of 1, including a four-carbon-bridged molecule (11), to constrain rotation of the diaryl rings. Computational analysis of 11 indicates approximately 20 kcal/mol energy barrier for rotation of the two aryl rings. NMR studies have determined the energy barrier to be approximately 18 kcal/mol and suggested atropisomers could exist. Receptor binding and functional studies with these compounds displayed reduced affinity and potency when compared to 1. This indicates that our structural modifications either constrain the ring systems in a suboptimal orientation for receptor interaction or the introduction of steric bulk leads to disfavored steric interactions with the receptor, and/or the relatively modest alterations in the molecular electrostatic potentials results in disfavored Coulombic interactions.

Binding mode prediction of conformationally restricted anandamide analogs within the CB1 receptor

BACKGROUND

CB1 cannabinoid receptors are G-protein coupled receptors for endocannabinoids including anandamide and 2-arachidonoylglycerol. Because these arachidonic acid metabolites possess a 20-carbon polyene chain as the alkyl terminal moiety, they are highly flexible with the potential to adopt multiple biologically relevant conformations, particularly those in a bent form. To better understand the molecular interactions associated with binding and steric trigger mechanisms of receptor activation, a series of conformationally-restricted anandamide analogs having a wide range of affinity and efficacy were evaluated.

RESULTS

A CB1 receptor model was constructed to include the extracellular loops, particularly extracellular loop 2 which possesses an internal disulfide linkage. Using both Glide (Schrödinger) and Affinity (Accelrys) docking programs, binding conformations of six anandamide analogs were identified that conform to rules applicable to the potent, efficacious and stereoselective non-classical cannabinoid CP55244. Calculated binding energies of the optimum structures from both procedures correlated well with the reported binding affinity values. The most potent and efficacious of the ligands adopted conformations characterized by interactions with both the helix-3 lysine and hydrophobic residues that interact with CP55244. The other five compounds formed fewer or less energetically favorable interactions with these critical residues. The flexibility of the tested anandamide analogs, measured by torsion angles around the benzene as well as the stretch between side chain moieties, could contribute to the differences in ability to interact with the CB1 receptor.

CONCLUSION

Analyses of multiple poses of conformationally-restricted anandamide analogs permitted identification of favored amino acid interactions within the CB1 receptor binding pocket. A ligand possessing both high affinity and cannabinoid agonist efficacy was able to interact with both polar and hydrophobic interaction sites utilized by the potent and efficacious non-classical cannabinoid CP55940. In contrast, other analogs characterized by reduced affinity or efficacy exhibited less favorable interactions with those key residues.

Targeting the cannabinoid CB2 receptor: modelling and structural determinants of CB2 selective ligands

Recent developments indicate that CB2 receptor ligands have the potential to become therapeutically important. To explore this potential, it is necessary to develop compounds with high affinity for the CB2 receptor and little affinity for the CB1 receptor. This review will discuss structure-activity relations at both receptors for classical cannabinoids and cannabimimetic indoles. Examples of CB2 selective ligands from both classes of compounds are presented and the structural features leading to selectivity are described. Two approaches, receptor mutations and molecular modelling, have been employed to investigate the interaction of ligands with both cannabinoid receptors. These results obtained from these techniques are discussed.

GPCR structure-based virtual screening approach for CB2 antagonist search

The potential for therapeutic specificity in regulating diseases has made cannabinoid (CB) receptors one of the most important G-protein-coupled receptor (GPCR) targets in search for new drugs. Considering the lack of related 3D experimental structures, we have established a structure-based virtual screening protocol to search for CB2 bioactive antagonists based on the 3D CB2 homology structure model. However, the existing homology-predicted 3D models often deviate from the native structure and therefore may incorrectly bias the in silico design. To overcome this problem, we have developed a 3D testing database query algorithm to examine the constructed 3D CB2 receptor structure model as well as the predicted binding pocket. In the present study, an antagonist-bound CB2 receptor complex model was initially generated using flexible docking simulation and then further optimized by molecular dynamic and mechanical (MD/MM) calculations. The refined 3D structural model of the CB2-ligand complex was then inspected by exploring the interactions between the receptor and ligands in order to predict the potential CB2 binding pocket for its antagonist. The ligand-receptor complex model and the predicted antagonist binding pockets were further processed and validated by FlexX-Pharm docking against a testing compound database that contains known antagonists. Furthermore, a consensus scoring (CScore) function algorithm was established to rank the binding interaction modes of a ligand on the CB2 receptor. Our results indicated that the known antagonists seeded in the testing database can be distinguished from a significant amount of randomly chosen molecules. Our studies demonstrated that the established GPCR structure-based virtual screening approach provided a new strategy with a high potential for in silico identifying novel CB2 antagonist leads based on the homology-generated 3D CB2 structure model.

A chemogenomic analysis of the transmembrane binding cavity of human G-protein-coupled receptors

The amino acid sequences of 369 human nonolfactory G-protein-coupled receptors (GPCRs) have been aligned at the seven transmembrane domain (TM) and used to extract the nature of 30 critical residues supposed--from the X-ray structure of bovine rhodopsin bound to retinal--to line the TM binding cavity of ground-state receptors. Interestingly, the clustering of human GPCRs from these 30 residues mirrors the recently described phylogenetic tree of full-sequence human GPCRs (Fredriksson et al., Mol Pharmacol 2003;63:1256-1272) with few exceptions. A TM cavity could be found for all investigated GPCRs with physicochemical properties matching that of their cognate ligands. The current approach allows a very fast comparison of most human GPCRs from the focused perspective of the predicted TM cavity and permits to easily detect key residues that drive ligand selectivity or promiscuity.

Cannabinoid CB2/CB1 selectivity. Receptor modeling and automated docking analysis

Three-dimensional models of the CB1 and CB2 cannabinoid receptors were constructed by means of a molecular 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 cannabinoid system was studied by means of docking techniques. An analysis of the interaction of WIN55212-2 with both receptors showed that CB2/CB1 selectivity is mainly determined by the interaction in the CB2 with the nonconserved residues S3.31 and F5.46, whose importance was suggested by site-directed mutagenesis data. We also carried out an automated docking of several ligands into the CB2 model, using the AUTODOCK 3.0 program; the good correlation obtained between the estimated free energy binding and the experimental binding data confirmed our binding hypothesis and the reliability of the model.

Alkylamides from Echinacea are a new class of cannabinomimetics. Cannabinoid type 2 receptor-dependent and -independent immunomodulatory effects

Alkylamides (alkamides) from Echinacea modulate tumor necrosis factor alpha mRNA expression in human monocytes/macrophages via the cannabinoid type 2 (CB2) receptor (Gertsch, J., Schoop, R., Kuenzle, U., and Suter, A. (2004) FEBS Lett. 577, 563-569). Here we show that the alkylamides dodeca-2E,4E,8Z,10Z-tetraenoic acid isobutylamide (A1) and dodeca-2E,4E-dienoic acid isobutylamide (A2) bind to the CB2 receptor more strongly than the endogenous cannabinoids. The Ki values of A1 and A2 (CB2 approximately 60 nM; CB1 >1500 nM) were determined by displacement of the synthetic high affinity cannabinoid ligand [3H]CP-55,940. Molecular modeling suggests that alkylamides bind in the solvent-accessible cavity in CB2, directed by H-bonding and pi-pi interactions. In a screen with 49 other pharmacologically relevant receptors, it could be shown that A1 and A2 specifically bind to CB2 and CB1. A1 and A2 elevated total intracellular Ca2+ in CB2-positive but not in CB2-negative promyelocytic HL60 cells, an effect that was inhibited by the CB2 antagonist SR144528. At 50 nM, A1, A2, and the endogenous cannabinoid anandamide (CB2 Ki >200 nM) up-regulated constitutive interleukin (IL)-6 expression in human whole blood in a seemingly CB2-dependent manner. A1, A2, anandamide, the CB2 antagonist SR144528 (Ki <10 nM), and also the non-CB2-binding alkylamide undeca-2E-ene,8,10-diynoic acid isobutylamide all significantly inhibited lipopolysaccharide-induced tumor necrosis factor alpha, IL-1beta, and IL-12p70 expression (5-500 nM) in a CB2-independent manner. Alkylamides and anandamide also showed weak differential effects on anti-CD3-versus anti-CD28-stimulated cytokine expression in human whole blood. Overall, alkylamides, anandamide, and SR144528 potently inhibited lipopolysaccharide-induced inflammation in human whole blood and exerted modulatory effects on cytokine expression, but these effects are not exclusively related to CB2 binding.

3D-QSAR Studies of Arylpyrazole Antagonists of Cannabinoid Receptor Subtypes CB1 and CB2. A Combined NMR and CoMFA Approach

The present work focuses on the study of the three-dimensional (3D) structural requirements for selective antagonist activity of arylpyrazole compounds at the cannabinoid CB1 and CB2 receptors. Initially, a combined high-resolution two-dimensional (2D) NMR and computer modeling approach was carried out to study the solution structure of the key pyrazole derivative N-(piperidin-1-yl)-5-phenyl-1-(n-pentyl)-4-methyl-1H-pyrazole-3-carboxamide (AM263). By using the NMR-determined molecular conformers as templates, the 3D quantitative structure-activity relationship (QSAR) studies were performed with the comparative molecular field analysis (CoMFA) approach on a set of arylpyrazole cannabinoid receptor antagonists. Molecular alignments suitable for deriving valuable pharmacophoric features for this series of compounds were determined. Such systematic 3D-QSAR/CoMFA analyses of 29 molecules and their receptor affinities gave guidance for understanding the binding affinities of arylpyrazoles at the CB1 and CB2 binding sites, respectively. Comparison of CoMFA steric and potential contour maps for affinity at the two cannabinoid receptor subtypes helps to differentiate structural requirements for each subtype and serves as a basis for the design of later-generation analogues.

Molecular biology of cannabinoid receptors

To date, two cannabinoid receptors have been isolated by molecular cloning. The CB1 and CB2 cannabinoid receptors are members of the G protein-coupled receptor family. There is also evidence for additional cannabinoid receptor subtypes. The CB1 and CB2 receptors recognize endogenous and exogenous cannabinoid compounds, which fall into five structurally diverse classes. Mutagenesis and molecular modeling studies have identified several key amino acid residues involved in the selective recognition of these ligands. Numerous residues involved in receptor activation have been elucidated. Regions of the CB1 receptor mediating desensitization and internalization have also been discovered. The known genetic structures of the CB1 and CB2 receptors indicate polymorphisms and multiple exons that maybe involved in tissue and species-specific regulation of these genes. The cannabinoid receptors are regulated during chronic agonist exposure, and gene expression is altered in disease states. There is a complex molecular architecture of the cannabinoid receptors that allows a single receptor to recognize multiple classes of compounds and produce an array of distinct downstream effects.

Chemical modification of the naphthoyl 3-position of JWH-015: in search of a fluorescent probe to the cannabinoid CB2 receptor

In silico modelling was used to guide the positioning of the fluorescent dye NBD-F on the cannabinoid CB2 receptor agonist JWH-015. While the ultimate fluorescent conjugate lost extensive binding affinity to the cannabinoid CB2 receptor, affinity and efficacy studies on the naphthoyl 3-position modified precursor molecules have provided new insight into structure-activity relationships associated with this position.

Cysteine residues in the human cannabinoid receptor: only C257 and C264 are required for a functional receptor, and steric bulk at C386 impairs antagonist SR141716A binding

The human neuronal cannabinoid receptor (CB1) is a G-protein-coupled receptor (GPCR) triggered by the psychoactive ingredients in marijuana, as well as endogenous cannabinoids produced in the brain. As with most GPCRs, the mechanism of CB1 activation is poorly understood. In this work, we have assessed the role of cysteine residues in CB1 ligand binding and activation, and demonstrate a method for mapping key determinants in CB1 structure and function. Through mutational analysis, we find that only two cysteines, C257 and C264, are required for high-level expression and receptor function. In addition, through cysteine reactivity studies, we find that a cysteine in transmembrane helix seven, C386 (C7.42), is reactive toward methanethiosulfonate (MTS) sulfhydryl labeling agents, and is thus solvent accessible. Interestingly, steric bulk introduced at this site, either through MTS labeling or by mutation, inhibits binding of the antagonist drug SR141716A (also known as Rimonabant or Accomplia), but does not affect the binding of the agonist CP55940. Our subsequent modeling studies suggest this effect is caused by steric clash of the modified C386 residue with the piperidine ring of SR141716A and/or disruption of an aromatic microdomain in the binding pocket. On the basis of these results, we hypothesize that bound SR141716A inhibits the ability of transmembrane helix 6 to move during formation of the functionally active receptor state.

Cysteine 2.59(89) in the Second Transmembrane Domain of Human CB2 Receptor Is Accessible within the Ligand Binding Crevice: Evidence for Possible CB2 Deviation from a Rhodopsin Template

In this study, the sensitivity of the CB2 receptor to methanethiosulfonate (MTS) derivatives was tested, and a native cysteine residue conferring the sensitivity was identified. By incubating human embryonic kidney 293 cells stably transfected with CB2 receptors and MTS derivatives such as MTS ethylammonium (MTSEA), [(3)H]HU-243 binding was inhibited. Pretreatment of the CB2 receptor with cannabinoid ligands prevented this inhibition, suggesting that MTSEA modification occurred within the binding crevice. To identify the cysteine(s) responsible for the MTSEA sensitivity, 10 CB2 mutants were prepared in which the eight cysteines in transmembrane domains or extracellular loop 2 were mutated to serine or alanine, one at a time or in combination. Five mutants exhibited specific [(3)H]HU-243 binding, with K(d) and B(max) values similar to those of wild-type CB2. However, five other mutants had no detectable ligand binding and were not detected on cell membranes by Western blot analysis. Among the five mutants with normal binding, only the sensitivity to MTSEA of the C2.59(89)S mutant was reduced significantly. These data demonstrate that C2.59(89) is the residue that mainly confers the inhibitory effect of MTSEA on ligand binding. Furthermore, the magnitude of the second-order rate constant (1.14 +/- 0.28 M(-1)s(-1)) for the MTSEA reaction with wild-type CB2 suggests that C2.59(89) resides at the margin of the CB2 binding site crevice. The accessibility of C2.59(89) to MTSEA provides experimental evidence for a possible conformational difference between TMH2 of CB2 versus Rho. Modeling studies undertaken to explore the origin of such differences suggest it is possibly caused by the conformational influence of S2.54(84).

Homology models of the cannabinoid CB1 and CB2 receptors. A docking analysis study

The 3D models of both CB1 and CB2 human receptors have been established by homology modeling using as template the X-ray structure of bovine Rhodopsin (code pdb: 1F88) a G-protein-coupled receptor (GPCR). A recursive approach comprising sequence alignment and model building was used to build both models, followed by the refinement of non-conserved regions. The cannabinoid system has been studied by means of docking techniques, using the 3D models of both CB1 and CB2 and well known reference inverse agonist/antagonist compounds. An approach based on the flexibility of the structures has been used to model the receptor-ligand complexes. The structural effects of ligand binding were studied and analyzed on the basis of hydrogen bond interactions, and binding energy calculations. Potential interaction sites of the receptor were determined from analysis of the difference accessible surface area (DASA) study of the protein with and without ligand.

NMR structural comparison of the cytoplasmic juxtamembrane domains of G-protein-coupled CB1 and CB2 receptors in membrane mimetic dodecylphosphocholine micelles

The fourth cytoplasmic domain, the so-called C-terminal juxtamembrane segment or helix VIII, has been identified in numerous G-protein-coupled receptors and exhibits unique functional characteristics. Efforts have been devoted to studying the juxtamembrane segment in order to understand the biological importance of the segment in G-protein activation of the cannabinoid CB1 and CB2 receptors. Recent biochemical data revealed that the CB1 C-terminal juxtamembrane peptide fragment CB1-(401-417) can directly activate the G-protein and also showed that the specificity of the signal transduction activation by the C-terminal juxtamembrane region is unique to the CB1 receptor but not to the CB2 receptor (Mukhopadhyay, S., and Howlett, A. C. (2001) Eur. J. Biochem. 268, 499-505). However, there is experimental work, not yet reported, on the conformational analyses and structural comparison between the respective helix VIII segments of the two receptors. In the present study, we have examined the conformational specificities of the cytoplasmic helical domains for both cannabinoid receptors. Three-dimensional structural features of two synthetic CB1 and CB2 peptides, CB1I397-G418 and CB2I298-K319, respectively, in membrane mimetic DPC micelles were studied using a combined high resolution NMR and computer modeling approach. Comparisons of the NMR-determined structures of the two peptides as well as their correspondent mutant peptides revealed their conformational properties and salt bridge dissimilarity, which might help us to understand the different structural roles of the fourth cytoplasmic helices in the function and regulation of CB1 and CB2 receptors.

Development of a 3D model for the human cannabinoid CB1 receptor

A novel comparison model of the human cannabinoid CB1 receptor has been constructed using the bovine rhodopsin X-ray structure as a template. The model was subjected to a 500-ps molecular dynamics simulation, and thereafter new conformers of the receptor model were produced in a simulated annealing procedure. Using an automated docking procedure, well-known cannabimimetic ligands were docked into six different model conformers, of which one was chosen for a detailed study of receptor-ligand interactions. The docking results confirm, for example, the importance of lysine K3.28(192) in the binding of these ligands. Also, other experimental data are fairly consistent with the present model, though there are some differences when compared to other recent CB1 comparison models. The present model will serve as a tool to investigate the receptor-ligand interactions and facilitate the design of novel cannabimimetic drugs.

Phylogenomic and chemotaxonomic analysis of the endocannabinoid system

The endocannabinoid system consists of two cannabinoid (CB) receptors, seven ligands, and ligand-catabolizing enzymes such as fatty acid amid hydrolase (FAAH) and monoglyceride lipase (MGL). The system's phylogenetic distribution is poorly known. The ligands cannot be molecularly investigated because they are not polypeptides and their specific synthetic enzymes have not been identified, so no sequences are available. Ligand phylogenetics can be inferred, nonetheless, by their presence in a range of extant organisms. Thus a meta-analysis of ligand extraction studies was performed (chemotaxonomy), and compared to a molecular search for homologs of CB receptors, vanilloid receptors (VR1), FAAH, and MGL in the genomes of sequenced organisms (phylogenomics). Putative homologs underwent functional mapping to ascertain the presence of critical amino acid motifs known to impart protein functionality. From an evolutionary perspective it appears that (1) endocannabinoid ligands evolved before CB receptors; (2) the ligands evolved independently multiple times; (3) CB receptors evolved prior to the metazoan-bilaterian divergence (ie, between extant Hydra and leech), but were secondarily lost in the Ecdysozoa; (4) VR1 may predate CB receptors but its affinity for endocannabinoids is a recent acquisition, appearing after the lower vertebrate-mammal divergence; (5) MGL may be as old as the ligands, whereas FAAH evolved recently, after the appearance of vertebrates. FAAH's emergence correlates with VR1's newly-found affinity for anandamide; this overlap in evolutionary time is recapitulated by complementary distribution patterns of FAAH, VR1, and anandamide in the brain. Linking FAAH, VR1, and anandamide implies a coupling among the remaining "older" parts of the endocannabinoid system, MGL, CB receptors, and 2-AG.

3D structural model of the G-protein-coupled cannabinoid CB2 receptor

The potential for therapeutic specificity in regulating diseases and for reduced side effects has made cannabinoid (CB) receptors one of the most important G-protein-coupled receptor (GPCR) targets for drug discovery. The cannabinoid (CB) receptor subtype CB2 is of particular interest due to its involvement in signal transduction in the immune system and its increased characterization by mutational and other studies. However, our understanding of their mode of action has been limited by the absence of an experimental receptor structure. In this study, we have developed a 3D model of the CB2 receptor based on the recent crystal structure of a related GPCR, bovine rhodopsin. The model was developed using multiple sequence alignment of homologous receptor sub-types in humans and mammals, and compared with other GPCRs. Alignments were analyzed with mutation scores, pairwise hydrophobicity profiles and Kyte-Doolittle plots. The 3D model of the transmembrane segment was generated by mapping the CB2 sequence onto the homologous residues of the rhodopsin structure. The extra- and intracellular loop regions of the CB2 were generated by searching for homologous C(alpha) backbone sequences in published structures in the Brookhaven Protein Databank (PDB). Residue side chains were positioned through a combination of rotamer library searches, simulated annealing and minimization. Intermediate models of the 7TM helix bundles were analyzed in terms of helix tilt angles, hydrogen-bond networks, conserved residues and motifs, possible disulfide bonds. The amphipathic cytoplasmic helix domain was also correlated with biological and site-directed mutagenesis data. Finally, the model receptor-binding cavity was characterized using solvent-accessible surface approach.

Functional mapping of cannabinoid receptor homologs in mammals, other vertebrates, and invertebrates

Over the past decade, several putative homologs of cannabinoid receptors (CBRs) have been identified by homology screening. Homology screening utilizes sequence alignment search engines to recognize homologs. We investigated these putative CBR homologs further by 'functional mapping' of their deduced amino acid sequences. The entire pharmacophore of a CBR has not yet been elucidated, but point-mutation studies have identified over 20 amino acid residues that impart CBR specificity for ligand recognition and/or signal transduction. Twenty point-mutation studies were used to construct a CBR functionality matrix. Sixteen putative CBR homologs were then mapped over the matrix. Several putative homologs did not hold up to this analysis: human GPR3, GPR6, GPR12, and Caenorhabditis elegans C02H7.2 expressed a series of crippling substitutions in the matrix, strongly suggesting they do not encode functional CBRs. Mapping the contested leech (Hirudo medicinalis) CBR sequence suggests that it encodes a functional CB1; it expresses fewer substitutions than the sea squirt (Ciona intestinalis) CB1 sequence. Mapping a putative CB2 ortholog in the puffer fish (Fugu rubripes T012234) suggests it may encode a CBR other than CB2. These findings are consistent with the lack of experimental data proving these putative CBRs have affinity for cannabinoid ligands. Matrix analysis also reveals that SR144528, a 'CB2-specific' synthetic antagonist, has affinity for non-mammalian CB1 receptors, and that L3.45 appears to be CB2-specific, its cognate in CB1 receptors is F3.45. In conclusion, functional mapping, utilizing point-mutation studies, may improve the specificity of homology screening performed by sequence alignment search engines.

Integrity of extracellular loop 1 of the human cannabinoid receptor 1 is critical for high-affinity binding of the ligand CP 55,940 but not SR 141716A

Like other G-protein coupled receptors with hydrophobic ligands, the human cannabinoid receptor 1 (CB1) is thought to bind its ligands within the transmembrane region of the receptor. However, for some of these receptors the extracellular loops (ECs) have also been shown to play a role in ligand recognition and selectivity. We have taken a mutagenesis approach to examine the role of the amino terminus, EC1, and EC3 of CB1 in ligand binding. Eight mutant receptors, each with a dipeptide insertion, were constructed, expressed, and evaluated for binding to the cannabinoid ligands (-)-cis-3[2-hydroxy-4-(1',1'-dimethylheptyl)phenyl]-trans-4-(3-hydroxypropyl)cyclohexanol (CP 55,940) and N-(piperidin-1-yl)-5-(4-chlorophenyl)-1-(2,4-dichlorophenyl)-4-methyl-1H-pyrazole-3-carboxamide hydrochloride (SR 141716A). Mutants with insertions in the membrane distal region of the amino terminus or EC3 maintained affinity for both ligands. Those with insertions in the membrane proximal region of the amino terminus or EC1 exhibited a loss of affinity for CP 55,940 while retaining wild-type affinity for SR 141716A. Representative mutants were analyzed for agonist-induced inhibition of cyclic AMP accumulation, and the results indicated that G-protein coupling remained intact. Alanine substitution mutants were made to address whether it was the perturbation of the overall structure of the region or the displacement of particular side chains that was responsible for the loss of CP 55,940 binding. We conclude that a structurally intact EC1, but not the comparably short EC3, is essential for high-affinity CP 55,940 binding.

A critical role for a tyrosine residue in the cannabinoid receptors for ligand recognition

Previous mutation and modeling studies have identified an aromatic cluster in the transmembrane helix (TMH) 3-4-5 region as important for ligand binding at the CB(1) and CB(2) cannabinoid receptors. Through novel mixed mode Monte Carlo/Stochastic Dynamics (MC/SD) calculations, we tested the importance of aromaticity at position 5.39(275) in CB(1). MC/SD calculations were performed on wild-type (WT) CB(1) and two mutants, Y5.39(275)F and Y5.39(275)I. Results indicated that while the CB(1) Y5.39(275)F mutant is very similar to WT, the Y5.39(275)I mutant shows pronounced topology changes in the TMH 3-4-5 region. Site-directed mutagenesis studies of tyrosine 5.39 to phenylalanine (Y-->F) or isoleucine (Y-->I) in both CB(1) and CB(2) were performed to determine the functional role of this amino acid in each receptor subtype. HEK 293 cells transfected with mutant receptor cDNAs were evaluated in radioligand binding and cyclic AMP assays. The CB(1) mutant and WT receptors were also co-expressed with G-protein-coupled inwardly rectifying channels (GIRK1 and GIRK4) in Xenopus oocytes to assess functional coupling. The Y-->F mutation resulted in cannnabinoid receptors with subtle differences in WT binding and signal transduction. In contrast, the Y-->I mutations produced receptors that could not produce signal transduction or bind to multiple cannabinoid compounds. However, immunofluorescence data indicate that the Y-->I mutation was compartmentalized and expressed at a level similar to that of the WT cannabinoid receptor. These results underscore the importance of aromaticity at position CB(1) 5.39(275) and CB(2) 5.39(191) for ligand recognition in the cannabinoid receptors.

Molecular biology of cannabinoid receptors

During the last decade, research on the molecular biology and genetics of cannabinoid receptors has led to a remarkable progress in understanding of the endogenous cannabinoid system, which functions in a plethora of physiological processes in the animal. At present, two types of cannabinoid receptors have been cloned from many vertebrates, and three endogenous ligands (the endocannabinoids arachidonoyl ethanolamide, 2-arachidonoyl glycerol and 2-arachidonoyl-glycerol ether) have been characterized. Cannabinoid receptor type 1 (CB(1)) is expressed predominantly in the central and peripheral nervous system, while cannabinoid receptor type 2 (CB(2)) is present almost exclusively in immune cells. Cannabinoid receptors have not yet been cloned from invertebrates, but binding proteins for endocannabinoids, endocannabinoids and metabolic enzyme activity have been described in a variety of invertebrates except for molting invertebrates such as Caenorhabditis elegans and Drosophila. In the central nervous system of mammals, there is strong evidence emerging that the CB(1) and its ligands comprise a neuromodulatory system functionally interacting with other neurotransmitter systems. Furthermore, the presynaptic localization of CB(1) together with the results obtained from electrophysiological experiments strengthen the notion that in cerebellum and hippocampus and possibly in other regions of the central nervous system, endocannabinoids may act as retrograde messengers to suppress neurotransmitter release at the presynaptic site. Many recent studies using genetically modified mouse lines which lack CB(1) and/or CB(2) finally could show the importance of cannabinoid receptors in animal physiology and will contribute to unravel the full complexity of the cannabinoid system.