3D modeling, ligand binding and activation studies of the cloned mouse delta, mu; and kappa opioid receptors

Refined 3D models of the transmembrane domains of the cloned delta, mu and kappa opioid receptors belonging to the superfamily of G-protein coupled receptors (GPCRs) were constructed from a multiple sequence alignment using the alpha carbon template of rhodopsin recently reported. Other key steps in the procedure were relaxation of the 3D helix bundle by unconstrained energy optimization and assessment of the stability of the structure by performing unconstrained molecular dynamics simulations of the energy optimized structure. The results were stable ligand-free models of the TM domains of the three opioid receptors. The ligand-free delta receptor was then used to develop a systematic and reliable procedure to identify and assess putative binding sites that would be suitable for similar investigation of the other two receptors and GPCRs in general. To this end, a non-selective, 'universal' antagonist, naltrexone, and agonist, etorphine, were used as probes. These ligands were first docked in all sites of the model delta opioid receptor which were sterically accessible and to which the protonated amine of the ligands could be anchored to a complementary proton-accepting residue. Using these criteria, nine ligand-receptor complexes with different binding pockets were identified and refined by energy minimization. The properties of all these possible ligand-substrate complexes were then examined for consistency with known experimental results of mutations in both opioid and other GPCRs. Using this procedure, the lowest energy agonist-receptor and antagonist-receptor complexes consistent with these experimental results were identified. These complexes were then used to probe the mechanism of receptor activation by identifying differences in receptor conformation between the agonist and the antagonist complex during unconstrained dynamics simulation. The results lent support to a possible activation mechanism of the mouse delta opioid receptor similar to that recently proposed for several other GPCRs. They also allowed the selection of candidate sites for future mutagenesis experiments.     

Desipramine modulation of sigma and opioid peptide receptor expression in glial cells.

Exposure of C6 glial cell cultures to desipramine induced the appearance of opioid receptors and up-regulated sigma receptors. Opioid binding was demonstrated with 3H-etorphine and 3H-dihydromorphine (DHM), but was not observed with the mu, delta and kappa ligands 3H-DAMGE, 3H-DADLE or 3H-(-)ethylketocyclazocine in the presence of specific blockers, respectively. Competition experiments with 3H-DHM and either (-)naloxone or (+)naloxone indicated the presence of authentic opioid receptors. In similar studies with beta-endorphin, its truncated form (1-27) or their N-acetyl derivatives, beta-endorphin proved to have the highest affinity. Opioid receptors in glial cell aggregates were primarily kappa, with few mu and delta sites. Desipramine increased Bmax values for kappa but not mu and delta.     

Molecular determinants of non-specific recognition of delta, mu, and kappa opioid receptors

Identification of the molecular determinants of recognition common to all three opioid receptors embedded in a single three-dimensional (3D) non-specific recognition pharmacophore has been carried out. The working hypothesis that underlies the computational study reported here is that ligands that bind with significant affinity to all three cloned opioid receptors, delta, mu, and kappa, but with different combinations of activation and inhibition properties at these receptors, could be promising behaviorally selective analgesics with diminished side effects. The study presented here represents the first step towards the rational design of such therapeutic agents. The common 3D pharmacophore developed for recognition of delta, mu, and kappa opioid receptors was based on the receptor affinities determined for 23 different opioid ligands that display no specificity for any of the receptor subtypes. The pharmacophore centers identified are a protonated amine, two hydrophobic groups, and the centroid of an aromatic group in a geometric arrangement common to all 23, non-specific, opioid ligands studied. Using this three-dimensional pharmacophore as a query for searching 3D structural databases, novel compounds potentially involved in non-specific recognition of delta, mu, and kappa opioid receptors were retrieved. These compounds can be valuable candidates for novel behaviorally selective analgesics with diminished or no side effects, and thus with potential therapeutic usefulness.     

Evidence for coupling of the kappa opioid receptor to brain GTPase

In membranes from guinea pig cerebellum, a tissue which predominantly contains kappa opioid receptors, exogenous and endogenous kappa-selective opioid agonists stimulated low-km GTPase activity by 11-20% with concentrations for half-maximal stimulation of 3-23 microM. Opioid ligands of the mu and delta type had no effect on GTPase in these membranes. Similar stimulation of GTPase by kappa opiates was obtained in rat and monkey brain membranes pretreated with beta-funaltrexamine (beta-FNA) and cis-(+/-)-3-methylfentanyl isothiocyanate (superfit) to alkylate the mu and delta receptors, respectively. The stimulation of brain GTPase by kappa opiates in both types of membranes was inhibited by naloxone with IC50's of 0.35 microM and 0.40 microM. The results demonstrate the coupling of the kappa opioid receptor to high affinity GTPase, the Ni regulatory protein of the adenylate cyclase complex.     

Pyrrolo- and pyridomorphinans: Non-selective opioid antagonists and delta opioid agonists/mu opioid partial agonists

Opioid ligands have found use in a number of therapeutic areas, including for the treatment of pain and opiate addiction (using agonists) and alcohol addiction (using antagonists such as naltrexone and nalmefene). The reaction of imines, derived from the opioid ligands oxymorphone and naltrexone, with Michael acceptors leads to pyridomorphinans with structures similar to known pyrrolo- and indolomorphinans. One of the synthesized compounds, 5e, derived from oxymorphone had substantial agonist activity at delta opioid receptors but not at mu and/or kappa opioid receptors and in that sense profiled as a selective delta opioid receptor agonist. The pyridomorphinans derived from naltrexone and naloxone were all found to be non-selective potent antagonists and as such could have utility as treatments for alcohol abuse.     

The occurrence of different kappa opioid receptors (K1 and K2) in frog (Rana esculenta) brain membranes

Opioid receptors of the frog (Rana esculenta) brain are characterized mainly by their relatively high ethylketocyclazocine (EKC) binding properties and by their low affinity to mu and delta ligands when D-Ala2-(Me)Phe4-Gly5-ol enkephalin (DAGO) and D-Ala2-Leu5-enkephalin (DALE) is used. In competition experiments it has been established that EKC and N-cyclopropylmethyl-norazidomorphine (CAM), which are non-selective kappa-ligands, have relatively high affinity to frog brain as well as the kappa 2 (which is DALE sensitive subpopulation of the kappa receptor) ligands etorphine and Metenkephalin-Arg6-Phe7 (1.). The kappa 2 subtype in frog brain resembles more to the mu subtype than to the delta subtype of opioid receptors, but it differs from the mu subtype in displaying low affinity toward beta-endorphin and DAGO.     

Homology modeling and molecular dynamics simulations of the active state of the nociceptin receptor reveal new insights into agonist binding and activation

The opioid receptor-like receptor, also known as the nociceptin receptor (NOP), is a class A G protein-coupled receptor (GPCR) in the opioid receptor family. Although NOP shares a significant homology with the other opioid receptors, it does not bind known opioid ligands and has been shown to have a distinct mechanism of activation compared to the closely related opioid receptors mu, delta, and kappa. Previously reported homology models of the NOP receptor, based on the inactive-state GPCR crystal structures, give limited information on the activation and selectivity features of this fourth member of the opioid receptor family. We report here the first active-state homology model of the NOP receptor based on the opsin GPCR crystal structure. An inactive-state homology model of NOP was also built using a multiple template approach. Molecular dynamics simulation of the active-state NOP model and comparison to the inactive-state model suggest that NOP activation involves movements of transmembrane (TM)3 and TM6 and several activation microswitches, consistent with GPCR activation. Docking of the selective nonpeptidic NOP agonist ligand Ro 64-6198 into the active-state model reveals active-site residues in NOP that play a role in the high selectivity of this ligand for NOP over the other opioid receptors. Docking the shortest active fragment of endogenous agonist nociceptin/orphaninFQ (residues 1-13) shows that the NOP extracellular loop 2 (EL2) loop interacts with the positively charged residues (8-13) of N/OFQ. Both agonists show extensive polar interactions with residues at the extracellular end of the TM domain and EL2 loop, suggesting agonist-induced reorganization of polar networks, during receptor activation.     

Morphine upregulates mu opioid receptors of human and monkey lymphocytes

Opioid receptors of subtypes delta, kappa, and mu similar to those found in brain cells have been identified in immune cells. The current study demonstrates by competitive polymerase chain reaction the treatment of human lymphocytic cells with morphine resulting in an increased amount of gene expression of mu opioid receptors. Antibodies against the MOR-1, the neuronal mu opioid receptor, were used in Western blot analysis of mu proteins and the results revealed a single band of approximately 50 kDa, the intensity of which was increased by morphine treatment. Similar results of mu opioid receptor activation were observed when monkey lymphocytes were treated with morphine. These studies suggest that in addition to causing an immune effect through communication with the neuroendocrine system, the psychoactive drug morphine may modulate immune functions by acting directly on the mu opioid receptors expressed on lymphocytes.     

Use of hydrophilic diamines for bridging of two opioid peptide pharmacophores. Synthesis and receptor binding of two new analogues.

The bivalent ligand approach, which assumes that two pharmacophores are connected by a spacer, was used to design receptor type-selective ligands for opioid receptors. The first two opioid peptide bivalent ligands with different spacer lengths containing different numbers of hydroxyl groups, (Tyr-D-Ala-Gly-Phe-NH-CH2-CHOH-)2 (Tyr-D-Ala-Gly-Phe-NH-CH2-CHOH-CHOH-)2, were synthesized and their binding to mu, delta, and kappa opioid receptors was characterized. Both analogues were found to possess high opioid in vitro activities. The length of the hydrophilic spacer does not affect the affinity for delta receptors, whereas shorter spacer length increases affinity for mu and even more so for kappa receptors. Thus receptor type-selective peptides for opioid receptors can be designed using the bivalent approach.     

Subunit Structure and Purification of Opioid Receptors

Abstract

Considerable evidence indicates the existence of multiple types of opioid receptors. The three major types have been named mu, delta and kappa. The earlier evidence was based on pharmacological as well as membrane binding experiments. This paper will emphasize more recent studies using solubilized opioid binding sites.

Several laboratories, including our own, have succeeded in separating kappa receptors from other types. A similar separation of mu from delta receptors has not yet been achieved. By crosslinking experiments with ¹²⁵I- human beta-endorphi we have been able to provide strong evidence for differences in molecular size between the major binding components of mu (65K) and delta (53K) receptors… It is not yet established whether the difference resides in the in the protein or carbohydrate portion of these glycoproteins. These results suggest that the three major types of opioid receptors represent distinct molecular entities.

An active opioid binding protein solubilized from bovine striatal membranes has been purified to apparent homogeneity. The major purification steps involve affinity chromatography and lectin chromatography on immobilized wheat germ agglutinin. The purified material gave a single band of molecular weight 65K Da on SDS-PAGE. Its specific activity for opioid binding was ca. 13000 pmol/mg protein and its properties are those of a component of the mu receptor.     

Synthesis of [D-Ala2,4'-125I-Phe3,Glu4]deltorphin and characterization of its delta opioid receptor binding properties.

Both [D-Ala2,Glu4]Deltorphin and [D-Ala2,4'-I-Phe3,Glu4]Deltorphin are highly selective ligands for delta, relative to mu, opioid receptors. Radiolabeled [D-Ala2, 4'-125I-Phe3,Glu4]Deltorphin ([125I]Deltorphin) was prepared with a specific activity of 2200 Ci/mmol from [D-Ala2, 4'-NH2-Phe3, Glu4]Deltorphin through a diazonium salt intermediate. The inhibition of [125I]Deltorphin binding to rat brain membranes by ligands selective for mu, delta, and kappa opioid receptors is consistent with binding by the radioligand to a single site having the properties of a delta opioid receptor. The results of these studies are in good agreement with those obtained by structurally different delta opioid receptor ligands. The similarity between the delta receptor site labeled by [125I]Deltorphin and those labeled by other delta receptor agonists, in contrast to differences seen by in vivo studies of their analgesic effects, is discussed.     

In-vitro investigation of oxazol and urea analogues of morphinan at opioid receptors

A series of 2-amino-oxazole (7 and 8) analogs and 2-one-oxazole analogs (9 and 10) were synthesized from cyclorphan (1) or butorphan (2) and evaluated in-vitro by their binding affinity at mu, delta, and kappa opioid receptors and compared with their 2-aminothiozole analogs 5 and 6. Ligands 7-10 showed decreased affinities at kappa and mu receptors. Urea analogs (11-14) were also prepared from 2-aminocyclorphan (3) or 2-aminobutorphan (4) and evaluated in-vitro by their binding affinity at mu, delta, and kappa opioid receptors. The urea derived opioids retained their affinities at mu receptors while showing increased affinities at delta receptors and decreased affinities at kappa receptors. Functional activities of these compounds were measured in the [35S]GTPgammaS binding assay, illustrating that all of these ligands were kappa agonists. At the mu receptor, compounds 11 and 12 were mu agonist/antagonists.     

Multiple templates-based homology modeling enhances structure quality of AT1 receptor: validation by molecular dynamics and antagonist docking

We present a comparative account on 3D-structures of human type-1 receptor (AT1) for angiotensin II (AngII), modeled using three different methodologies. AngII activates a wide spectrum of signaling responses via the AT1 receptor that mediates physiological control of blood pressure and diverse pathological actions in cardiovascular, renal, and other cell types. Availability of 3D-model of AT1 receptor would significantly enhance the development of new drugs for cardiovascular diseases. However, templates of AT1 receptor with low sequence similarity increase the complexity in straightforward homology modeling, and hence there is a need to evaluate different modeling methodologies in order to use the models for sensitive applications such as rational drug design. Three models were generated for AT1 receptor by, (1) homology modeling with bovine rhodopsin as template, (2) homology modeling with multiple templates and (3) threading using I-TASSER web server. Molecular dynamics (MD) simulation (15 ns) of models in explicit membrane-water system, Ramachandran plot analysis and molecular docking with antagonists led to the conclusion that multiple template-based homology modeling outweighs other methodologies for AT1 modeling.     

Exploring the unique pharmacology of a novel opioid receptor, ZFOR1, using molecular modeling and the 'message-address' concept.

Previous studies have probed the structural basis of ligand selectivity in the mu, delta and kappa opioid receptors through the application of molecular modeling techniques in concert with the 'message-address' concept. Here, this approach was used in an attempt to rationalize the unique pharmacological profile of a recently cloned novel opioid receptor, ZFOR1 (ZebraFish Opioid Receptor 1). Specifically, a model of the transmembrane domains of ZFOR1 was constructed and used to explore the binding modes of various prototypical opioid ligands. The results show that the 'message' portion of the binding pocket of ZFOR1 is highly conserved; hence, the binding modes of non-selective opioid ligands are well preserved. In contrast, a small number of variant residues at the extracellular end of the binding pocket, particularly Lys288 (VI:26) and Trp304 (VII:03), are shown to create adverse steric interactions with all delta and kappa selective ligands examined, thereby disrupting their binding modes. These results are consistent with, and serve as an explanation for, the observed pharmacology of this receptor, lending support to both the validity of the 'message-address' concept itself and to the use of molecular modeling approaches in its application.     

Conformation-activity relationships of opioid peptides with selective activities at opioid receptors

The discovery of endogenous opioid peptides 25 years ago opened up a new chapter in efforts to understand the origins and control of pain, its relationships to other biological functions, including inflammatory and other immune responses, and the relationships of opioid peptides and their receptors to a variety of undesirable or toxic side effects often associated with the nonpeptide opiates such as morphine including addiction, constipation, a variety of neural toxicities, tolerance, and respiratory depression. For these investigations the need for potent and highly receptor selective agonists and antagonists has been crucial since they in principle allow one to distinguish unequivocally the roles of the different opioid receptors (mu, delta, and kappa) in the various biological and pathological roles of the opioid peptides and their receptors. Conformational and topographical constraint of the linear natural endogenous opioid peptides has played a major role in developing peptide ligands with high selectivity for mu, delta, and kappa receptors, and in understanding the conformational, topographical, and stereoelectronic structural requirements of the opioid peptides for their interactions with opioid receptors. In turn, this had led to insights into the three-dimensional pharmacophore for opioid receptors. In this article we review and discuss some of the developments that have led to potent, selective, and stable peptide and peptidomimetic ligands that are highly potent and selective, and that have delta agonist, mu antagonist, and kappa agonist biological activities (other authors in this issue will discuss the development of other types of activities and selectivities). These have led to ligands that provide unique insight into opioid pharmacophores and the critical roles opioid ligands and receptor scan play in pain, addiction, and other human maladies.     

(4-Carboxamido)phenylalanine is a surrogate for tyrosine in opioid receptor peptide ligands

(S)-4-(Carboxamido)phenylalanine (Cpa) is examined as a bioisosteric replacement for the terminal tyrosine (Tyr) residue in a variety of known peptide ligands for the mu, delta and kappa opioid receptors. The Cpa-containing peptides, assayed against cloned human opioid receptors, display comparable binding affinity (Ki), and agonist potency (EC50) to the parent ligands at the three receptors. Cpa analogs of delta selective peptides show an increase in delta selectivity relative to the mu receptor. Cpa is the first example of an amino acid that acts as a surrogate for Tyr in opioid peptide ligands, challenging the long-standing belief that a phenolic residue is required for high affinity binding.     

Modeling and simulation of the human delta opioid receptor

A model for the human delta opioid receptor has been generated via sequence alignment, structure building using the crystal structure of bovine rhodopsin as a template, and refinement by molecular dynamics simulation. The model building suggested that, in addition to the previously postulated interaction between D128 and Y308, an internal salt bridge also exists between residues D128 and R192, both of which are conserved in all the opioid receptors. The model and salt bridge were then shown to be stable during a 20-nsec simulation in a lipid bilayer. It is therefore proposed that both of these interactions play a role in stabilizing the inactive state of the receptor. The model is also used in an effort to rationalize many of the mutational studies performed on delta opioid receptors, and to suggest a plausible explanation for the differences between known delta opioid agonists and antagonists.     

Comparison of the amino acid residues in the sixth transmembrane domains accessible in the binding-site crevices of mu, delta, and kappa opioid receptors.

We have mapped the residues in the sixth transmembrane domains (TMs 6) of the mu, delta, and kappa opioid receptors that are accessible in the binding-site crevices by the substituted cysteine accessibility method (SCAM). We previously showed that ligand binding to the C7.38S mutants of the mu and kappa receptors and the wild-type delta receptor was relatively insensitive to methanethiosulfonate ethylammonium (MTSEA), a positively charged sulfhydryl-specific reagent. These MTSEA-insensitive constructs were used as the templates, and 22 consecutive residues in TM6 (excluding C6.47) of each receptor were mutated to cysteine, 1 at a time. Most mutants retained binding affinities for [3H]diprenorphine, a nonselective opioid antagonist, similar to that of the template receptors. Treatment with MTSEA significantly inhibited [3H]diprenorphine binding to 11 of 22 mutants of the rat mu receptor and 9 of 22 mutants of the human delta receptor and 10 of 22 mutants of the human kappa receptor. Naloxone or diprenorphine protected all sensitive mutants, except the A6.42(287)C mu mutant. Thus, V6.40, F6.44, W6.48, I6.51, Y6.54, V6.55, I6.56, I6.57, K6.58, and A6.59 of the mu receptor; F6.44, I6.51, F6.54, V6.55, I6.56, V6.57, W6.58, T6.59, and L6.60 of the delta receptor; and F6.44, W6.48, I6.51, F6.54, I6.55, L6.56, V6.57, E6.58, A6.59, and L6.60 of the kappa receptor are on the water-accessible surface of the binding-site crevices. The accessibility patterns of residues in the TMs 6 of the mu, delta, and kappa opioid receptors are consistent with the notion that the TMs 6 are in alpha-helical conformations with a narrow strip of accessibility on the intracellular side of 6.54 and a wider area of accessibility on the extracellular side of 6.54, likely due to a proline kink at 6.50 that bends the helix in toward the binding pocket and enables considerable motion in this region. The wider exposure of residues 6.55-6.60 to the binding-site crevice, combined with the divergent amino acid sequences, is consistent with the inferred role of residues in this region in determining ligand binding selectivity. The conservation of the accessibility pattern on the cytoplasmic side of 6.54 suggests that this region may be important for receptor activation. This accessibility pattern is similar to that of the D2 dopamine receptor, the only other GPCR in which TM6 has been mapped by SCAM. That opioid receptors and the remotely related D2 dopamine receptor have similar accessibility patterns in TM6 suggest that these segments of GPCRs in the rhodopsin-like subfamily not only share secondary structure but also are packed similarly into the transmembrane bundle and thus have similar tertiary structure.     

Molecular modeling of the human serotonin(1A) receptor: role of membrane cholesterol in ligand binding of the receptor

Serotonin(1A) receptors are important neurotransmitter receptors and belong to the superfamily of G-protein coupled receptors (GPCRs). Although it is an important drug target, the crystal structure of the serotonin(1A) receptor has not been solved yet. Earlier homology models of the serotonin(1A) receptor were generated using rhodopsin as a template. We have used two recent crystal structures of the human β(2)-adrenergic receptor, one of which shows specific cholesterol binding site(s), as templates to model the human serotonin(1A) receptor. Since the sequence similarity between the serotonin(1A) receptor and β(2)-adrenergic receptor is considerably higher than the similarity between the serotonin(1A) receptor and rhodopsin, our model is more reliable. Based on these templates, we generated models of the serotonin(1A) receptor in the absence and presence of cholesterol. The receptor model appears more compact in the presence of cholesterol. We validated the stability of 'compactness' using coarse-grain MD simulation. Importantly, all ligands exhibit higher binding energies when docked to the receptor in the presence of cholesterol, thereby implying that membrane cholesterol facilitates ligand binding to the serotonin(1A) receptor. To the best of our knowledge, this is one of the first reports in which lipid-specific receptor conformations have been modeled by homology modeling.