Molecular Modeling of µ Opioid Receptor Ligands with Various Functional Properties: PZM21, SR-17018, Morphine, and Fentanyl-Simulated Interaction Patterns Confronted with Experimental Data

Implementing physics-based digital patient twins to tailor the switch of oral morphine to transdermal fentanyl patches based on patient physiology

Fentanyl transdermal patches are widely implemented for cancer-induced pain treatment due to the high potency of fentanyl and gradual drug release. However, transdermal fentanyl up-titration for opioid-naïve patients is difficult, which is why opioid treatment is often started with oral/iv morphine. Based on the daily dose of morphine, the initial dose of the fentanyl patch is decided upon. After reaching a stable level of pain, the switch is made from oral/iv morphine to transdermal fentanyl. There are standard calculation tools for transferring from oral/iv morphine to transdermal fentanyl, which is the same for all patients. By considering the variations in the physiology of the patients, a unique switching strategy cannot meet the needs of different patients. This study explores the outcome in terms of pain relief and minute ventilation during opioid therapy. For this, we used physics-based simulations on a virtually-generated population of patients, and we applied the same therapy to all patients. We could show that patients' physiology, such as gender, age, and weight, greatly impact the outcome of the therapy; as such, the correlation coefficient between pain intensity and age is 0.89, and the correlation coefficient between patient's weight and maximum plasma concentration of morphine and fentanyl is -0.98 and -0.97. Additionally, a different combination of the duration of overlap between morphine and fentanyl therapy with different doses of fentanyl was considered for the virtual patients to find the best opioid-switching strategy for each patient. We explored the impact of combining physiological features to determine the best-suited strategy for virtual patients. Our findings suggest that tailoring morphine and fentanyl therapy only based on a limited number of features is insufficient, and increasing the number of impactful physiological features positively influences the outcome of the therapy.

Carfentanil is a β-arrestin-biased agonist at the μ opioid receptor

BACKGROUND AND PURPOSE

The illicit use of fentanyl-like drugs (fentanyls), which are μ opioid receptor agonists, and the many overdose deaths that result, has become a major problem. Fentanyls are very potent in vivo, leading to respiratory depression and death. However, the efficacy and possible signalling bias of different fentanyls is not clearly known. Here, we compared the relative efficacy and bias of a series of fentanyls.

EXPERIMENTAL APPROACH

For agonist signalling bias and efficacy measurements, Bioluminescence Resonance Energy Transfer experiments were undertaken in HEK293T cells transiently transfected with μ opioid receptors, to assess Gi protein activation and β-arrestin 2 recruitment. Agonist-induced cell surface receptor loss was assessed using an enzyme-linked immunosorbent assay, whilst agonist-induced G protein-coupled inwardly rectifying potassium channel current activation was measured electrophysiologically from rat locus coeruleus slices. Ligand poses in the μ opioid receptor were determined in silico using molecular dynamics simulations.

KEY RESULTS

Relative to the reference ligand DAMGO, carfentanil was β-arrestin-biased, whereas fentanyl, sufentanil and alfentanil did not display bias. Carfentanil induced potent and extensive cell surface receptor loss, whilst the marked desensitisation of G protein-coupled inwardly rectifying potassium channel currents in the continued presence of carfentanil in neurones was prevented by a GRK2/3 inhibitor. Molecular dynamics simulations suggested unique interactions of carfentanil with the orthosteric site of the receptor that could underlie the bias.

CONCLUSIONS AND IMPLICATIONS

Carfentanil is a β-arrestin-biased opioid drug at the μ receptor. It is uncertain how such bias influences in vivo effects of carfentanil relative to other fentanyls.

Dynamics of the Apo µ-Opioid Receptor in Complex with Gi Protein

Opioid receptors, particularly the µ-opioid receptor (μOR), play a pivotal role in mediating the analgesic and addictive effects of opioid drugs. G protein signaling is an important pathway of μOR function, usually associated with painkilling effects. However, the molecular mechanisms underlying the interaction between the μOR and G protein remain poorly understood. In this study, we employed classical all-atom molecular dynamics simulations to investigate the structural changes occurring with the μOR-G protein complex under two different conditions: with the G protein in the apo form (open) and with the GDP bound G protein (closed, holo form). The receptor was in the apo form and active conformation in both cases, and the simulation time comprised 1µs for each system. In order to assess the effect of the G protein coupling on the receptor activation state, three parameters were monitored: the correlation of the distance between TM3 and TM6 and the RMSD of the NPxxYA motif; the universal activation index (A100); and the χ2 dihedral distribution of residue W2936.48. When complexed with the open G protein, receptor conformations with intermediate activation state prevailed throughout the molecular dynamics, whereas in the condition with the closed G protein, mostly inactive conformations of the receptor were observed. The major effect of the G protein in the receptor conformation comes from a steric hindrance involving an intracellular loop of the receptor and a β-sheet region of the G protein. This suggests that G-protein precoupling is essential for receptor activation, but this fact is not sufficient for complete receptor activation.

N-substituted tetrahydro-beta-carboline as mu-opioid receptors ligands: in silico study; molecular docking, ADMET and molecular dynamics approach

Manipulating intracellular signals by interaction with transmembranal G-protein-coupled receptors (GPCRs) is the way of action of more than 30% of available medicines. Designing molecules against GPCRs is most challenging due to their flexible binding orthosteric and allosteric pockets, a property that lead to different mode and extent of activation of intracellular mediators. Here, in the current study we aimed to design N-substituted tetrahydro-beta-carbolines (THβC's) targeting Mu Opioid Receptors (MORs). We performed ligand docking study for reference and designed compounds against active and inactive states of MOR, as well as the active state bound to intracellular mediator of Gi. The reference compounds include 40 known agonists and antagonists, while the designed compounds include 25,227 N-substituted THβC analogues. Out of the designed compounds, 15 compounds were comparatively having better extra precision (XP) Gscore and were analyzed for absorption, distribution, metabolism, and excretion-toxicity (ADMET) properties, drug-likness, and molecular dynamic (MD) simulation. The results showed that N-substituted tetrahydro-beta-carbolines with and without C6-methoxy group substitutions (THBC/6MTHBC) analogues of A1/B1 and A9/B9 have relatively acceptable affinity and within pocket-stability toward MOR compared to the reference compounds of morphine (agonist) and naloxone (antagonist). Moreover, the designed analogues interact with key residue within the binding pocket of Asp 147 that is reported to be involved in receptor activation. In conclusion, the designed THBC analogues represent a good starting point for designing opioid receptor ligands other than morphinan scaffold, that have good synthetic accessibility which promotes feasible structural manipulation to tailor pharmacological effects with minimal side effects. Workflow rational in the discovery of potential Mu Opioid Receptor ligands.

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

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.

The anomalous pharmacology of fentanyl

Fentanyl is a key therapeutic, used in anaesthesia and pain management. It is also increasingly used illicitly and is responsible for a large and growing number of opioid overdose deaths, especially in North America. A number of factors have been suggested to contribute to fentanyl's lethality, including rapid onset of action, in vivo potency, ligand bias, induction of muscle rigidity and reduced sensitivity to reversal by naloxone. Some of these factors can be considered to represent 'anomalous' pharmacological properties of fentanyl when compared with prototypical opioid agonists such as morphine. In this review, we examine the nature of fentanyl's 'anomalous' properties, to determine whether there is really a pharmacological basis to support the existence of such properties, and also discuss whether such properties are likely to contribute to overdose deaths involving fentanyls. LINKED ARTICLES: This article is part of a themed issue on Advances in Opioid Pharmacology at the Time of the Opioid Epidemic. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v180.7/issuetoc.

Binding preference at the μ-opioid receptor underlies distinct pharmacology of cyclopropyl versus valeryl analogs of fentanyl

Illicitly manufactured fentanyl is driving the current opioid crisis, and various fentanyl analogs are appearing in recreational drug markets worldwide. To assess the potential health risks posed by fentanyl analogs, it is necessary to understand structure-activity relationships for these compounds. Here we compared the pharmacology of two structurally related fentanyl analogs implicated in opioid overdose: cyclopropylfentanyl and valerylfentanyl. Cyclopropylfentanyl has a three-carbon ring attached to the carbonyl group on the fentanyl scaffold, whereas valerylfentanyl has a four-carbon chain at the same position. In vitro assays examining μ-opioid receptor (MOR) coupling to G proteins in CHO cells showed that cyclopropylfentanyl is a full agonist (EC50 = 8.6 nM, %Emax = 113%), with potency and efficacy similar to fentanyl (EC50 = 10.3 nM, %Emax = 113%). By contrast, valerylfentanyl is a partial agonist at MOR (EC50 = 179.8 nM, %Emax = 60%). Similar results were found in assays assessing MOR-mediated β-arrestin recruitment in HEK cells. In vivo studies in male CD-1 mice demonstrated that both fentanyl analogs induce naloxone-reversible antinociception and respiratory suppression, but cyclopropylfentanyl is 100-times more potent as an antinociceptive agent (ED50 = 0.04 mg/kg, s. c.) than valerylfentanyl (ED50 = 4.0 mg/kg, s. c.). Molecular simulation results revealed that the alkyl chain of valerylfentanyl cannot be well accommodated by the active state of MOR and may transition the receptor toward an inactive state, converting the fentanyl scaffold to a partial agonist. Taken together, our results suggest that cyclopropylfentanyl presents much greater risk of adverse effects when compared to valerylfentanyl. Moreover, the summed findings may provide clues to the design of therapeutic opioids with reduced adverse side effects.

Isolation, characterization and in silico study of propenamide alkaloids from Hymenoepmecis bicolor poison against active μ-opioid receptor

Some insects produce venoms to defend against predators and directly interact with opioid receptors. In the present study, it was identified two alkaloids in the wasp venom species Hymenoepimecis bicolor. It was demonstrated that these could act as potential inhibitors of opioid receptors through their robust affinity to the receptors. The interaction profile was given to opioid receptors (μOR), with 60% of targets similar to alkaloid 1, with 0.25 probability, and 46.7% of targets similar to alkaloid 2, with a probability 0.17 of affinity as a target, which is considered signaling macromolecules and can mediate the most potent analgesic and addictive properties of opiate alkaloids. Notably, both alkaloids showed -7.6 kcal/mol affinity to the morphine agonies through six residues, Gly124, Asp147, Trp293, Ile296, Ile322, and Tyr326. These observations suggest further research on opioid receptors using in vitro studies of possible therapeutic applications.Communicated by Ramaswamy H. Sarma.

Demonstration of the Sequence Alignment to Predict Across Species Susceptibility Tool for Rapid Assessment of Protein Conservation

The US Environmental Protection Agency Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS) tool is a fast, freely available, online screening application that allows researchers and regulators to extrapolate toxicity information across species. For biological targets in model systems such as human cells, mice, rats, and zebrafish, toxicity data are available for a variety of chemicals. Through the evaluation of protein target conservation, this tool can be used to extrapolate data generated from such model systems to thousands of other species lacking toxicity data, yielding predictions of relative intrinsic chemical susceptibility. The latest releases of the tool (versions 2.0-6.1) have incorporated new features that allow for the rapid synthesis, interpretation, and use of the data for publication plus presentation-quality graphics. Among these features are customizable data visualizations and a comprehensive summary report designed to summarize SeqAPASS data for ease of interpretation. This paper describes the protocol to guide users through submitting jobs, navigating the various levels of protein sequence comparisons, and interpreting and displaying the resulting data. New features of SeqAPASS v2.0-6.0 are highlighted. Furthermore, two use-cases focused on transthyretin and opioid receptor protein conservation using this tool are described. Finally, SeqAPASS' strengths and limitations are discussed to define the domain of applicability for the tool and highlight different applications for cross-species extrapolation.

Biased Agonism: Lessons from Studies of Opioid Receptor Agonists

In ligand bias different agonist drugs are thought to produce distinct signaling outputs when activating the same receptor. If these signaling outputs mediate therapeutic versus adverse drug effects, then agonists that selectively activate the therapeutic signaling pathway would be extremely beneficial. It has long been thought that μ-opioid receptor agonists that selectively activate G protein- over β-arrestin-dependent signaling pathways would produce effective analgesia without the adverse effects such as respiratory depression. However, more recent data indicate that most of the therapeutic and adverse effects of agonist-induced activation of the μ-opioid receptor are actually mediated by the G protein-dependent signaling pathway, and that a number of drugs described as G protein biased in fact may not be biased, but instead may be low-intrinsic-efficacy agonists. In this review we discuss the current state of the field of bias at the μ-opioid receptor and other opioid receptor subtypes.

Recent Molecular Insights into Agonist-specific Binding to the Mu-Opioid Receptor

Opioid agonists produce their analgesic effects primarily by acting at the µ-opioid receptor (µOR). µOR agonists with different efficacies exert diverse molecular changes in the µOR which dictate the faith of the receptor’s signaling pathway and possibly it’s the degree of desensitization. Since the development of the active conformations of the µOR, growing data have been published in relation to ligand-specific changes in µOR activation. In this regard, this review summarizes recent data regarding the most studied opioid agonists in in silico µOR activation, including how these ligands are recognized by the µOR, how their binding signal is transmitted toward the intracellular parts of the µOR, and finally, what type of large-scale movements do these changes trigger in the µOR’s domains.

Endogenous opiates and behavior: 2020

This paper is the forty-third consecutive installment of the annual anthological review of research concerning the endogenous opioid system, summarizing articles published during 2020 that studied the behavioral effects of molecular, pharmacological and genetic manipulation of opioid peptides and receptors as well as effects of opioid/opiate agonists and antagonists. The review is subdivided into the following specific topics: molecular-biochemical effects and neurochemical localization studies of endogenous opioids and their receptors (1), the roles of these opioid peptides and receptors in pain and analgesia in animals (2) and humans (3), opioid-sensitive and opioid-insensitive effects of nonopioid analgesics (4), opioid peptide and receptor involvement in tolerance and dependence (5), stress and social status (6), learning and memory (7), eating and drinking (8), drug abuse and alcohol (9), sexual activity and hormones, pregnancy, development and endocrinology (10), mental illness and mood (11), seizures and neurologic disorders (12), electrical-related activity and neurophysiology (13), general activity and locomotion (14), gastrointestinal, renal and hepatic functions (15), cardiovascular responses (16), respiration and thermoregulation (17), and immunological responses (18).

A comprehensive evaluation of the potential binding poses of fentanyl and its analogs at the µ-opioid receptor

Fentanyl and its analogs are selective agonists of the µ-opioid receptor (MOR). Among novel synthetic opioids (NSOs), they dominate the recreational drug market and are the main culprits for the opioid crisis, which has been exacerbated by the COVID-19 pandemic. By taking advantage of the crystal structures of the MOR, several groups have investigated the binding mechanism of fentanyl, but have not reached a consensus, in terms of both the binding orientation and the fentanyl conformation. Thus, the binding mechanism of fentanyl at the MOR remains an unsolved and challenging question. Here, we carried out a systematic computational study to investigate the preferred fentanyl conformations, and how these conformations are being accommodated in the MOR binding pocket. We characterized the free energy landscape of fentanyl conformations with metadynamics simulations, and compared and evaluated several possible fentanyl binding conditions in the MOR with long-timescale molecular dynamics simulations. Our results indicate that the most preferred binding pose in the MOR binding pocket corresponds well with the global minimum on the energy landscape of fentanyl in the absence of the receptor, while the energy landscape can be reconfigured by modifying the fentanyl scaffold. The interactions with the receptor may stabilize a slightly unfavored fentanyl conformation in an alternative binding pose. By extending similar investigations to fentanyl analogs, our findings establish a structure–activity relationship of fentanyl binding at the MOR. In addition to providing a structural basis to understand the potential toxicity of the emerging NSOs, such insights will contribute to developing new, safer analgesics.

G protein signaling-biased mu opioid receptor agonists that produce sustained G protein activation are noncompetitive agonists

The ability of a ligand to preferentially promote engagement of one signaling pathway over another downstream of GPCR activation has been referred to as signaling bias, functional selectivity, and biased agonism. The presentation of ligand bias reflects selectivity between active states of the receptor, which may result in the display of preferential engagement with one signaling pathway over another. In this study, we provide evidence that the G protein-biased mu opioid receptor (MOR) agonists SR-17018 and SR-14968 stabilize the MOR in a wash-resistant yet antagonist-reversible G protein-signaling state. Furthermore, we demonstrate that these structurally related biased agonists are noncompetitive for radiolabeled MOR antagonist binding, and while they stimulate G protein signaling in mouse brains, partial agonists of this class do not compete with full agonist activation. Importantly, opioid antagonists can readily reverse their effects in vivo. Given that chronic treatment with SR-17018 does not lead to tolerance in several mouse pain models, this feature may be desirable for the development of long-lasting opioid analgesics that remain sensitive to antagonist reversal of respiratory suppression.

Low Intrinsic Efficacy Alone Cannot Explain the Improved Side Effect Profiles of New Opioid Agonists

In a recent report in Science Signaling (Gillis, A., et al. Low intrinsic efficacy for G protein activation can explain the improved side effect profiles of new opioid agonists. Sci. Signaling 2020, 13, eaaz3140 10.1126/scisignal.aaz3140), it was suggested that low intrinsic agonism, and not biased agonism, leads to an improvement in the separation of potency in opioid-induced respiratory suppression versus antinociception. Although many of the compounds that were tested have been shown to display G protein signaling bias in prior publications, the authors conclude that because they cannot detect biased agonism in their cellular signaling studies the compounds are therefore not biased agonists. Rather, they conclude that it is low intrinsic efficacy that leads to the therapeutic window improvement. Intrinsic efficacy is the extent to which an agonist can stimulate a G protein-coupled receptor response in a system, while biased agonism takes into consideration not only the intrinsic efficacy but also the potency of an agonist in an assay. Herein, we have reanalyzed the data presented in the published work (10.1126/scisignal.aaz3140) [including the recent Erratum (10.1126/scisignal.abf9803)] to derive intrinsic efficacy and bias factors as ΔΔlog(τ/K_A) and ΔΔlog(Emax/EC₅₀), respectively. On the basis of this reanalysis, the data support the conclusion that biased agonism, favoring G protein signaling, was observed. Moreover, a conservation of rank order intrinsic efficacy was not observed upon comparing responses in each assay, further suggesting that multiple active receptor states were present. These observations agree with prior studies in which oliceridine, PZM21, and SR-17018 were first described as biased agonists with improvement in antinociception over respiratory suppression in mice. Therefore, the data in the Science Signaling paper provide strong corroborating evidence that G protein signaling bias may be a means of improving opioid analgesia while avoiding certain undesirable side effects.

Opioids and Their Receptors: Present and Emerging Concepts in Opioid Drug Discovery

The interest in opioids such as morphine, the prototypical opioid ligand, has been maintained throughout the years [...].