Next Article in Journal
Hyperbaric Oxygen Therapy Reduces Oxidative Stress and Inflammation, and Increases Growth Factors Favouring the Healing Process of Diabetic Wounds
Previous Article in Journal
Analysis of Ku70 S155 Phospho-Specific BioID2 Interactome Identifies Ku Association with TRIP12 in Response to DNA Damage
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Three-Dimensional Structural Insights Have Revealed the Distinct Binding Interactions of Agonists, Partial Agonists, and Antagonists with the µ Opioid Receptor

1
National Center for Toxicological Research, US Food and Drug Administration, Jefferson, AR 72079, USA
2
Center for Drug Evaluation and Research, US Food and Drug Administration, Silver Spring, MD 20903, USA
3
National Center for Advancing Translational Sciences, National Institutes of Health, Bethesda, MD 20892, USA
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(8), 7042; https://doi.org/10.3390/ijms24087042
Submission received: 14 March 2023 / Revised: 9 April 2023 / Accepted: 9 April 2023 / Published: 11 April 2023
(This article belongs to the Section Molecular Biology)

Abstract

:
The United States is experiencing the most profound and devastating opioid crisis in history, with the number of deaths involving opioids, including prescription and illegal opioids, continuing to climb over the past two decades. This severe public health issue is difficult to combat as opioids remain a crucial treatment for pain, and at the same time, they are also highly addictive. Opioids act on the opioid receptor, which in turn activates its downstream signaling pathway that eventually leads to an analgesic effect. Among the four types of opioid receptors, the µ subtype is primarily responsible for the analgesic cascade. This review describes available 3D structures of the µ opioid receptor in the protein data bank and provides structural insights for the binding of agonists and antagonists to the receptor. Comparative analysis on the atomic details of the binding site in these structures was conducted and distinct binding interactions for agonists, partial agonists, and antagonists were observed. The findings in this article deepen our understanding of the ligand binding activity and shed some light on the development of novel opioid analgesics which may improve the risk benefit balance of existing opioids.

1. Introduction

The opioid crisis is one of the most prominent and severe public health issues in U.S. history owing to its rapidly changing nature. It has affected all types of communities across the country and indeed factors including diverse geography and demography, high multiplicity and complexity of its causes, and devastating consequences from opioid misuse and opioid use disorder (OUD) have all contributed to the difficulty of combating this crisis. According to the Centers for Disease Control and Prevention (CDC) Wonder Multiple Cause of Death database (https://wonder.cdc.gov/, accessed on 13 February 2023), the number of deaths caused by opioids increased annually from 1999 to 2020. In 1999, there were 6984 opioid-related deaths, while in 2020, the number of fatalities significantly increased by approximately 14-fold to 94,371. Opioid-related deaths have been separated by the CDC Wonder database (https://wonder.cdc.gov/wonder/help/mcd.html#UCD%20-%20ICD-10%20Codes, accessed on 13 February 2023) into four specific categories: natural or semisynthetic opioids, synthetic opioids, heroin, and mental and behavioral disorders. Natural opioids such as morphine and codeine and semisynthetic opioids such as oxycodone and hydrocodone are frequently prescribed for pain reduction. Synthetic opioids such as fentanyl are often prescribed for treating severe pain, typically advanced cancer pain, although the most recent cases of harm are linked to illicitly manufactured synthetic products. Since 2013, synthetic opioids have contributed the most to the rapidly increasing trend of all opioid-related deaths and are responsible for 42% of all mortalities. The opioid crisis has not only resulted in increasing fatalities throughout the years but also financially impacts society. In November 2017, the White House suggested that previous estimations of the economic impact of the opioid crisis had been extremely underestimated in the Council of Economic Advisers’ report [1]. The annual cost of the opioid crisis increased by approximately six-fold when all losses were considered [1]. Florence et al. estimated $1.02 trillion as the economic cost of the opioid crisis, in which the number of lives lost due to opioid overdose ($480.8 billion) and the reduced quality of life due to OUD ($390.0 billion) have accounted for over 85% of the total economic burden [2]. In addition, the healthcare spending cost accounted for almost $35 billion and moreover, the criminal justice-related spending accounted for $14.8 billion [2].
Such high costs, as well as the dangers of opioid addiction, have raised public awareness of opioid misuse. Many efforts and interventions have been devoted to targeting the opioid crisis including the development of opioid prescribing guidelines; government-funded programs for OUD prevention, treatment, and recovery; and research into various aspects associated with the crisis. A recent opioid guideline was developed by CDC to explicitly address the concerns of both opioid overdose and chronic pain (https://www.cdc.gov/mmwr/volumes/71/rr/rr7103a1.htm, accessed on 14 February 2023). Government programs such as the Opioid Abuse and Overdose Prevention program of the CDC and the State Targeted Response to the Opioid Abuse Crisis program of the Substance Abuse and Mental Health Services Administration (SAMHSA) have received increased funding in response to the continuing opioid crisis. Opioid-related research has also grown, as evidenced by the surge in the number of article titles in PubMed containing the word “opioid” since 2017.
Most opioid-related research has focused on developing new drugs for pain management and OUD treatment. Unlike drugs that are only used illicitly, access to opioids cannot be overly restrictive as they remain some of the most potent painkillers in modern medicine. To reduce the potential for harm caused by opioids for pain management, alternative drugs that alleviate pain with reduced or no addiction risk are under development. Another strategy is to study the mechanism of opioid addiction and optimize current treatment regimens by combining it with other drugs and reducing doses of the opioid component. Both strategies have been investigated by researchers from a variety of perspectives and continue to be explored.
In order to elucidate the mechanisms of action of opioids, researchers have used pharmacological and biochemical approaches to investigate the signaling pathways of opioids. The experimental data were able to verify that opioids elicit their primary effects by binding to the opioid receptor, which in turn inhibits pain transmission upon activation [3,4]. Therefore, the opioid receptor plays a crucial role in the pathway involving pain transmission, making it an important drug target for pain therapy. Although opioids are highly potent painkillers, all opioids can not only cause addiction but can also manifest in dependence, tolerance, withdrawal effects, and adverse effects including respiratory depression, sedation, and constipation [4]. A major goal for researchers targeting opioid receptors is the development of new compounds that have high analgesic efficacy with fewer side effects and a lower risk of tolerance, dependence, and addiction.
The structure–activity relationship between the opioid receptor and its ligands has been extensively studied, and researchers have been able to develop derivatives of naturally occurring opioids and synthetic opioids. However, these compounds did not have much success in overcoming the adverse effects of existing opioids. With advancements in X-ray crystallography and cryogenic electron microscopy (cryo-EM), more protein structures have been resolved in recent years, including the opioid receptors. The availability of 3D protein structures enables new perspectives and allows new approaches in drug design and development, including, but not limited to, the discovery of cryptic pockets on the receptor as potential binding sites, allosteric binding sites leading to allosteric modulation of the receptor, protein–protein interaction interfaces emerging as drug targets, and polypharmacology design. Overall, 3D protein structures provide valuable insights into the development of more specific and effective treatments [5,6,7,8]. In this article, we review recently solved structures of the μ subtype of the opioid receptor. In particular, we analyze the binding site of the μ opioid receptor and ligand–receptor interactions, as well as the allosteric modulation of the μ opioid receptor.

2. Biology of Opioid Receptors

The hypothesis of opioid receptor existence dates back to the 1950s, and it was based on the rigid structural activity relationship of opioids [4]. The concept of selective recognition sites has led to extensive studies on the family of opioid receptors, including receptor expression, cloning, and ligand binding [9,10,11]. However, the discovery of the opioid receptor was not experimentally demonstrated until 1973 when three laboratories reported the stereoselectivity of the opioid receptor in binding assays utilizing different radioligands: [3H]naloxone [12], [3H]dihydromorphine [13], and [3H]etorphine [14]. Further investigation using pharmacologic experiments revealed the possibility of multiple subtypes of opioid receptor: μ for morphine type, δ for enkephalin type, κ for ketocyclazocine type, σ for SKF10047 (N-allylnormetazocine) type, and ε for β-endorphin type [15].
The classification of the opioid receptor subtypes has evolved as our understanding of the molecular aspects of opioid receptors has profoundly grown. After two decades of the biochemical demonstration of opioid receptors, four distinct complementary DNAs (cDNAs) were isolated and identified as members of the opioid receptor family [4]. Three of the cDNAs were later correlated to the pharmacologically defined μ-(MOR), δ-(DOR), and κ-opioid receptors (KOR), while the fourth receptor did not bind opioid ligands with high affinity [15]. Instead, a novel peptide, nociceptin/orphanin FQ was identified as the endogenous ligand for the fourth receptor in the opioid receptor family, hence the name nociceptin opioid receptor (NOR, originally named ORL-1) [16]. The amino acid sequence similarity across these four subtypes of opioid receptors is approximately 50% [15].
Opioid receptors are primarily localized in pain-modulating descending pathways, including the medulla locus coeruleus and the periaqueductal gray area [17]. They are also expressed in limbic, midbrain, and cortical structures [17,18]. Opioid receptor activation at these sites directly inhibits neurons, which leads to the inhibition of spinal cord pain transmission [17]. Meanwhile, both opioid receptors and opioid peptides can be found throughout the nociceptive neural circuitry and critical regions of the central nervous system that are involved in reward and emotion [17,18]. The endogenous opioid peptides bind to the opioid receptor with different specificity. Endomorphin-1 and -2 selectively bind to the MOR, and dynorphin-A and -B and α-neoendorphin only bind to KOR. In contrast, β-endorphin and enkephalins bind to both MOR and DOR (Table 1) [3,19,20].
In order to study the mechanism of action of the opioid receptor, elucidating its structure becomes crucial. Molecular cloning and hydropathy analysis predicted that the primary structure of the opioid receptor is composed of seven transmembrane segments, which is a structural characteristic of G protein-coupled receptors (GPCR) [15]. This prediction was validated in 2012 by X-ray crystallography [21]. The overall architecture of the opioid receptor resembles the structure of the GPCRs with seven transmembrane α-helices connected by three intracellular loops (ICL1-3) and three extracellular loops (ECL1-3) [21].
As members of the non-olfactory class A GPCRs, all four known subtypes of opioid receptors are coupled to pertussis toxin-sensitive G proteins (Gi or Go or both) [17,19,20,22]. Upon agonist binding to the orthosteric binding site, the conformational change of the receptor prompts the activation of intracellular G proteins. Subsequently, the Gαi and Gβγ subunits dissociate from each other and then act on various downstream intracellular signaling pathways [17,19,20,22]. For instance, translocation of the Gα subunit is followed by the inhibition of the adenylyl cyclase (AC) activity, thereby inhibiting the formation of cAMP [23,24,25]. cAMP activates protein kinase A (PKA), a protein that is responsible for the phosphorylation of various ion channels, proteins, and enzymes, which eventually leads to their activation or inhibition [26]. The release of the Gβγ subunit leads to the inhibition of the voltage-gated Ca2+ channels (VGCC, L-type, and N-type) in addition to K+ channel activation (Figure 1) [17,24,27].
GPCRs are known for their sophisticated signaling pathways and conformational landscapes. Ligands may preferentially activate or inhibit a particular pathway by altering the conformations of the receptor, and this is referred to as functional selectivity or ligand bias [28,29]. The functional selectivity of the opioid receptor has been studied extensively due to potential therapeutic benefits of biased agonism and partial agonism. The opioid receptor can have G protein dependent signaling, β-arrestin dependent signaling, and G protein β-arrestin complex dependent signaling [30,31,32]. The G protein dependent signaling leads to a reduction of cAMP, a decreased Ca2+ response, and the activation of G. protein-coupled inwardly rectifying potassium (GIRK) channels [33]. In contrast, the phosphorylation of opioid receptors by G protein-coupled receptor kinases (GRK) leads to the recruitment of β-arrestin [34]. The binding of β-arrestin, in turn, regulates receptor desensitization, sequestration, sorting, internalization, and degradation [34]. In this case, the β-arrestin dependent signaling terminates the G protein signaling as the opioid receptor undergoes internalization and degradation [35]. Phosphorylated β-arrestin opioid receptor complex can also recruit other signal transduction cascades such as the mitogen-activated protein kinase (MAPK) pathway and the p38 pathway [36,37]. Moreover, the formation of the Gαi subunit and the β-arrestin complex also mediate various signaling pathways, including the activation of extracellular signal-regulated kinase (ERK) (Figure 2) [37].

3. Structures of μ Opioid Receptor (MOR)

Several structures of MOR with atomic resolution have been reported since 2012 and are listed in Table 2.
These include two structures in complex with antagonists β-funaltrexamine (β-FNA) (PDB ID: 4DKL) [21] and alvimopan (PDB ID: 7UL4) [38]; five structures in complex with partial agonists PZM21 (PDB ID: 7SBF, 8EFO) [39,40], FH210 (PDB ID: 7SCG) [39], oliceridine (TRV130) (PDB ID: 8EFB) [40], and SR17018 (PDB ID: 8EFL) [40]; and six structures bound with full agonists BU72 (PDB ID: 5C1M) [41], [D-Ala2, N-MePhe4, Gly-ol5] enkephalin (DAMGO, a synthetic peptide) (PDB ID: 6DDE, 6DDF, and 8EFQ) [40,42], fentanyl (PDB ID: 8EF5) [40], morphine (PDB ID:8EF6) [40], lofentanil (PDB ID: 7T2H) [43], and mitragynine pseudoindoxyl (PDB ID: 7T2G) [43]. Another recently solved structure is bound with a bitopic ligand (PDB ID: 7U2L) [8]. The two structures reported earlier in 2012 [21] and 2015 [41] utilized an X-ray diffraction method and the stabilization of the receptor alone was difficult, therefore the T4 lysozyme (T4L) fusion protein strategy and nanobody were used to help stabilize the protein for crystallization. The remaining structures were from cryo-EM experiments and the receptors were solved in complex with the G proteins and/or megabodies (Table 2).
The MOR consists of seven transmembrane helices that are connected by loops ECL1-3 and ICL1-3. The orthosteric binding site is located near the extracellular side of the helix bundle. The alignment of the structures of MOR revealed that all ligands were bound within the same site (Figure 3). Intriguingly, this binding site is largely exposed to the extracellular surface, whereas in other GPCRs the ligands were often buried deeper within the helix bundle. This structural characteristic may explain the rapid dissociation kinetics of extremely potent opioids. For instance, the inhibition constant of alvimopan is 350 pM with a half-life of 30 min [21].
Both agonists and antagonists bind to the same orthosteric site of the MOR, and the structural differences presented in these structures were relatively subtle. In particular, the binding modes of the antagonist β-FNA and full agonist BU72 are very similar. Both ligands have a morphinan scaffold with a similar binding orientation (Figure 4). Two conserved interactions between the ligands and the receptor were observed: (1) a water-mediated hydrogen bond between the phenolic hydroxyl of the morphinan group and H2976.52 (superscripts are Ballesteros–Weinstein numbering in GPCRs, all residue numbering scheme in this manuscript refers to the mouse MOR, there are also human MOR cryo-EM structures available with different residue number scheme) and (2) an ionic interaction between the tertiary amine of the morphinan group and D1473.32. Meanwhile, the side chain of D1473.32 also forms a hydrogen bond with Y3267.43. Despite the structural differences in BU72 and DAMGO, the orientation of the residues that interacted with the full agonists are highly similar. The ionic interaction with D1473.32 was preserved as the N-terminus of DAMGO formed a salt bridge with the amine group of D1473.32, and the same amine group also formed a hydrogen bond with Y3267.43 in this structure (Figure 4). The same water-mediated hydrogen bond network between H2976.52 and the phenol group of DAMGO was observed in the molecular dynamics (MD) simulation [37]. Although β-FNA and alvimopan are both antagonists, unlike the morphinan scaffold of β-FNA, alvimopan possesses phenol-piperidine moieties which share similarities with the structure of fentanyl and its analogs. The protonated piperidine forms a salt bridge with the D1473.32 (Figure 4). The recently solved fentanyl-bound MOR structure has shown a direct ᴨ–ᴨ stacking interaction. between the benzene ring of fentanyl and the side chain of W2936.48 and Y3267.43, which is not presented in the morphine-bound MOR structure [40]. Additionally, a hydrophobic interaction between the phenylethyl moiety of fentanyl and a minor pocket located between TM2 and TM3 of the MOR is also absent in the morphine-bound MOR structure [40]. Mutations of the residues within this minor pocket including residues Q1242.60, W13323.50, and I1443.27 to alanine affected the potency of fentanyl more than that of morphine [40]. These atomic details revealed by cryo-EM structures may explain the higher potency of fentanyl to the receptor. The 2D structures of all ligands are shown in Table 3.
Although the binding pocket and binding interaction are highly similar between full agonist bound structures and antagonist bound structures, the partial agonist bound structures maintained the same interactions while also showing a distinctive feature. Both partial agonists PZM21 and FH210 retained the polar interactions between their basic amine and D1473.32 and water mediated hydrogen bond between their phenol hydroxy group and H2976.52 (Figure 4). Intriguingly, both partial agonists also showed high complementarity to a lipophilic vestibule formed by the extracellular surface of transmembrane helices (TM2 and TM3) and ECL1. The interactions between the thiolphenylalkyl moiety of PZM21 and naphthyl moiety of FH210 with this lipophilic vestibule may help to explain their partial agonism and the increase in G protein biased signaling over β-arrestin recruitment.
Despite the overall structural similarities shared between the structures of MOR, the distinctive binding feature of the partial agonists raised the question of whether there are other differences between the binding pockets. To answer this question, we first analyzed the composition of the surrounding residues within 5 Å of the ligands in all structures. The number of residues within 5 Å of the ligands is between 15 to 20, half of which are hydrophobic, and the remaining half are either charged residues or polar residues (Table 4). This analysis suggested that the residue composition surrounding the ligands is very similar. The structure in complex with the bitopic ligand (PDB ID: 7U2L) is not included in our analyses due to the undetermined pharmacological effect of the bitopic ligand. Next, instead of focusing on the type of residues surrounding the ligands, we categorized the residues by ligand types. In other words, we investigated if the residues are shared by all ligands, between two types of ligands, or unique to one type of ligand. We also narrowed down the search by limiting the distance between the ligand and the receptor to 4 Å. A total of 11 residues were found in all structures. W13323.50, W2936.48, H2976.52, K3036.58, and W3187.35 were shared between structures in complex with either full agonists or antagonists. V1433.28 and C21745.50 were shared between structures in complex with either full agonists or partial agonists. Y751.39, T21845.51, H3197.35, and G3257.41 were the residues unique to the binding of full agonists. T1202.56 and K2335.39 were the residues found within 4 Å of the partial agonist that are distinctive, and L2195.52 and E2295.35 were unique to the antagonists (Figure 5). We computed the contact area for the ligands, revealing that DAMGO has the largest contact area of 147 Å2 and BU72 has the smallest contact area of 40 Å2. As mentioned above, water molecules play an important role in mediating hydrogen bonding interactions between the receptor and the ligand. We also examined the water molecules in the structures. Out of the 15 structures, 10 structures did not solve any water molecules (PDB ID: 6DDE, 6DDF, 7SBF, 7T2G, 8EF5, 8EF6, 8EFB, 8EFL, 8EFQ, and 8EFO), while the other 5 structures (PDB ID: 4DKL, 7UL4, 5C1M, 7SCG, and 7T2H) have 2, 1, 7, 2, and 8 water molecules solved, respectively.
To further analyze the binding site, we measured the solvent-accessible surface area and solvent-accessible volume for each structure. Although residues near the ligand exhibit similar conformations, binding of the ligand may induce a conformational change in residues of the second interaction shell. Therefore, the binding of different ligands in the same binding site may still lead to subtle differences in the solvent-accessible surface area and volume of nearby residues. Table 5 lists the solvent-accessible surface area and volume of the binding site for all structures. CASTp 3.0 was used for SASA and SAV calculation [47]. There was no significant difference found in the solvent-accessible surface area among all structures, but the full agonist bound structures have the smallest solvent-accessible volume compared to that of the antagonist bound structures or partial agonist bound structures. This characteristic may be explained by the induced-fit mechanism. Although both agonists and antagonists bind to the same binding site, the size of the ligands varies, leading to differences in the shape complementarity of the binding pocket.
Molecular dynamics simulation is a technique often employed to study protein dynamics with atomic details, which provides a deeper understanding of small molecule binding and protein functional selectivity mechanisms [48,49,50,51]. Multiple studies have endeavored to investigate ligand binding and selectivity of the MOR, the role of water molecules in the MOR, the activation mechanism of the MOR, and biased signaling of the MOR using MD simulations [52,53,54,55,56,57,58,59,60,61,62,63]. Here, we summarized the findings focused on the interactions involved in ligand binding and how water molecules were involved in ligand binding to the MOR. Both Liao et al. and Podlewska et al. discovered that the partial agonist PZM21 interacts strongly with D1473.32, Y1483.33, and Y3267.43, which is congruent with the interactions found in the cryo-EM structures [64,65]. Intriguingly, in the MD simulation trajectories Podlewska et al. observed that the increased intensity of the interactions between PZM21 with I2966.51 and H2976.52 prompted the increased intensity of interaction with W3187.34 [64]. Valerylfentanyl, a fentanyl analog often implicated in opioid overdose, is a partial agonist of the MOR [66]. MD simulations have revealed that the alkyl chain of this compound is not well accommodated by the active state of the MOR, which may shift the receptor toward an inactive state [66]. Cheng et al. observed that oliceridine (TRV130) formed a hydrogen bond with D1473.32 and a direct interaction with Y3267.43, which maintained the close distance between Y3267.43 and W2936.48, thus enhancing the stability of W2936.48 [56]. Shim et al. showed that the salt bridge between β-FNA and D1473.32 was consistently maintained in MD simulations [52]. Overall, the interactions between the ligands and the MOR found in MD simulations and cryo-EM structures were consistent, while the MD simulations provided additional atomic details that may help to explain the potential mechanism of action of the partial agonism.
Fourier-transform infrared spectroscopy and UV-visible spectroscopy data have shown that internal water molecules are essential for the function of various membrane proteins including GPCRs [67,68]. It has been proposed that internal water molecules play an important role in GPCR activation and signaling, but the mechanism behind these findings remains elusive [69]. Yuan et al. performed MD simulations on the MOR to study the role of water molecules in MOR activation. Two systems were built: an agonist- (morphine) bound system and an antagonist- (levallorphan) bound system. They discovered that the agonist-bound system has more intrinsic water molecules compared to the antagonist-bound system [53]. Closer examination revealed that besides the hydrogen bond with D1473.32, levallorphan also formed a σ–ᴨ stacking interaction with Y3207.43 that prohibits water penetration [53]. These findings suggest that the size of the ligand, as well as the degree of “openness” of the binding pocket, may play a role in determining whether a ligand acts as an agonist or antagonist. Consistent with this idea, the agonist-bound structures have a smaller solvent accessible volume due to a smaller ligand size and greater exposure of the binding pocket, whereas the antagonists are larger in size, requiring the binding pocket to expand, resulting in a larger solvent accessible volume. Later, the same group published a second paper with five different systems: the apo system, agonist (morphine) bound system with and without sodium, and antagonist (naltrexone) bound system with and without sodium [55]. The analysis showed that water molecules within the binding site frequently exchange with the extracellular water molecules in the agonist-bound systems, which was not observed in the antagonist-bound systems [55]. The water molecule movement near the conserved NPxxY motifs has been associated with the switches of GPCR activation, and its importance has been demonstrated in the MD simulations [53,70,71,72]. The importance of the internal water environment near the NPxxY region has also been demonstrated by Cheng et al., in that water molecules in the agonist-bound system are much more active compared to that of the antagonist-bound system [56]. All these findings are consistent with the crystal structure data where the structure of the MOR in complex with agonist BU72 (PDB ID: 5C1M) [41] was solved with water molecules around the NPxxY region, proving the importance of water in MOR activation.

4. Allosteric Modulation of the MOR

With the ongoing opioid crisis, the urgency for understanding the allosteric modulation of opioid receptors to guide the development of effective and safe drugs for pain management has increased. Allosteric modulators are ligands that bind to the non-orthosteric binding site of the receptor, either enhancing (positive allosteric modulator or PAM) or diminishing (negative allosteric modulator or NAM) the signaling activities of the receptor. Allosteric modulators can alter the affinity, potency, and efficacy of the orthosteric ligands. PAMs to MOR are considered as new candidates for a safer approach to pain management. Ideally, PAMs would indirectly activate the MOR by enhancing the effects of endogenous opioid peptides, leading to the desired analgesic effect with a lower risk of severe side effects or addiction liability. Meanwhile, PAMs might increase the efficacy or potency of opioid drugs such as morphine, so the same therapeutic effect can be achieved at a lower dose with fewer side effects. It is worth noting that allosteric modulators can sometimes act as biased ligands. Bias refers to the phenomenon in which a ligand selectively activates certain signaling pathways downstream of the receptor. This phenomenon can occur for both orthosteric and allosteric ligands. In this case, the binding of PAMs may alter the preferred signaling pathway of the receptor, leading to biased downstream signaling. The advantage is the avoidance of the pathway that can potentially lead to undesired effects. For example, the β-arrestin-biased signaling of the MOR has been connected to respiratory depression and constipation [73,74,75]. Additional studies found controversy surrounding the connection between β-arrestin signaling and respiratory depression [76]. Therefore, the connection between β-arrestin-biased signaling and respiratory depression remains an area of ongoing investigation, and further studies are needed to fully understand the role of this pathway in the development of these side effects. However, PAMs may still have the potential advantage of steering the downstream signaling toward G protein-biased signaling without activating the β-arrestin-biased signaling pathway. Additionally, allosteric modulators can sometimes activate the receptor directly. Therefore, the allosteric activity is not only dependent on the binding affinity of the modulator itself and the allosteric cooperativity, which characterizes the capacity of the modulator to alter the affinity and/or efficacy of the orthosteric ligand, but also dependent on the orthosteric ligand as well [77]. The phenomenon that the same allosteric modulator has different effects depending on the orthosteric ligand is referred to as probe dependence [77]. On the other hand, two NAMs were discovered prior to the development of MOR PAMs. One of these NAMs is cannabidiol, which is a CB1 receptor agonist and was found to be a NAM of both MOR and DOR agonists [78]. Cannabidiol was shown to accelerate the dissociation of DAMGO from the MOR by a factor of at least 12 and also accelerated the dissociation of [3H]naltrindole from the DOR by a factor of at least 2 [78]. In a high-throughput screen using a β-arrestin recruitment assay, salvinorin A, a potent KOR agonist, was also identified as a NAM of MOR [79]. The potential use of NAMs may be as a treatment for opioid use disorder (OUD). As mentioned earlier, NAMs bind to the allosteric binding site of the receptor, weakening or reducing signaling activity of the receptor by slowing or inhibiting the binding of the ligands to the orthosteric binding site. In the case of the MOR, NAMs can reduce the signaling activity of the receptor, leading to a reduction in the effects of opioid agonists. This may have potential therapeutic benefits for OUD treatment. Currently, there are no MOR NAMs approved for the treatment of OUD, and more research is needed to fully explore their potential therapeutic benefits and safety profile.
Monovalent and divalent cations, in particular the monovalent sodium cation (Na+) and the divalent magnesium cation (Mg2+), are well-known PAMs and NAMs for the MOR [61]. Interestingly, Na+ and Mg2+ have opposite effects, with Na+ acting as a NAM by reducing agonist binding and Mg2+ acting as a PAM by affecting MOR signaling [61]. The idea of Na+ as a NAM has been established for decades, and the hypothesis was that Na+ stabilizes the inactive conformation of the receptor, which later was extended to several GPCRs including the MOR [53,58]. Both high-resolution X-ray structures of GPCRs and multiple MD simulations have supported this hypothesis [58,80,81,82,83,84,85,86]. The recently solved structure of inactive DOR revealed that the allosteric binding site of Na+ is coordinated by two water molecules and residues N1313.35, S1353.39, and D952.50 (PDB ID: 4N6H [80]) (Figure 6). Binding events of Na+ to this allosteric binding site were also observed in MD simulations in both inactive MOR and KOR systems [84]. This allosteric binding site was found to be collapsed in the active GPCR structures [81]. A recent study of a prototypic GPCR and adenosine A2A receptor revealed the unexpected results that Na+ also stabilizes the intermediate state that may correlate to partial agonism, which may explain why Na+ was also found to promote agonist-induced MOR activation [87]. The same study also revealed that Mg2+ promotes the opening of the G protein binding cleft, leading to the bridging between the acidic residues located on the extracellular region of TM5 and TM6, and resulting in receptor activation [87]. MD simulation studies have discovered that the preferred binding site of Mg2+ is located near the extracellular region of the MOR, in particular in proximity to residues S214ECL2, D216ECL2, and E310ECL3 (Figure 6) [61]. Additionally, Mg2+ exhibits a higher binding affinity to the active conformation of the MOR compared to the inactive conformation [61].
The conformational state of the MOR is altered not only by cations but also by other endogenous molecules such as lipids. As a membrane-embedded protein, the translocation and internalization of MOR are controlled by lipids while its function is also modulated by lipids [88], especially cholesterol, a sterol-type of lipid that comprises approximately 30% of the cell membrane [89]. Experimental data have shown that the removal of cholesterol reduces cAMP signaling of the MOR [90]. One of the ways that cholesterol can modulate the activity of the MOR is by binding directly to the receptor which leads to allosteric effects; the other is by changing the fluidity of the membrane, which in turn affects the signaling of the receptor [89]. Cholesterol was found co-crystallized with many GPCRs, including the MOR. In both active (PDB ID: 5C1M [41]) and inactive (PDB ID: 4DKL [21]) MOR structures, cholesterol was found to be co-crystallized at the same binding site. This binding site was also found in the active structures of KOR (PDB ID: 6PT2 [91], 6B73 [92]), the structure of the adenosine A2A receptor (PDB ID: 5IU4 [93]), and the active structure of the CXCR3 receptor (PDB ID: 5WB2 [94]). This binding site is a groove between TM2, TM3, and TM4 or TM6 at the extracellular region of the receptor (Figure 6).
Endogenous molecules are not the only allosteric modulators of MOR; many small molecules were either serendipitously found or intentionally designed to act as PAMs for the MOR. The first selective small molecule PAMs of the MOR are BMS-986121 and BMS-986122, both of which were identified by high-throughput screening in 2013 [95]. Both molecules were first screened for the ability to enhance the binding of endogenous agonist endomorphin-1 that recruits β-arrestin to MOR [95]. However, further studies found these compounds enhanced the affinity and/or efficacy of various opioid agonists. Consequently, more analogs of these two compounds were developed along with SAM (silent allosteric modulator) BMS-986124 [95]. PAM BMS-986122 exhibits probe dependence: it increases the potency and affinity of methadone and DAMGO but only enhanced the efficacy of morphine and nalbuphine, with no effect on the binding of antagonists [77,96]. To date, the published BMS series of compounds only act as PAMs or SAMs, and the structure–activity relationships remain unclear. However, there have been several attempts to identify the binding site of the BMS compoundseries by docking and MD simulations. Bartuzi et al. revealed that the potential binding site of BMS-986122 may be slightly above the orthosteric binding site and towards the extracellular region of TM2 and TM7 (Figure 6) [97]. Recently, Bartuzi et al. developed a novel small molecule PAM of the MOR that can enhance the efficacy of morphine [98]. They also predicted this novel compound binds in proximity to TM2 and the binding pocket is composed of W13523.50, C1423.25, V1453.28, I1463.29, and C21945.50 (Figure 6) [98].
Figure 6. The in silico prediction of potential allosteric binding sites for PAMs and NAMs including cholesterol, small molecules, and ions. Each color depicts a potential site, and five sites are illustrated [61,80,94,97,98].
Figure 6. The in silico prediction of potential allosteric binding sites for PAMs and NAMs including cholesterol, small molecules, and ions. Each color depicts a potential site, and five sites are illustrated [61,80,94,97,98].
Ijms 24 07042 g006

5. Future Perspectives

Although medicinal chemists spent almost two centuries developing drugs targeting the MOR, many of the compounds still share the same scaffold as morphine or have similar synthetic scaffolds, which in turn limits the diversity of the active compounds. Moreover, this scaffold was reported to induce undesirable adverse effects, including respiratory depression and constipation. More recent approaches to developing new drugs targeting the MOR involve high-throughput screening and defined signaling assays in stable MOR transfected cell lines. However, these approaches are time-consuming and expensive. The advancement of X-ray crystallography and cryo-EM helped to solve more opioid receptor structures in recent years, which provided an unprecedented opportunity in computer-aided drug design. High-resolution structures of the MOR in complex with agonists or antagonists can confirm the binding site of the ligands and provide insights into the chemical recognition of the ligands. The molecular details of the binding site itself provide valuable information for structure-based drug design. Furthermore, some of the MOR structures have been crystallized with G proteins, which not only provides hints regarding the conformational changes that drive the downstream signaling pathway, but also in terms of the atomic details of the protein–protein interaction interface, which can potentially serve as a binding site for protein–protein interaction disruptors or molecular glues that stabilize the interface. Finally, these high-resolution structures have enabled researchers to explore novel allosteric binding sites. Recent improvements in computational chemistry have increased the accuracy of docking calculations, and with rapidly increasing computational power, researchers can explore the infinite chemical space and screen hundreds of millions of compounds. Structurally divergent compounds are likely to have unique biological properties that may potentially act as safer analgesics. Overall, the availability of high-resolution structures and a better understanding of the MOR structures may hold the key to successful opioid analogue discovery.

Author Contributions

Data curation and analysis, Z.L., J.L. and F.D.; writing—original draft preparation, Z.L.; writing—review and editing, all authors; supervision, T.A.P. and H.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This work was supported in part by an appointment to the Research Participation Program at the National Center for Toxicological Research (Zoe Li, Fan Dong), administered by the Oak Ridge Institute for Science and Education through an interagency agreement between the US Department of Energy and the US Food and Drug Administration. This work was also supported in part by the Intramural Research Program of the National Center for Advancing Translational Sciences, National Institutes of Health.

Conflicts of Interest

The authors declare no conflict of interest. This article reflects the views of the authors and does not necessarily reflect those of the National Institutes of Health or the U.S. Food and Drug Administration.

References

  1. Ropero-Miller, J.D.; Speaker, P.J. The Hidden Costs of the Opioid Crisis and the Implications for Financial Management in the Public Sector. Forensic Sci. Int. Synerg. 2019, 1, 227–238. [Google Scholar] [CrossRef] [PubMed]
  2. Florence, C.; Luo, F.; Rice, K. The Economic Burden of Opioid Use Disorder and Fatal Opioid Overdose in the United States, 2017. Drug Alcohol. Depend. 2021, 218, 108350. [Google Scholar] [CrossRef] [PubMed]
  3. Darcq, E.; Kieffer, B.L. Opioid Receptors: Drivers to Addiction? Nat. Rev. Neurosci. 2018, 19, 499–514. [Google Scholar] [CrossRef] [PubMed]
  4. Pasternak, G.W.; Pan, Y.-X. Mu Opioids and Their Receptors: Evolution of a Concept. Pharm. Rev. 2013, 65, 1257–1317. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Manglik, A.; Lin, H.; Aryal, D.K.; McCorvy, J.D.; Dengler, D.; Corder, G.; Levit, A.; Kling, R.C.; Bernat, V.; Hübner, H.; et al. Structure-Based Discovery of Opioid Analgesics with Reduced Side Effects. Nature 2016, 537, 185–190. [Google Scholar] [CrossRef] [Green Version]
  6. Zheng, Z.; Huang, X.-P.; Mangano, T.J.; Zou, R.; Chen, X.; Zaidi, S.A.; Roth, B.L.; Stevens, R.C.; Katritch, V. Structure-Based Discovery of New Antagonist and Biased Agonist Chemotypes for the Kappa Opioid Receptor. J. Med. Chem. 2017, 60, 3070–3081. [Google Scholar] [CrossRef] [Green Version]
  7. Poli, G.; Dimmito, M.P.; Mollica, A.; Zengin, G.; Benyhe, S.; Zador, F.; Stefanucci, A. Discovery of Novel Μ-Opioid Receptor Inverse Agonist from a Combinatorial Library of Tetrapeptides through Structure-Based Virtual Screening. Molecules 2019, 24, 3872. [Google Scholar] [CrossRef] [Green Version]
  8. Faouzi, A.; Wang, H.; Zaidi, S.A.; DiBerto, J.F.; Che, T.; Qu, Q.; Robertson, M.J.; Madasu, M.K.; El Daibani, A.; Varga, B.R.; et al. Structure-Based Design of Bitopic Ligands for the µ-Opioid Receptor. Nature 2023, 613, 767–774. [Google Scholar] [CrossRef]
  9. Portoghese, P.S. Relationships Between Stereostructure and Pharmacological Activities. Annu. Rev. Pharmacol. 1970, 10, 51–76. [Google Scholar] [CrossRef]
  10. Portoghese, P.S. A New Concept on the Mode of Interaction of Narcotic Analgesics with Receptors. J. Med. Chem. 1965, 8, 609–616. [Google Scholar] [CrossRef]
  11. Beckett, A.H.; Casy, A.F. 5 Analgesics and Their Antagonists: Biochemical Aspects and Structure-Activity Relationships. In Progress in Medicinal Chemistry; Elsevier: Amsterdam, The Netherlands, 1965; Volume 4, pp. 171–218. ISBN 978-0-444-53323-4. [Google Scholar]
  12. Pert, C.B.; Snyder, S.H. Opiate Receptor: Demonstration in Nervous Tissue. Science 1973, 179, 1011–1014. [Google Scholar] [CrossRef] [PubMed]
  13. Terenius, L. Stereospecific Interaction Between Narcotic Analgesics and a Synaptic Plasma Membrane Fraction of Rat Cerebral Cortex. Acta Pharmacol. Toxicol. 2009, 32, 317–320. [Google Scholar] [CrossRef] [PubMed]
  14. Simon, E.J.; Hiller, J.M.; Edelman, I. Stereospecific Binding of the Potent Narcotic Analgesic [ 3 H]Etorphine to Rat-Brain Homogenate. Proc. Natl. Acad. Sci. USA 1973, 70, 1947–1949. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Gropper, M.A.; Miller, R.D.; Cohen, N.H. (Eds.) Miller’s Anesthesia, 9th ed.; Elsevier: Philadelphia, PA, USA, 2020; ISBN 978-0-323-59604-6. [Google Scholar]
  16. Mollereau, C.; Parmentier, M.; Mailleux, P.; Butour, J.-L.; Moisand, C.; Chalon, P.; Caput, D.; Vassart, G.; Meunier, J.-C. ORL1, a Novel Member of the Opioid Receptor Family: Cloning, Functional Expression and Localization. FEBS Lett. 1994, 341, 33–38. [Google Scholar] [CrossRef] [PubMed]
  17. Al-Hasani, R.; Bruchas, M.R. Molecular Mechanisms of Opioid Receptor-Dependent Signaling and Behavior. Anesthesiology 2011, 115, 1363–1381. [Google Scholar] [CrossRef] [Green Version]
  18. Le Merrer, J.; Becker, J.A.J.; Befort, K.; Kieffer, B.L. Reward Processing by the Opioid System in the Brain. Physiol. Rev. 2009, 89, 1379–1412. [Google Scholar] [CrossRef]
  19. Pathan, H.; Williams, J. Basic Opioid Pharmacology: An Update. Br. J. Pain 2012, 6, 11–16. [Google Scholar] [CrossRef] [Green Version]
  20. Snyder, S.H.; Pasternak, G.W. Historical Review: Opioid Receptors. Trends Pharmacol. Sci. 2003, 24, 198–205. [Google Scholar] [CrossRef]
  21. Manglik, A.; Kruse, A.C.; Kobilka, T.S.; Thian, F.S.; Mathiesen, J.M.; Sunahara, R.K.; Pardo, L.; Weis, W.I.; Kobilka, B.K.; Granier, S. Crystal Structure of the Μ-Opioid Receptor Bound to a Morphinan Antagonist. Nature 2012, 485, 321–326. [Google Scholar] [CrossRef] [Green Version]
  22. Feng, Y.; He, X.; Yang, Y.; Chao, D.; H Lazarus, L.; Xia, Y. Current Research on Opioid Receptor Function. CDT 2012, 13, 230–246. [Google Scholar] [CrossRef] [Green Version]
  23. Wiesenfeld-Hallin, Z.; de Araúja Lucas, G.; Alster, P.; Xu, X.-J.; Hökfelt, T. Cholecystokinin/Opioid Interactions. Brain Res. 1999, 848, 78–89. [Google Scholar] [CrossRef] [PubMed]
  24. Pasternak, G.W. Multiple Opiate Receptors: Déjà vu All over Again. Neuropharmacology 2004, 47, 312–323. [Google Scholar] [CrossRef] [PubMed]
  25. Waldhoer, M.; Bartlett, S.E.; Whistler, J.L. Opioid Receptors. Annu. Rev. Biochem. 2004, 73, 953–990. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Sassone-Corsi, P. The Cyclic AMP Pathway. Cold Spring Harb. Perspect. Biol. 2012, 4, a011148. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, L.; Zhang, J.-T.; Hang, L.; Liu, T. Mu Opioid Receptor Heterodimers Emerge as Novel Therapeutic Targets: Recent Progress and Future Perspective. Front. Pharmacol. 2020, 11, 1078. [Google Scholar] [CrossRef]
  28. Tan, L.; Yan, W.; McCorvy, J.D.; Cheng, J. Biased Ligands of G Protein-Coupled Receptors (GPCRs): Structure–Functional Selectivity Relationships (SFSRs) and Therapeutic Potential. J. Med. Chem. 2018, 61, 9841–9878. [Google Scholar] [CrossRef]
  29. Yang, D.; Zhou, Q.; Labroska, V.; Qin, S.; Darbalaei, S.; Wu, Y.; Yuliantie, E.; Xie, L.; Tao, H.; Cheng, J.; et al. G Protein-Coupled Receptors: Structure- and Function-Based Drug Discovery. Signal Transduct. Target. Ther. 2021, 6, 7. [Google Scholar] [CrossRef]
  30. Cong, X.; Maurel, D.; Déméné, H.; Vasiliauskaité-Brooks, I.; Hagelberger, J.; Peysson, F.; Saint-Paul, J.; Golebiowski, J.; Granier, S.; Sounier, R. Molecular Insights into the Biased Signaling Mechanism of the μ-Opioid Receptor. Mol. Cell 2021, 81, 4165–4175.e6. [Google Scholar] [CrossRef]
  31. Che, T.; Dwivedi-Agnihotri, H.; Shukla, A.K.; Roth, B.L. Biased Ligands at Opioid Receptors: Current Status and Future Directions. Sci. Signal. 2021, 14, eaav0320. [Google Scholar] [CrossRef]
  32. Kelly, B.; Hollingsworth, S.A.; Blakemore, D.C.; Owen, R.M.; Storer, R.I.; Swain, N.A.; Aydin, D.; Torella, R.; Warmus, J.S.; Dror, R.O. Delineating the Ligand–Receptor Interactions That Lead to Biased Signaling at the μ-Opioid Receptor. J. Chem. Inf. Model. 2021, 61, 3696–3707. [Google Scholar] [CrossRef]
  33. Nagi, K.; Pineyro, G. Kir3 Channel Signaling Complexes: Focus on Opioid Receptor Signaling. Front. Cell Neurosci. 2014, 8, 186. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Mafi, A.; Kim, S.-K.; Goddard, W.A. Mechanism of β-Arrestin Recruitment by the μ-Opioid G Protein-Coupled Receptor. Proc. Natl. Acad. Sci. USA 2020, 117, 16346–16355. [Google Scholar] [CrossRef] [PubMed]
  35. Noor, N.; Patel, C.B.; Rockman, H.A. β-Arrestin: A Signaling Molecule and Potential Therapeutic Target for Heart Failure. J. Mol. Cell. Cardiol. 2011, 51, 534–541. [Google Scholar] [CrossRef] [Green Version]
  36. Macé, G.; Miaczynska, M.; Zerial, M.; Nebreda, A.R. Phosphorylation of EEA1 by P38 MAP Kinase Regulates μ Opioid Receptor Endocytosis. EMBO J. 2005, 24, 3235–3246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Belcheva, M.M.; Clark, A.L.; Haas, P.D.; Serna, J.S.; Hahn, J.W.; Kiss, A.; Coscia, C.J. μ and κ Opioid Receptors Activate ERK/MAPK via Different Protein Kinase C Isoforms and Secondary Messengers in Astrocytes. J. Biol. Chem. 2005, 280, 27662–27669. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Robertson, M.J.; Papasergi-Scott, M.M.; He, F.; Seven, A.B.; Meyerowitz, J.G.; Panova, O.; Peroto, M.C.; Che, T.; Skiniotis, G. Structure Determination of Inactive-State GPCRs with a Universal Nanobody. Nat. Struct. Mol. Biol. 2022, 29, 1188–1195. [Google Scholar] [CrossRef]
  39. Wang, H.; Hetzer, F.; Huang, W.; Qu, Q.; Meyerowitz, J.; Kaindl, J.; Hübner, H.; Skiniotis, G.; Kobilka, B.K.; Gmeiner, P. Structure-Based Evolution of G Protein-Biased Μ-Opioid Receptor Agonists. Angew. Chem. Int. Ed. 2022, 61, e202200269. [Google Scholar] [CrossRef]
  40. Zhuang, Y.; Wang, Y.; He, B.; He, X.; Zhou, X.E.; Guo, S.; Rao, Q.; Yang, J.; Liu, J.; Zhou, Q.; et al. Molecular Recognition of Morphine and Fentanyl by the Human μ-Opioid Receptor. Cell 2022, 185, 4361–4375.e19. [Google Scholar] [CrossRef]
  41. Huang, W.; Manglik, A.; Venkatakrishnan, A.J.; Laeremans, T.; Feinberg, E.N.; Sanborn, A.L.; Kato, H.E.; Livingston, K.E.; Thorsen, T.S.; Kling, R.C.; et al. Structural Insights into Μ-Opioid Receptor Activation. Nature 2015, 524, 315–321. [Google Scholar] [CrossRef] [Green Version]
  42. Koehl, A.; Hu, H.; Maeda, S.; Zhang, Y.; Qu, Q.; Paggi, J.M.; Latorraca, N.R.; Hilger, D.; Dawson, R.; Matile, H.; et al. Structure of the Μ-Opioid Receptor–Gi Protein Complex. Nature 2018, 558, 547–552. [Google Scholar] [CrossRef]
  43. Qu, Q.; Huang, W.; Aydin, D.; Paggi, J.M.; Seven, A.B.; Wang, H.; Chakraborty, S.; Che, T.; DiBerto, J.F.; Robertson, M.J.; et al. Insights into Distinct Signaling Profiles of the ΜOR Activated by Diverse Agonists. Nat. Chem. Biol. 2022, 19, 423–430. [Google Scholar] [CrossRef]
  44. Goddard, T.D.; Huang, C.C.; Meng, E.C.; Pettersen, E.F.; Couch, G.S.; Morris, J.H.; Ferrin, T.E. UCSF ChimeraX: Meeting Modern Challenges in Visualization and Analysis: UCSF ChimeraX Visualization System. Protein Sci. 2018, 27, 14–25. [Google Scholar] [CrossRef] [Green Version]
  45. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Meng, E.C.; Couch, G.S.; Croll, T.I.; Morris, J.H.; Ferrin, T.E. UCSF ChimeraX: Structure Visualization for Researchers, Educators, and Developers. Protein Sci. 2021, 30, 70–82. [Google Scholar] [CrossRef] [PubMed]
  46. Maestro; Maestro, Schrödinger, LLC: New York, NY, USA, 2021.
  47. Tian, W.; Chen, C.; Lei, X.; Zhao, J.; Liang, J. CASTp 3.0: Computed Atlas of Surface Topography of Proteins. Nucleic Acids Res. 2018, 46, W363–W367. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Sakkiah, S.; Guo, W.; Pan, B.; Ji, Z.; Yavas, G.; Azevedo, M.; Hawes, J.; Patterson, T.A.; Hong, H. Elucidating Interactions Between SARS-CoV-2 Trimeric Spike Protein and ACE2 Using Homology Modeling and Molecular Dynamics Simulations. Front. Chem. 2021, 8, 622632. [Google Scholar] [CrossRef] [PubMed]
  49. Sakkiah, S.; Kusko, R.; Pan, B.; Guo, W.; Ge, W.; Tong, W.; Hong, H. Structural Changes Due to Antagonist Binding in Ligand Binding Pocket of Androgen Receptor Elucidated Through Molecular Dynamics Simulations. Front. Pharmacol. 2018, 9, 492. [Google Scholar] [CrossRef]
  50. Selvaraj, C.; Sakkiah, S.; Tong, W.; Hong, H. Molecular Dynamics Simulations and Applications in Computational Toxicology and Nanotoxicology. Food Chem. Toxicol. 2018, 112, 495–506. [Google Scholar] [CrossRef]
  51. Shen, J.; Zhang, W.; Fang, H.; Perkins, R.; Tong, W.; Hong, H. Homology Modeling, Molecular Docking, and Molecular Dynamics Simulations Elucidated α-Fetoprotein Binding Modes. BMC Bioinform. 2013, 14, S6. [Google Scholar] [CrossRef] [Green Version]
  52. Shim, J.; Coop, A.; MacKerell, A.D. Molecular Details of the Activation of the μ Opioid Receptor. J. Phys. Chem. B 2013, 117, 7907–7917. [Google Scholar] [CrossRef] [Green Version]
  53. Yuan, S.; Vogel, H.; Filipek, S. The Role of Water and Sodium Ions in the Activation of the μ-Opioid Receptor. Angew. Chem. Int. Ed. 2013, 52, 10112–10115. [Google Scholar] [CrossRef]
  54. Bartuzi, D.; Kaczor, A.A.; Matosiuk, D. Activation and Allosteric Modulation of Human μ Opioid Receptor in Molecular Dynamics. J. Chem. Inf. Model. 2015, 55, 2421–2434. [Google Scholar] [CrossRef]
  55. Yuan, S.; Palczewski, K.; Peng, Q.; Kolinski, M.; Vogel, H.; Filipek, S. The Mechanism of Ligand-Induced Activation or Inhibition of μ- and κ-Opioid Receptors. Angew. Chem. Int. Ed. 2015, 54, 7560–7563. [Google Scholar] [CrossRef] [PubMed]
  56. Cheng, J.; Cheng, T.; Li, W.; Liu, G.; Zhu, W.; Tang, Y. Computational Insights into the G-Protein-Biased Activation and Inactivation Mechanisms of the μ Opioid Receptor. Acta Pharm. Sin. 2018, 39, 154–164. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Marino, K.A.; Shang, Y.; Filizola, M. Insights into the Function of Opioid Receptors from Molecular Dynamics Simulations of Available Crystal Structures. Br. J. Pharmacol. 2018, 175, 2834–2845. [Google Scholar] [CrossRef] [Green Version]
  58. Hu, X.; Wang, Y.; Hunkele, A.; Provasi, D.; Pasternak, G.W.; Filizola, M. Kinetic and Thermodynamic Insights into Sodium Ion Translocation through the μ-Opioid Receptor from Molecular Dynamics and Machine Learning Analysis. PLoS Comput. Biol. 2019, 15, e1006689. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Lipiński, P.F.J.; Jarończyk, M.; Dobrowolski, J.C.; Sadlej, J. Molecular Dynamics of Fentanyl Bound to μ-Opioid Receptor. J. Mol. Model. 2019, 25, 144. [Google Scholar] [CrossRef] [Green Version]
  60. Dumitrascuta, M.; Bermudez, M.; Ballet, S.; Wolber, G.; Spetea, M. Mechanistic Understanding of Peptide Analogues, DALDA, [Dmt1]DALDA, and KGOP01, Binding to the Mu Opioid Receptor. Molecules 2020, 25, 2087. [Google Scholar] [CrossRef]
  61. Hu, X.; Provasi, D.; Ramsey, S.; Filizola, M. Mechanism of μ-Opioid Receptor-Magnesium Interaction and Positive Allosteric Modulation. Biophys. J. 2020, 118, 909–921. [Google Scholar] [CrossRef]
  62. Mondal, D.; Kolev, V.; Warshel, A. Exploring the Activation Pathway and G i -Coupling Specificity of the μ-Opioid Receptor. Proc. Natl. Acad. Sci. USA 2020, 117, 26218–26225. [Google Scholar] [CrossRef]
  63. Zhao, Z.; Huang, T.; Li, J. Molecular Dynamics Simulations to Investigate How PZM21 Affects the Conformational State of the μ-Opioid Receptor Upon Activation. IJMS 2020, 21, 4699. [Google Scholar] [CrossRef]
  64. Podlewska, S.; Bugno, R.; Kudla, L.; Bojarski, A.J.; Przewlocki, R. Molecular Modeling of µ Opioid Receptor Ligands with Various Functional Properties: PZM21, SR-17018, Morphine, and Fentanyl—Simulated Interaction Patterns Confronted with Experimental Data. Molecules 2020, 25, 4636. [Google Scholar] [CrossRef]
  65. Liao, S.; Tan, K.; Floyd, C.; Bong, D.; Pino, M.J.; Wu, C. Probing Biased Activation of Mu-Opioid Receptor by the Biased Agonist PZM21 Using All Atom Molecular Dynamics Simulation. Life Sci. 2021, 269, 119026. [Google Scholar] [CrossRef]
  66. Xie, B.; Le Rouzic, V.P.; Goldberg, A.; Tsai, M.-H.M.; Chen, L.; Zhang, T.; Sinha, A.; Pan, Y.-X.; Baumann, M.H.; Shi, L. Binding Preference at the μ-Opioid Receptor Underlies Distinct Pharmacology of Cyclopropyl versus Valeryl Analogs of Fentanyl. Neuropharmacology 2023, 227, 109442. [Google Scholar] [CrossRef]
  67. Freier, E.; Wolf, S.; Gerwert, K. Proton Transfer via a Transient Linear Water-Molecule Chain in a Membrane Protein. Proc. Natl. Acad. Sci. USA 2011, 108, 11435–11439. [Google Scholar] [CrossRef] [Green Version]
  68. Fried, S.D.E.; Hewage, K.S.K.; Eitel, A.R.; Struts, A.V.; Weerasinghe, N.; Perera, S.M.D.C.; Brown, M.F. Hydration-Mediated G-Protein–Coupled Receptor Activation. Proc. Natl. Acad. Sci. USA 2022, 119, e2117349119. [Google Scholar] [CrossRef]
  69. Venkatakrishnan, A.J.; Deupi, X.; Lebon, G.; Tate, C.G.; Schertler, G.F.; Babu, M.M. Molecular Signatures of G-Protein-Coupled Receptors. Nature 2013, 494, 185–194. [Google Scholar] [CrossRef] [PubMed]
  70. Yuan, S.; Filipek, S.; Palczewski, K.; Vogel, H. Activation of G-Protein-Coupled Receptors Correlates with the Formation of a Continuous Internal Water Pathway. Nat. Commun. 2014, 5, 4733. [Google Scholar] [CrossRef] [Green Version]
  71. Weis, W.I.; Kobilka, B.K. The Molecular Basis of G Protein–Coupled Receptor Activation. Annu. Rev. Biochem. 2018, 87, 897–919. [Google Scholar] [CrossRef]
  72. Angel, T.E.; Chance, M.R.; Palczewski, K. Conserved Waters Mediate Structural and Functional Activation of Family A (Rhodopsin-like) G Protein-Coupled Receptors. Proc. Natl. Acad. Sci. USA 2009, 106, 8555–8560. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Raehal, K.M.; Walker, J.K.L.; Bohn, L.M. Morphine Side Effects in β-Arrestin 2 Knockout Mice. J. Pharm. Exp. 2005, 314, 1195–1201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Bohn, L.M.; Gainetdinov, R.R.; Lin, F.-T.; Lefkowitz, R.J.; Caron, M.G. μ-Opioid Receptor Desensitization by β-Arrestin-2 Determines Morphine Tolerance but Not Dependence. Nature 2000, 408, 720–723. [Google Scholar] [CrossRef]
  75. Bohn, L.M.; Lefkowitz, R.J.; Gainetdinov, R.R.; Peppel, K.; Caron, M.G.; Lin, F.-T. Enhanced Morphine Analgesia in Mice Lacking β-Arrestin 2. Science 1999, 286, 2495–2498. [Google Scholar] [CrossRef]
  76. Kliewer, A.; Gillis, A.; Hill, R.; Schmiedel, F.; Bailey, C.; Kelly, E.; Henderson, G.; Christie, M.J.; Schulz, S. Morphine-induced Respiratory Depression Is Independent of Β-arrestin2 Signalling. Br. J. Pharm. 2020, 177, 2923–2931. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Livingston, K.E.; Traynor, J.R. Allostery at Opioid Receptors: Modulation with Small Molecule Ligands. Br. J. Pharmacol. 2018, 175, 2846–2856. [Google Scholar] [CrossRef] [Green Version]
  78. Kathmann, M.; Flau, K.; Redmer, A.; Tränkle, C.; Schlicker, E. Cannabidiol Is an Allosteric Modulator at Mu- and Delta-Opioid Receptors. Naunyn. Schmied. Arch. Pharm. 2006, 372, 354–361. [Google Scholar] [CrossRef] [PubMed]
  79. Rothman, R.B.; Murphy, D.L.; Xu, H.; Godin, J.A.; Dersch, C.M.; Partilla, J.S.; Tidgewell, K.; Schmidt, M.; Prisinzano, T.E. Salvinorin A: Allosteric Interactions at the μ-Opioid Receptor. J. Pharm. Exp. 2007, 320, 801–810. [Google Scholar] [CrossRef] [Green Version]
  80. Fenalti, G.; Giguere, P.M.; Katritch, V.; Huang, X.-P.; Thompson, A.A.; Cherezov, V.; Roth, B.L.; Stevens, R.C. Molecular Control of δ-Opioid Receptor Signalling. Nature 2014, 506, 191–196. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Liu, W.; Chun, E.; Thompson, A.A.; Chubukov, P.; Xu, F.; Katritch, V.; Han, G.W.; Roth, C.B.; Heitman, L.H.; IJzerman, A.P.; et al. Structural Basis for Allosteric Regulation of GPCRs by Sodium Ions. Science 2012, 337, 232–236. [Google Scholar] [CrossRef] [Green Version]
  82. Miller-Gallacher, J.L.; Nehmé, R.; Warne, T.; Edwards, P.C.; Schertler, G.F.X.; Leslie, A.G.W.; Tate, C.G. The 2.1 Å Resolution Structure of Cyanopindolol-Bound Β1-Adrenoceptor Identifies an Intramembrane Na+ Ion That Stabilises the Ligand-Free Receptor. PLoS ONE 2014, 9, e92727. [Google Scholar] [CrossRef] [Green Version]
  83. Zhang, C.; Srinivasan, Y.; Arlow, D.H.; Fung, J.J.; Palmer, D.; Zheng, Y.; Green, H.F.; Pandey, A.; Dror, R.O.; Shaw, D.E.; et al. High-Resolution Crystal Structure of Human Protease-Activated Receptor 1. Nature 2012, 492, 387–392. [Google Scholar] [CrossRef] [Green Version]
  84. Shang, Y.; LeRouzic, V.; Schneider, S.; Bisignano, P.; Pasternak, G.W.; Filizola, M. Mechanistic Insights into the Allosteric Modulation of Opioid Receptors by Sodium Ions. Biochemistry 2014, 53, 5140–5149. [Google Scholar] [CrossRef] [Green Version]
  85. Vickery, O.N.; Carvalheda, C.A.; Zaidi, S.A.; Pisliakov, A.V.; Katritch, V.; Zachariae, U. Intracellular Transfer of Na+ in an Active-State G-Protein-Coupled Receptor. Structure 2018, 26, 171–180.e2. [Google Scholar] [CrossRef] [Green Version]
  86. Vickery, O.N.; Machtens, J.-P.; Tamburrino, G.; Seeliger, D.; Zachariae, U. Structural Mechanisms of Voltage Sensing in G Protein-Coupled Receptors. Structure 2016, 24, 997–1007. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Ye, L.; Neale, C.; Sljoka, A.; Lyda, B.; Pichugin, D.; Tsuchimura, N.; Larda, S.T.; Pomès, R.; García, A.E.; Ernst, O.P.; et al. Mechanistic Insights into Allosteric Regulation of the A2A Adenosine G Protein-Coupled Receptor by Physiological Cations. Nat. Commun. 2018, 9, 1372. [Google Scholar] [CrossRef]
  88. Qiu, Y.; Wang, Y.; Law, P.-Y.; Chen, H.-Z.; Loh, H.H. Cholesterol Regulates μ-Opioid Receptor-Induced β-Arrestin 2 Translocation to Membrane Lipid Rafts. Mol. Pharm. 2011, 80, 210–218. [Google Scholar] [CrossRef] [PubMed]
  89. Jakubík, J.; El-Fakahany, E.E. Allosteric Modulation of GPCRs of Class A by Cholesterol. IJMS 2021, 22, 1953. [Google Scholar] [CrossRef]
  90. Levitt, E.S.; Clark, M.J.; Jenkins, P.M.; Martens, J.R.; Traynor, J.R. Differential Effect of Membrane Cholesterol Removal on μ- and δ-Opioid Receptors. J. Biol. Chem. 2009, 284, 22108–22122. [Google Scholar] [CrossRef] [Green Version]
  91. Claff, T.; Yu, J.; Blais, V.; Patel, N.; Martin, C.; Wu, L.; Han, G.W.; Holleran, B.J.; Van der Poorten, O.; White, K.L.; et al. Elucidating the Active δ-Opioid Receptor Crystal Structure with Peptide and Small-Molecule Agonists. Sci. Adv. 2019, 5, eaax9115. [Google Scholar] [CrossRef] [Green Version]
  92. Che, T.; Majumdar, S.; Zaidi, S.A.; Ondachi, P.; McCorvy, J.D.; Wang, S.; Mosier, P.D.; Uprety, R.; Vardy, E.; Krumm, B.E.; et al. Structure of the Nanobody-Stabilized Active State of the Kappa Opioid Receptor. Cell 2018, 172, 55–67.e15. [Google Scholar] [CrossRef] [Green Version]
  93. Segala, E.; Guo, D.; Cheng, R.K.Y.; Bortolato, A.; Deflorian, F.; Doré, A.S.; Errey, J.C.; Heitman, L.H.; IJzerman, A.P.; Marshall, F.H.; et al. Controlling the Dissociation of Ligands from the Adenosine A 2A Receptor through Modulation of Salt Bridge Strength. J. Med. Chem. 2016, 59, 6470–6479. [Google Scholar] [CrossRef] [PubMed]
  94. Miles, T.F.; Spiess, K.; Jude, K.M.; Tsutsumi, N.; Burg, J.S.; Ingram, J.R.; Waghray, D.; Hjorto, G.M.; Larsen, O.; Ploegh, H.L.; et al. Viral GPCR US28 Can Signal in Response to Chemokine Agonists of Nearly Unlimited Structural Degeneracy. eLife 2018, 7, e35850. [Google Scholar] [CrossRef]
  95. Burford, N.T.; Clark, M.J.; Wehrman, T.S.; Gerritz, S.W.; Banks, M.; O’Connell, J.; Traynor, J.R.; Alt, A. Discovery of Positive Allosteric Modulators and Silent Allosteric Modulators of the μ-Opioid Receptor. Proc. Natl. Acad. Sci. USA 2013, 110, 10830–10835. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Livingston, K.E.; Traynor, J.R. Disruption of the Na + Ion Binding Site as a Mechanism for Positive Allosteric Modulation of the Mu-Opioid Receptor. Proc. Natl. Acad. Sci. USA 2014, 111, 18369–18374. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. Bartuzi, D.; Kaczor, A.A.; Matosiuk, D. Interplay between Two Allosteric Sites and Their Influence on Agonist Binding in Human μ Opioid Receptor. J. Chem. Inf. Model. 2016, 56, 563–570. [Google Scholar] [CrossRef]
  98. Bartuzi, D.; Kędzierska, E.; Kaczor, A.A.; Schmidhammer, H.; Matosiuk, D. Novel Positive Allosteric Modulators of µ Opioid Receptor—Insight from In Silico and In Vivo Studies. Int. J. Mol. Sci. 2020, 21, 8463. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The signal transduction pathway of the opioid receptor. Upon agonist binding, the opioid receptor couples with the heterotrimeric G protein that later dissociates into Gα and Gβγ subunits, followed by the translocation of the Gα subunit that leads to the inhibition of the adenylyl cyclase (AC) activity. The release of Gβγ subunit inhibits voltage-gated Ca2+ channels (VGCC, L-type, and N-type) and causes activation of K+ channels.
Figure 1. The signal transduction pathway of the opioid receptor. Upon agonist binding, the opioid receptor couples with the heterotrimeric G protein that later dissociates into Gα and Gβγ subunits, followed by the translocation of the Gα subunit that leads to the inhibition of the adenylyl cyclase (AC) activity. The release of Gβγ subunit inhibits voltage-gated Ca2+ channels (VGCC, L-type, and N-type) and causes activation of K+ channels.
Ijms 24 07042 g001
Figure 2. Biased signaling of MOR. The figure on the left shows the coupling between G proteins and MOR after ligand binding, leading to a reduction of cAMP, decreased Ca2+ signaling, and the activation of GIRK channels. In the middle, GRK phosphorylated the receptor, resulting in β-arrestin recruitment. This leads to receptor internalization and degradation. Meanwhile, β-arrestin mediates other signaling pathways, including those of MAPK and p38. The figure on the right shows the interaction between the Gαi subunit and β-arrestin. These two proteins form a complex that mediates various signaling pathways such as the activation of ERK.
Figure 2. Biased signaling of MOR. The figure on the left shows the coupling between G proteins and MOR after ligand binding, leading to a reduction of cAMP, decreased Ca2+ signaling, and the activation of GIRK channels. In the middle, GRK phosphorylated the receptor, resulting in β-arrestin recruitment. This leads to receptor internalization and degradation. Meanwhile, β-arrestin mediates other signaling pathways, including those of MAPK and p38. The figure on the right shows the interaction between the Gαi subunit and β-arrestin. These two proteins form a complex that mediates various signaling pathways such as the activation of ERK.
Ijms 24 07042 g002
Figure 3. Alignment of eight MOR structures: the two MOR structures in a complex with an antagonist are in different shades of blue and six MOR structures in a complex with an agonist are in different shades of yellow. The right panel is an enlarged illustration of all ligands (both agonist and antagonist) aligned. This figure was generated with ChimeraX 1.4. [44,45].
Figure 3. Alignment of eight MOR structures: the two MOR structures in a complex with an antagonist are in different shades of blue and six MOR structures in a complex with an agonist are in different shades of yellow. The right panel is an enlarged illustration of all ligands (both agonist and antagonist) aligned. This figure was generated with ChimeraX 1.4. [44,45].
Ijms 24 07042 g003
Figure 4. Binding pose comparison of antagonists β-FNA (top row left) and alvimopan (top row right), agonists BU72 (row 2 left), fentanyl (row 2 right), morphine (row 3 left), DAMGO (row 4 right), mitragynine pseudoindoxyl (row 5 left), and lofentanil (row 5 right), and partial agonists PZM21 (bottom row left), FH210 (bottom row right), oliceridine (TRV130) (row 3 right), and SR17018 (row 4 left). This figure was generated using Maestro. [46].
Figure 4. Binding pose comparison of antagonists β-FNA (top row left) and alvimopan (top row right), agonists BU72 (row 2 left), fentanyl (row 2 right), morphine (row 3 left), DAMGO (row 4 right), mitragynine pseudoindoxyl (row 5 left), and lofentanil (row 5 right), and partial agonists PZM21 (bottom row left), FH210 (bottom row right), oliceridine (TRV130) (row 3 right), and SR17018 (row 4 left). This figure was generated using Maestro. [46].
Ijms 24 07042 g004aIjms 24 07042 g004b
Figure 5. The common residues found within 4 Å for the antagonist (yellow), agonist (red), and partial agonist (light blue) are shown in licorice. The common residues found within 4 Å between the antagonist and agonist are colored orange, between the agonist and partial agonist are colored purple, and between the partial agonist and antagonist are colored green. The molecular image in this figure was generated with ChimeraX 1.4. [44,45].
Figure 5. The common residues found within 4 Å for the antagonist (yellow), agonist (red), and partial agonist (light blue) are shown in licorice. The common residues found within 4 Å between the antagonist and agonist are colored orange, between the agonist and partial agonist are colored purple, and between the partial agonist and antagonist are colored green. The molecular image in this figure was generated with ChimeraX 1.4. [44,45].
Ijms 24 07042 g005
Table 1. Four subtypes of opioid receptor and a list of their endogenous ligands.
Table 1. Four subtypes of opioid receptor and a list of their endogenous ligands.
Opioid Receptor SubtypeEndogenous Ligand
μβ-endorphin, enkephalins, endomorphin-1, endomorphin-2
δβ-endorphin, enkephalins
κdynorphin A, dynorphin B, α-neoendorphin
NOR 1N/OFQ (nociceptin/orphanin FQ)
1 nociceptin/orphanin FQ opioid peptide receptor.
Table 2. Available structures of μ opioid receptor.
Table 2. Available structures of μ opioid receptor.
PDB IDLigandLigand TypeIn Complex withStructure TypeResolution
4DKLβ-FNAAntagonistLysozyme chimeraX-RAY2.8
7UL4AlvimopanAntagonistMegabody 6Cryo-EM2.8
5C1MBU72Full agonistNanobody 39X-RAY2.07
8EF5FentanylFull agonistGαi-1, Gβ-1, Gγ-2Cryo-EM3.3
8EF6MorphineFull agonistGαi-1, Gβ-1, Gγ-2, scFv16Cryo-EM3.2
6DDEDAMGOFull agonistGαi-1, Gβ-1, Gγ-2, scFv16Cryo-EM3.5
6DDFDAMGOFull agonistGαi-1, Gβ-1, Gγ-2Cryo-EM3.5
8EFQDAMGOFull agonistGαi-1, Gβ-1, Gγ-2Cryo-EM3.3
7T2HLofentanilFull agonistGαi-1, Gβ-1, Gγ-2, scFv16Cryo-EM3.2
7T2GMitragynine pseudoindoxylFull agonistGαi-1, Gβ-1, Gγ-2Cryo-EM2.5
8EFBOliceridine
(TRV130)
Partial agonistGαi-1, Gβ-1, Gγ-2, scFv16Cryo-EM3.2
8EFLSR17018Partial agonistGαi-1, Gβ-1, Gγ-2, scFv16Cryo-EM3.2
8EFOPZM21Partial agonistGαi-1, Gβ-1, Gγ-2, scFv16Cryo-EM2.8
7SBFPZM21Partial agonistGαi-1, Gβ-1, Gγ-2, scFv16Cryo-EM2.9
7SCGFH210Partial agonistGαi-1, Gβ-1, Gγ-2, scFv16Cryo-EM3.0
7U2LNNPPP 1Bitopic ligandGαi-1, Gβ-1, Gγ-2, scFv16Cryo-EM3.2
1 N-(5-carbamimidamidopentyl)-N-[1-(2-phenylethyl) piperidin-4-yl] propanamide.
Table 3. The 2D structures of the ligands.
Table 3. The 2D structures of the ligands.
Ligand Type2D Structure
AgonistIjms 24 07042 i001
BU72
Ijms 24 07042 i002
DAMGO
Ijms 24 07042 i003
Morphine
Ijms 24 07042 i004
Fentanyl
Ijms 24 07042 i005
Lofentanil
Ijms 24 07042 i006
Mitragynine pseudoindoxyl
Partial agonistIjms 24 07042 i007
PZM21
Ijms 24 07042 i008
FH210
Ijms 24 07042 i009
Oliceridine (TRV130)
Ijms 24 07042 i010
SR17018
AntagonistIjms 24 07042 i011
β-FNA
Ijms 24 07042 i012
Alvimopan
Table 4. Binding pocket residue composition of each structure (residues within 5 Å of the ligand are included).
Table 4. Binding pocket residue composition of each structure (residues within 5 Å of the ligand are included).
PDB IDLigandCharged
Residues
Polar
Residues
Hydrophobic
Residues
4DKLβ-FNA528
7UL4Alvimopan3611
5C1MBU724410
8EF5Fentanyl1313
8EF6Morphine219
6DDEDAMGO4510
6DDF6DDEDAMGODAMGO44551010
6DDF8EFQDAMGO24551410
7T2HLofentanil2412
7T2GMitragynine pseudoindoxyl3510
8EFBOliceridine (TRV130)1315
8EFLSR170181413
8EFOPZM212413
7SBFPZM213511
7SCGFH210359
Table 5. Binding site solvent accessible surface area (SASA) and solvent accessible volume (SAV).
Table 5. Binding site solvent accessible surface area (SASA) and solvent accessible volume (SAV).
PDB IDLigandSASA (Å2)SAV (Å3)
4DKLβ-FNA724926
7UL4Alvimopan10161296
5C1MBU721212619
8EF5Fentanyl10621324
8EF6Morphine9601097
6DDEDAMGO716682
6DDFDAMGO631588
8EFQDAMGO827938
7T2HLofentanil8211024
7T2GMitragynine pseudoindoxyl612674
8EFBOliceridine
(TRV130)
10521067
8EFLSR1701812841875
8EFOPZM21949989
7SBFPZM2113412126
7SCGFH2109901330
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, Z.; Liu, J.; Dong, F.; Chang, N.; Huang, R.; Xia, M.; Patterson, T.A.; Hong, H. Three-Dimensional Structural Insights Have Revealed the Distinct Binding Interactions of Agonists, Partial Agonists, and Antagonists with the µ Opioid Receptor. Int. J. Mol. Sci. 2023, 24, 7042. https://doi.org/10.3390/ijms24087042

AMA Style

Li Z, Liu J, Dong F, Chang N, Huang R, Xia M, Patterson TA, Hong H. Three-Dimensional Structural Insights Have Revealed the Distinct Binding Interactions of Agonists, Partial Agonists, and Antagonists with the µ Opioid Receptor. International Journal of Molecular Sciences. 2023; 24(8):7042. https://doi.org/10.3390/ijms24087042

Chicago/Turabian Style

Li, Zoe, Jie Liu, Fan Dong, Nancy Chang, Ruili Huang, Menghang Xia, Tucker A. Patterson, and Huixiao Hong. 2023. "Three-Dimensional Structural Insights Have Revealed the Distinct Binding Interactions of Agonists, Partial Agonists, and Antagonists with the µ Opioid Receptor" International Journal of Molecular Sciences 24, no. 8: 7042. https://doi.org/10.3390/ijms24087042

APA Style

Li, Z., Liu, J., Dong, F., Chang, N., Huang, R., Xia, M., Patterson, T. A., & Hong, H. (2023). Three-Dimensional Structural Insights Have Revealed the Distinct Binding Interactions of Agonists, Partial Agonists, and Antagonists with the µ Opioid Receptor. International Journal of Molecular Sciences, 24(8), 7042. https://doi.org/10.3390/ijms24087042

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop