How μ-opioid receptor recognizes fentanyl

Discovery of Novel, Selective, and Nonbasic Agonists for the Kappa-Opioid Receptor Determined by Salvinorin A-Based Virtual Screening

Modulating the kappa-opioid receptor (KOR) is a promising strategy for treating various human diseases. KOR agonists show potential for treating pain, pruritus, and epilepsy, while KOR antagonists show potential for treating depression, anxiety, and addiction. The diterpenoid Salvinorin A (SalA), a secondary metabolite of Salvia divinorum, is a potent and selective KOR agonist. Unlike typical opioids, SalA lacks a basic nitrogen, which encouraged us to search for nonbasic KOR ligands. Through structure-based virtual screening using 3D pharmacophore models based on the binding mode of SalA, we identified novel, nonbasic, potent, and selective KOR agonists. In vitro studies confirmed two virtual hits, SalA-VS-07 and SalA-VS-08, as highly selective for the KOR and showing G protein-biased KOR agonist activity. Both KOR ligands share a novel spiro-moiety and a nonbasic scaffold. Our findings provide novel starting points for developing therapeutics aimed at treating pain and other conditions in which KOR is a central player.

Balancing acts: The dual faces of fentanyl in medicine and public health

Fentanyl is a potent synthetic opioid widely used in medicine for its effective analgesic properties, particularly in surgical procedures and in the treatment of severe, chronic pain. In recent decades, however, there has been a worrying increase in the illicit use of fentanyl, particularly in North America. This rise in illicit use is concerning because fentanyl is associated with polydrug abuse, which adds layers of complexity and dangerous. This review provides a comprehensive examination of fentanyl, focusing on its synthesis and medical use. It also discusses the significance of the piperidine ring in medicinal chemistry as well as the critical role of fentanyl in pain management and anesthesia. Furthermore, it addresses the challenges associated with the abuse potential of fentanyl and the resulting public health concerns. The study aims to strike a balance between the clinical benefits and risks of fentanyl by advocating for innovative uses while addressing public health issues. It examines the chemistry, pharmacokinetics and pharmacodynamics of fentanyl and highlights the importance of personalized medicine in the administration of opioids. The review underscores the necessity of continuous research and adaptation in both clinical use and public health strategies.

Effects of pH on opioid receptor activation and implications for drug design

G-protein-coupled receptors are integral membrane proteins that transduce chemical signals from the extracellular matrix into the cell. Traditional drug design has considered ligand-receptor interactions only under normal conditions. However, studies on opioids indicate that such interactions are very different in diseased tissues. In such microenvironments, protons play an important role in structural and functional alterations of both ligands and receptors. The pertinent literature strongly suggests that future drug design should take these aspects into account in order to reduce adverse side effects while preserving desired effects of novel compounds.

Fentanyl activates opposing opioid and non-opioid receptor systems that control breathing

Fentanyl elicits profound disturbances in ventilatory control processes in humans and experimental animals. The traditional viewpoint with respect to fentanyl-induced respiratory depression is that once the effects on the frequency of breathing (Freq), tidal volume (TV), and minute ventilation (MV = Freq × TV) are resolved, then depression of breathing is no longer a concern. The results of the present study challenge this concept with findings, as they reveal that while the apparent inhibitory effects of fentanyl (75 μg/kg, IV) on Freq, TV, and MV in adult male rats were fully resolved within 15 min, many other fentanyl-induced responses were in full effect, including opposing effects on respiratory timing parameters. For example, although the effects on Freq were resolved at 15 min, inspiratory duration (Ti) and end inspiratory pause (EIP) were elevated, whereas expiratory duration (Te) and end expiratory pause (EEP) were diminished. Since the effects of fentanyl on TV had subsided fully at 15 min, it would be expected that the administration of an opioid receptor (OR) antagonist would have minimal effects if the effects of fentanyl on this and other parameters had resolved. We now report that the intravenous injection of a 1.0 mg/kg dose of the peripherally restricted OR antagonist, methyl-naloxone (naloxone methiodide, NLXmi), did not elicit arousal but elicited some relatively minor changes in Freq, TV, MV, Te, and EEP but pronounced changes in Ti and EIP. In contrast, the injection of a 2.5 mg/kg dose of NLXmi elicited pronounced arousal and dramatic changes in many variables, including Freq, TV, and MV, which were not associated with increases in non-apneic breathing events such as apneas. The two compelling conclusions from this study are as follows: 1) the blockade of central ORs produced by the 2.5 mg/kg dose of NLXmi elicits pronounced increases in Freq, TV, and MV in rats in which the effects of fentanyl had apparently resolved, and 2) it is apparent that fentanyl had induced the activation of two systems with counter-balancing effects on Freq and TV: one being an opioid receptor inhibitory system and the other being a non-OR excitatory system.

Identification of Unique Fragmentation Patterns of Fentanyl Analog Protomers Using Structures for Lossless Ion Manipulations Ion Mobility-Orbitrap Mass Spectrometry

The opioid crisis in the United States is being fueled by the rapid emergence of new fentanyl analogs and precursors that can elude traditional library-based screening methods, which require data from known reference compounds. Since reference compounds are unavailable for new fentanyl analogs, we examined if fentanyls (fentanyl + fentanyl analogs) could be identified in a reference-free manner using a combination of electrospray ionization (ESI), high-resolution ion mobility (IM) spectrometry, high-resolution mass spectrometry (MS), and higher-energy collision-induced dissociation (MS/MS). We analyzed a mixture containing nine fentanyls and W-15 (a structurally similar molecule) and found that the protonated forms of all fentanyls exhibited two baseline-separated IM distributions that produced different MS/MS patterns. Upon fragmentation, both IM distributions of all fentanyls produced two high intensity fragments, resulting from amine site cleavages. The higher mobility distributions of all fentanyls also produced several low intensity fragments, but surprisingly, these same fragments exhibited much greater intensities in the lower mobility distributions. This observation demonstrates that many fragments of fentanyls predominantly originate from one of two different gas-phase structures (suggestive of protomers). Furthermore, increasing the water concentration in the ESI solution increased the intensity of the lower mobility distribution relative to the higher mobility distribution, which further supports that fentanyls exist as two gas-phase protomers. Our observations on the IM and MS/MS properties of fentanyls can be exploited to positively differentiate fentanyls from other compounds without requiring reference libraries and will hopefully assist first responders and law enforcement in combating new and emerging fentanyls.

Role of Protonation States in the Stability of Molecular Dynamics Simulations of High-Resolution Membrane Protein Structures

Classical molecular dynamics (MD) simulations provide unmatched spatial and time resolution of protein structure and function. However, the accuracy of MD simulations often depends on the quality of force field parameters and the time scale of sampling. Another limitation of conventional MD simulations is that the protonation states of titratable amino acid residues remain fixed during simulations, even though protonation state changes coupled to conformational dynamics are central to protein function. Due to the uncertainty in selecting protonation states, classical MD simulations are sometimes performed with all amino acids modeled in their standard charged states at pH 7. Here, we performed and analyzed classical MD simulations on high-resolution cryo-EM structures of two large membrane proteins that transfer protons by catalyzing protonation/deprotonation reactions. In simulations performed with titratable amino acids modeled in their standard protonation (charged) states, the structure diverges far from its starting conformation. In comparison, MD simulations performed with predetermined protonation states of amino acid residues reproduce the structural conformation, protein hydration, and protein-water and protein-protein interactions of the structure much better. The results support the notion that it is crucial to perform basic protonation state calculations, especially on structures where protonation changes play an important functional role, prior to the launch of any conventional MD simulations. Furthermore, the combined approach of fast protonation state prediction and MD simulations can provide valuable information about the charge states of amino acids in the cryo-EM sample. Even though accurate prediction of protonation states in proteinaceous environments currently remains a challenge, we introduce an approach of combining pKa prediction with cryo-EM density map analysis that helps in improving not only the protonation state predictions but also the atomic modeling of density data.

Fentanyl dysregulates neuroinflammation and disrupts blood-brain barrier integrity in HIV-1 Tat transgenic mice

Opioid overdose deaths have dramatically increased by 781% from 1999 to 2021. In the setting of HIV, opioid drug abuse exacerbates neurotoxic effects of HIV in the brain, as opioids enhance viral replication, promote neuronal dysfunction and injury, and dysregulate an already compromised inflammatory response. Despite the rise in fentanyl abuse and the close association between opioid abuse and HIV infection, the interactive comorbidity between fentanyl abuse and HIV has yet to be examined in vivo. The HIV-1 Tat-transgenic mouse model was used to understand the interactive effects between fentanyl and HIV. Tat is an essential protein produced during HIV that drives the transcription of new virions and exerts neurotoxic effects within the brain. The Tat-transgenic mouse model uses a glial fibrillary acidic protein (GFAP)-driven tetracycline promoter which limits Tat production to the brain and this model is well used for examining mechanisms related to neuroHIV. After 7 days of fentanyl exposure, brains were harvested. Tight junction proteins, the vascular cell adhesion molecule, and platelet-derived growth factor receptor-β were measured to examine the integrity of the blood brain barrier. The immune response was assessed using a mouse-specific multiplex chemokine assay. For the first time in vivo, we demonstrate that fentanyl by itself can severely disrupt the blood-brain barrier and dysregulate the immune response. In addition, we reveal associations between inflammatory markers and tight junction proteins at the blood-brain barrier.

Molecular Pharmacology Profiling of Phenylfentanil and Its Analogues to Understand the Putative Involvement of an Adrenergic Mechanism in Fentanyl-Induced Respiratory Depression

While there are approved therapeutics to treat opioid overdoses, the need for treatments to reverse overdoses due to ultrapotent fentanyls remains unmet. This may be due in part to an adrenergic mechanism of fentanyls in addition to their stereotypical mu-opioid receptor (MOR) effects. Herein, we report our efforts to further understanding of the functions these distinct mechanisms impart. Employing the known MOR neutral antagonist phenylfentanil as a lead, 17 analogues were designed based on the concept of isosteric replacement. To probe mechanisms of action, these analogues were pharmacologically evaluated in vitro and in vivo, while in silico modeling studies were also conducted on phenylfentanil. While it did not indicate MOR involvement in vivo, phenylfentanil yielded respiratory minute volumes similar to those caused by fentanyl. Taken together with molecular modeling studies, these results indicated that respiratory effects of fentanyls may also correlate to inhibition of both α1A- and α1B-adrenergic receptors.

Molecular basis of opioid receptor signaling

Opioids are used for pain management despite the side effects that contribute to the opioid crisis. The pursuit of non-addictive opioid analgesics remains unattained due to the unresolved intricacies of opioid actions, receptor signaling cascades, and neuronal plasticity. Advancements in structural, molecular, and computational tools illuminate the dynamic interplay between opioids and opioid receptors, as well as the molecular determinants of signaling pathways, which are potentially interlinked with pharmacological responses. Here, we review the molecular basis of opioid receptor signaling with a focus on the structures of opioid receptors bound to endogenous peptides or pharmacological agents. These insights unveil specific interactions that dictate ligand selectivity and likely their distinctive pharmacological profiles. Biochemical analysis further unveils molecular features governing opioid receptor signaling. Simultaneously, the synergy between computational biology and medicinal chemistry continues to expedite the discovery of novel chemotypes with the promise of yielding more efficacious and safer opioid compounds.

Differential interaction patterns of opioid analgesics with µ opioid receptors correlate with ligand-specific voltage sensitivity

The µ opioid receptor (MOR) is the key target for analgesia, but the application of opioids is accompanied by several issues. There is a wide range of opioid analgesics, differing in their chemical structure and their properties of receptor activation and subsequent effects. A better understanding of ligand-receptor interactions and the resulting effects is important. Here, we calculated the respective binding poses for several opioids and analyzed interaction fingerprints between ligand and receptor. We further corroborated the interactions experimentally by cellular assays. As MOR was observed to display ligand-induced modulation of activity due to changes in membrane potential, we further analyzed the effects of voltage sensitivity on this receptor. Combining in silico and in vitro approaches, we defined discriminating interaction patterns responsible for ligand-specific voltage sensitivity and present new insights into their specific effects on activation of the MOR.

Comparing spinal anesthesia approaches for transurethral lithotripsy in patients with proximal ureteral stones: A randomized clinical trial of bupivacaine alone versus bupivacaine with fentanyl

Despite the benefits of spinal anesthesia and the preference of anesthesiologists for this technique, it is less accepted by urologists due to the proximity of the stone place in the ureter and the possibility of pain, restlessness, and occasional movements of the patient during surgery. The current study investigated the success of bupivacaine plus intrathecal fentanyl in patients undergoing transurethral lithotripsy (TUL). In this randomized clinical trial, from April 2021 to September 2021, 54 patients with proximal urolithiasis candidates for TUL were enrolled. Patients were randomly categorized into two groups: group A received bupivacaine 10 mg and 0.5 ml of normal saline, while group B received bupivacaine 10 mg with 0.5 ml (25µg) of intrathecal fentanyl. According to our findings, about 74% of the patients were men, and the mean age of the patients was 66.14±22.46 years. The onset time of the sensory block, sensory block level, pain score, degree of relaxation, depth of the motor block, occurrence of anesthesia complications, oxygen saturation, and mean arterial blood pressure were not significantly different between the two groups. However, the duration of the motor block in group B was longer than in group A (p<0.001). Also, retropulsion was observed only in five patients (18.5%) in group A, significantly higher than in group B (p=0.019). Bupivacaine with fentanyl 25µg provided adequate spinal anesthesia with lower retropulsion in patients with urolithiasis who are candidates for TUL.

Protons taken hostage: Dynamic H-bond networks of the pH-sensing GPR68

Proton-sensing G Protein Coupled Receptors (GPCRs) sense changes in the extracellular pH to effect cell signaling for cellular homeostasis. They tend to be overexpressed in solid tumors associated with acidic extracellular pH, and are of direct interest as drug targets. How proton-sensing GPCRs sense extracellular acidification and activate upon protonation change is important to understand, because it may guide the design of therapeutics. Lack of publicly available experimental structures make it challenging to discriminate between conflicting mechanisms proposed for proton-binding, as main roles have been assigned to either an extracellular histidine cluster or to an internal carboxylic triad. Here we present a protocol to derive and evaluate structural models of the proton-sensing GPR68. This approach integrates state-of-the-art homology modeling with microsecond-timescale atomistic simulations, and with a detailed assessment of the compatibility of the structural models with known structural features of class A GPCRs. To decipher structural elements of potential interest for protonation-coupled conformational changes of GPR68, we used the best-compatible model as a starting point for independent atomistic simulations of GPR68 with different protonation states, and graph computations to characterize the response of GPR68 to changes in protonation. We found that GPR68 hosts an extended hydrogen-bond network that inter-connects the extracellular histidine cluster to the internal carboxylic triad, and which can even reach groups at the cytoplasmic G-protein binding site. Taken together, results suggest that GPR68 relies on dynamic, hydrogen-bond networks to inter-connect extracellular and internal proton-binding sites, and to elicit conformational changes at the cytoplasmic G-protein binding site.

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.

Fentanyl and the Fluorinated Fentanyl Derivative NFEPP Elicit Distinct Hydrogen-Bond Dynamics of the Opioid Receptor

The development of safe therapeutics to manage pain is of central interest for biomedical applications. The fluorinated fentanyl derivative N-(3-fluoro-1-phenethylpiperidin-4-yl)-N-phenylpropionamide (NFEPP) is potentially a safer alternative to fentanyl because unlike fentanyl─which binds to the μ-opioid receptor (MOR) at both physiological and acidic pH─NFEPP might bind to the MOR only at acidic pH typical of inflamed tissue. Knowledge of the protonation-coupled dynamics of the receptor-drug interactions is thus required to understand the molecular mechanism by which receptor activation initiates cell signaling to silence pain. To this end, here we have carried out extensive atomistic simulations of the MOR in different protonation states, in the absence of opioid drugs, and in the presence of fentanyl vs NFEPP. We used graph-based analyses to characterize internal hydrogen-bond networks that could contribute to the activation of the MOR. We find that fentanyl and NFEPP prefer distinct binding poses and that, in their binding poses, fentanyl and NFEPP partake in distinct internal hydrogen-bond networks, leading to the cytoplasmic G-protein-binding region. Moreover, the protonation state of functionally important aspartic and histidine side chains impacts hydrogen-bond networks that extend throughout the receptor, such that the ligand-bound MOR presents at its cytoplasmic G-protein-binding side, a hydrogen-bonding environment where dynamics depend on whether fentanyl or NFEPP is bound, and on the protonation state of specific MOR groups. The exquisite sensitivity of the internal protein-water hydrogen-bond network to the protonation state and to details of the drug binding could enable the MOR to elicit distinct pH- and opioid-dependent responses at its cytoplasmic G-protein-binding site.

Involvement of the Opioid Peptide Family in Cancer Progression

Peptides mediate cancer progression favoring the mitogenesis, migration, and invasion of tumor cells, promoting metastasis and anti-apoptotic mechanisms, and facilitating angiogenesis/lymphangiogenesis. Tumor cells overexpress peptide receptors, crucial targets for developing specific treatments against cancer cells using peptide receptor antagonists and promoting apoptosis in tumor cells. Opioids exert an antitumoral effect, whereas others promote tumor growth and metastasis. This review updates the findings regarding the involvement of opioid peptides (enkephalins, endorphins, and dynorphins) in cancer development. Anticancer therapeutic strategies targeting the opioid peptidergic system and the main research lines to be developed regarding the topic reviewed are suggested. There is much to investigate about opioid peptides and cancer: basic information is scarce, incomplete, or absent in many tumors. This knowledge is crucial since promising anticancer strategies could be developed alone or in combination therapies with chemotherapy/radiotherapy.

Towards a scalable and controllable preparation of highly-uniform surface-enhanced Raman scattering substrates: Defect-free nanofilms as templates

The uniformity and reproducibility of substrates highly determine the applicability of surface-enhanced Raman scattering (SERS). Production of them, however, remains a challenge. Herein, we report a template-based strategy for the strictly controllable and handily scalable preparation of a very uniform SERS substrate, Ag nanoparticles (AgNPs)/nanofilm, where the template used is a flexible, transparent, self-standing, defect-free and robust nanofilm. Importantly, the obtained AgNPs/nanofilm is self-adhesive to surfaces of different properties and morphologies, ensuring in-situ and at real-time SERS detection. The enhancement factor (EF) of the substrate for rhodamine 6G (R6G) could reach 5.8 × 10¹⁰ with a detection limit (DL) of 1.0 × 10^-15 mol L^-1. Moreover, 500 bending tests and one-month storage showed no observable performance degradation, and up to 50.0 cm² scaled-up preparation depicted negligible effect upon the structure and the sensing performance. The real-life applicability of AgNPs/nanofilm was demonstrated by the sensitive detection of tetramethylthiuram disulfide on cherry tomato and fentanyl in methanol with a routine handheld Raman spectrometer. This work thus provides a reliable strategy for large area wet-chemical preparation of high-quality SERS substrates.

Catalytic Antibody Blunts Carfentanil-Induced Respiratory Depression

Carfentanil, the most potent of the fentanyl analogues, is at the forefront of synthetic opioid-related deaths, second to fentanyl. Moreover, the administration of the opioid receptor antagonist naloxone has proven inadequate for an increasing number of opioid-related conditions, often requiring higher/additional doses to be effective, as such interest in alternative strategies to combat more potent synthetic opioids has intensified. Increasing drug metabolism would be one strategy to detoxify carfentanil; however, carfentanil's major metabolic pathways involve N-dealkylation or monohydroxylation, which do not lend themselves readily to exogenous enzyme addition. Herein, we report, to our knowledge, the first demonstration that carfentanil's methyl ester when hydrolyzed to its acid was found to be 40,000 times less potent than carfentanil in activating the μ-opioid receptor. Physiological consequences of carfentanil and its acid were also examined through plethysmography, and carfentanil's acid was found to be incapable of inducing respiratory depression. Based upon this information, a hapten was chemically synthesized and immunized, allowing the generation of antibodies that were screened for carfentanil ester hydrolysis. From the screening campaign, three antibodies were found to accelerate the hydrolysis of carfentanil's methyl ester. From this series of catalytic antibodies, the most active underwent extensive kinetic analysis, allowing us to postulate its mechanism of hydrolysis against this synthetic opioid. In the context of potential clinical applications, the antibody, when passively administered, was able to reduce respiratory depression induced by carfentanil. The data presented supports further development of antibody catalysis as a biologic strategy to complement carfentanil overdose reversal.

A Molecular Toolbox of Positron Emission Tomography Tracers for General Anesthesia Mechanism Research

With appropriate radiotracers, positron emission tomography (PET) allows direct or indirect monitoring of the spatial and temporal distribution of anesthetics, neurotransmitters, and biomarkers, making it an indispensable tool for studying the general anesthesia mechanism. In this Perspective, PET tracers that have been recruited in general anesthesia research are introduced in the following order: 1) ¹¹C/¹⁸F-labeled anesthetics, i.e., PET tracers made from inhaled and intravenous anesthetics; 2) PET tracers targeting anesthesia-related receptors, e.g., neurotransmitters and voltage-gated ion channels; and 3) PET tracers for studying anesthesia-related neurophysiological effects and neurotoxicity. The radiosynthesis, pharmacodynamics, and pharmacokinetics of the above PET tracers are mainly discussed to provide a practical molecular toolbox for radiochemists, anesthesiologists, and those who are interested in general anesthesia.

Machine-learning repurposing of DrugBank compounds for opioid use disorder

Opioid use disorder (OUD) is a chronic and relapsing condition that involves the continued and compulsive use of opioids despite harmful consequences. The development of medications with improved efficacy and safety profiles for OUD treatment is urgently needed. Drug repurposing is a promising option for drug discovery due to its reduced cost and expedited approval procedures. Computational approaches based on machine learning enable the rapid screening of DrugBank compounds, identifying those with the potential to be repurposed for OUD treatment. We collected inhibitor data for four major opioid receptors and used advanced machine learning predictors of binding affinity that fuse the gradient boosting decision tree algorithm with two natural language processing (NLP)-based molecular fingerprints and one traditional 2D fingerprint. Using these predictors, we systematically analyzed the binding affinities of DrugBank compounds on four opioid receptors. Based on our machine learning predictions, we were able to discriminate DrugBank compounds with various binding affinity thresholds and selectivities for different receptors. The prediction results were further analyzed for ADMET (absorption, distribution, metabolism, excretion, and toxicity), which provided guidance on repurposing DrugBank compounds for the inhibition of selected opioid receptors. The pharmacological effects of these compounds for OUD treatment need to be tested in further experimental studies and clinical trials. Our machine learning studies provide a valuable platform for drug discovery in the context of OUD treatment.

Deciphering the intrinsic dynamics of unphosphorylated IRAK4 kinase bound to type I and type II inhibitors

Interleukin-1 receptor-associated kinase 4 (IRAK4) is a vital protein involved in Toll-like and interleukin-1 receptor signal transduction. Several studies have reported regarding the crystal structure, dynamic properties, and interactions with inhibitors of the phosphorylated form of IRAK4. However, no dynamic properties of inhibitor-bound unphosphorylated IRAK4 have been previously studied. Herein, we report the intrinsic dynamics of unphosphorylated IRAK4 (uIRAK4) bound to type I and type II inhibitors. The corresponding apo and inhibitor-bound forms of uIRAK4 were subjected to three independent simulations of 500 ns (total 1.5 μs) each, and their trajectories were analyzed. The results indicated that all three systems were relatively stable, except for the type II inhibitor-bound form of uIRAK4, which exhibited less compact folding and higher solvent surface area. The intra-hydrogen bonds corroborated the structural deformation of the type-II inhibitor-bound complex, which could be attributed to the long molecular structure of the type-II inhibitor. Moreover, the type II inhibitor bound to uIRAK4 showed higher binding free energy with uIRAK4 than the type I inhibitor. The free energy landscape analysis showed a reorientation of Phe330 side chain from the DFG motif at different metastable states for all the systems. The intra-residual distance between residues Lys213, Glu233, Tyr262, and Phe330 suggests a functional interplay when the inhibitors are bound to uIRAK4, thereby hinting at their crucial role in the inhibition mechanism. Ultimately, the intrinsic dynamics study observed between type I/II inhibitor-bound forms of uIRAK4 may assist in better understanding the enzyme and designing therapeutic compounds.

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.

Fentanyl inhibits cell invasion and migration by modulating NF-κB activation in glioma

Fentanyl is widely used for anesthesia and analgesia in cancer patients. Recent studies have revealed its anti-growth effect in several categories of cancer. Gliomas are the most common primary tumors in the central nervous system with poor prognosis. To investigate the effects of fentanyl on gliomas, glioma cells were treated with different concentrations of fentanyl both in vitro and in vivo. Consequences of proliferation and invasive phenotypes, and related protein expression were evaluated in two human glioma cell lines (U251 and U87). Naloxone, Mu Opioid Receptor (MOR) antagonist, was introduced into culture media to assess the involvement of MOR in Fentanyl-mediated changes. When compared with control group, it could be found that Fentanyl inhibited function of glioma cells only at high concentrations. Western blot and immunofluorescence results revealed that Fentanyl exerted its action via modulating NF-κB (P65) activation which is likely independent of MOR . Moreover, overexpression of P65 by transfection P65-expressing vector restored the invasion and migration of glioma cells, which were inhibited by Fentanyl. In summary, this study showed that opioid pain medication Fentanyl was capable of decreasing invasiveness of glioma cells at a high concentration both in vitro and in vivo, likely via modulating P65 activation.

Structure-Kinetics Relationships of Opioids from Metadynamics and Machine Learning Analysis

The nation's opioid overdose deaths reached an all-time high in 2021. The majority of deaths are due to synthetic opioids represented by fentanyl. Naloxone, which is a FDA-approved reversal agent, antagonizes opioids through competitive binding at the μ-opioid receptor (mOR). Thus, knowledge of the opioid's residence time is important for assessing the effectiveness of naloxone. Here, we estimated the residence times (τ) of 15 fentanyl and 4 morphine analogs using metadynamics and compared them with the most recent measurement of the opioid kinetic, dissociation, and naloxone inhibitory constants (Mann et al. Clin. Pharmacol. Therapeut. 2022, 120, 1020-1232). Importantly, the microscopic simulations offered a glimpse at the common binding mechanism and molecular determinants of dissociation kinetics for fentanyl analogs. The insights inspired us to develop a machine learning approach to analyze the kinetic impact of fentanyl's substituents based on the interactions with mOR residues. This proof-of-concept approach is general; for example, it may be used to tune ligand residence times in computer-aided drug discovery.

Endogenous opiates and behavior: 2021

This paper is the forty-fourth consecutive installment of the annual anthological review of research concerning the endogenous opioid system, summarizing articles published during 2021 that studied the behavioral effects of molecular, pharmacological and genetic manipulation of opioid peptides and receptors as well as effects of opioid/opiate agonizts 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).

Modified Akuamma Alkaloids with Increased Potency at the Mu-opioid Receptor

Akuammine (1) and pseudoakuammigine (2) are indole alkaloids found in the seeds of the akuamma tree (Picralima nitida). Both alkaloids are weak agonists of the mu opioid receptor (μOR); however, they produce minimal effects in animal models of antinociception. To probe the interactions of 1 and 2 at the opioid receptors, we have prepared a collection of 22 semisynthetic derivatives. Evaluation of this collection at the μOR and kappa opioid receptor (κOR) revealed structural-activity relationship trends and derivatives with improved potency at the μOR. Most notably, the introduction of a phenethyl moiety to the N1 of 2 produces a 70-fold increase in potency and a 7-fold increase in selectivity for the μOR. The in vitro potency of this compound resulted in increased efficacy in the tail-flick and hot-plate assays of antinociception. The improved potency of these derivatives highlights the promise of exploring natural product scaffolds to probe the opioid receptors.

Off-target activity of NBOMes and NBOMe analogs at the µ opioid receptor

New psychoactive substances (NPS) are introduced on the illicit drug market at a rapid pace. Their molecular targets are often inadequately elucidated, which contributes to the delayed characterization of their pharmacological effects. Inspired by earlier findings, this study set out to investigate the µ opioid receptor (MOR) activation potential of a large set of psychedelics, substances which typically activate the serotonin (5-HT_2A) receptor as their target receptor. We observed that some substances carrying the N-benzyl phenethylamine (NBOMe) structure activated MOR, as confirmed by both the NanoBiT® βarr2 recruitment assay and the G protein-based AequoScreen® Ca²⁺ release assay. The use of two orthogonal systems proved beneficial as some aspecific, receptor independent effects were found for various analogs when using the Ca²⁺ release assay. The specific 'off-target' effects at MOR could be blocked by the opioid antagonist naloxone, suggesting that these NBOMes occupy the same common opioid binding pocket as conventional opioids. This was corroborated by molecular docking, which revealed the plausibility of multiple interactions of 25I-NBOMe with MOR, similar to those observed for opioids. Additionally, structure-activity relationship findings seen in vitro were rationalized in silico for two 25I-NBOMe isomers. Overall, as MOR activity of these psychedelics was only noticed at high concentrations, we consider it unlikely that for the tested compounds there will be a relevant opioid toxicity in vivo at physiologically relevant concentrations. However, small modifications to the original NBOMe structure may result in a panel of more efficacious and potent MOR agonists, potentially exhibiting a dual MOR/5-HT_2A activation potential.

New science, drug regulation, and emergent public health issues: The work of FDA's division of applied regulatory science

The U.S. Food and Drug Administration (FDA) Division of Applied Regulatory Science (DARS) moves new science into the drug review process and addresses emergent regulatory and public health questions for the Agency. By forming interdisciplinary teams, DARS conducts mission-critical research to provide answers to scientific questions and solutions to regulatory challenges. Staffed by experts across the translational research spectrum, DARS forms synergies by pulling together scientists and experts from diverse backgrounds to collaborate in tackling some of the most complex challenges facing FDA. This includes (but is not limited to) assessing the systemic absorption of sunscreens, evaluating whether certain drugs can convert to carcinogens in people, studying drug interactions with opioids, optimizing opioid antagonist dosing in community settings, removing barriers to biosimilar and generic drug development, and advancing therapeutic development for rare diseases. FDA tasks DARS with wide ranging issues that encompass regulatory science; DARS, in turn, helps the Agency solve these challenges. The impact of DARS research is felt by patients, the pharmaceutical industry, and fellow regulators. This article reviews applied research projects and initiatives led by DARS and conducts a deeper dive into select examples illustrating the impactful work of the Division.

Solving an Old Puzzle: Elucidation and Evaluation of the Binding Mode of Salvinorin A at the Kappa Opioid Receptor

The natural product Salvinorin A (SalA) was the first nitrogen-lacking agonist discovered for the opioid receptors and exhibits high selectivity for the kappa opioid receptor (KOR) turning SalA into a promising analgesic to overcome the current opioid crisis. Since SalA's suffers from poor pharmacokinetic properties, particularly the absence of gastrointestinal bioavailability, fast metabolic inactivation, and subsequent short duration of action, the rational design of new tailored analogs with improved clinical usability is highly desired. Despite being known for decades, the binding mode of SalA within the KOR remains elusive as several conflicting binding modes of SalA were proposed hindering the rational design of new analgesics. In this study, we rationally determined the binding mode of SalA to the active state KOR by in silico experiments (docking, molecular dynamics simulations, dynophores) in the context of all available mutagenesis studies and structure-activity relationship (SAR) data. To the best of our knowledge, this is the first comprehensive evaluation of SalA's binding mode since the determination of the active state KOR crystal structure. SalA binds above the morphinan binding site with its furan pointing toward the intracellular core while the C2-acetoxy group is oriented toward the extracellular loop 2 (ECL2). SalA is solely stabilized within the binding pocket by hydrogen bonds (C210ECL2, Y3127.35, Y3137.36) and hydrophobic contacts (V1182.63, I1393.33, I2946.55, I3167.39). With the disruption of this interaction pattern or the establishment of additional interactions within the binding site, we were able to rationalize the experimental data for selected analogs. We surmise the C2-substituent interactions as important for SalA and its analogs to be experimentally active, albeit with moderate frequency within MD simulations of SalA. We further identified the non-conserved residues 2.63, 7.35, and 7.36 responsible for the KOR subtype selectivity of SalA. We are confident that the elucidation of the SalA binding mode will promote the understanding of KOR activation and facilitate the development of novel analgesics that are urgently needed.

Impact of fentanyl analgesia on the accuracy of HVPG measurements in patients with cirrhosis: a prospective, multicenter study

BACKGROUND

HVPG measurement is the gold standard for assessing portal hypertension. Many patients decline HVPG measurements due to associated pain. According to previous studies, propofol sedation during HVPG measurements potentially alters HVPG readings. However, opioid analgesics' effects on HVPG await full elucidation. This study aimed to evaluate fentanyl analgesia's effects on HVPG measurement accuracy in patients with cirrhosis.

METHODS

This prospective, multicenter study included patients with cirrhosis undergoing HVPG measurements, which were performed preanalgesia and under analgesia with fentanyl injection (1.0 or 1.5 μg/kg).

RESULTS

Of the 48 enrolled patients with cirrhosis, 23 were administered 1.0 μg/kg fentanyl analgesia during HVPG measurement. The HVPG was 13.4±4.9 mm Hg in preanalgesia and 13.5±5.2 mm Hg under analgesia. HVPG measurement accuracy was not altered after fentanyl analgesia (p = 0.801). The following measures also did not change: heart rate (p = 0.132), mean arterial pressure (p = 0.348), and blood oxygen saturation (p = 0.748); however, respiratory rate (p = 0.001) changes occurred. The Verbal Numerical Rating Score for comfort under analgesia was higher than that in preanalgesia (p = 0.001). Twenty-five patients were administered 1.5 μg/kg fentanyl analgesia during HVPG measurement. The HVPG was 19.5±5.7 mm Hg in preanalgesia and 19.6±5.6 mm Hg under analgesia. HVPG measurement accuracy did not alter after fentanyl analgesia (p = 0.469). Similarly, the following measures did not change: mean arterial pressure (p = 0.871) and oxygen saturation (p = 0.327); nevertheless, respiratory rate (p = 0.015) and heart rate (p = 0.019) changes occurred. The Verbal Numerical Rating Score for comfort under analgesia was higher than that in preanalgesia (p < 0.001).

CONCLUSION

Fentanyl analgesia did not alter HVPG measurement accuracy, and fentanyl improved comfort in patients with cirrhosis during HVPG measurements.

Metadynamics simulations leveraged by statistical analyses and artificial intelligence-based tools to inform the discovery of G protein-coupled receptor ligands

G Protein-Coupled Receptors (GPCRs) are a large family of membrane proteins with pluridimensional signaling profiles. They undergo ligand-specific conformational changes, which in turn lead to the differential activation of intracellular signaling proteins and the consequent triggering of a variety of biological responses. This conformational plasticity directly impacts our understanding of GPCR signaling and therapeutic implications, as do ligand-specific kinetic differences in GPCR-induced transducer activation/coupling or GPCR-transducer complex stability. High-resolution experimental structures of ligand-bound GPCRs in the presence or absence of interacting transducers provide important, yet limited, insights into the highly dynamic process of ligand-induced activation or inhibition of these receptors. We and others have complemented these studies with computational strategies aimed at characterizing increasingly accurate metastable conformations of GPCRs using a combination of metadynamics simulations, state-of-the-art algorithms for statistical analyses of simulation data, and artificial intelligence-based tools. This minireview provides an overview of these approaches as well as lessons learned from them towards the identification of conformational states that may be difficult or even impossible to characterize experimentally and yet important to discover new GPCR ligands.

GPU-Accelerated All-Atom Particle-Mesh Ewald Continuous Constant pH Molecular Dynamics in Amber

Constant pH molecular dynamics (MD) simulations sample protonation states on the fly according to the conformational environment and user specified pH conditions; however, the current accuracy is limited due to the use of implicit-solvent models or a hybrid solvent scheme. Here, we report the first GPU-accelerated implementation, parametrization, and validation of the all-atom continuous constant pH MD (CpHMD) method with particle-mesh Ewald (PME) electrostatics in the Amber22 pmemd.cuda engine. The titration parameters for Asp, Glu, His, Cys, and Lys were derived for the CHARMM c22 and Amber ff14sb and ff19sb force fields. We then evaluated the PME-CpHMD method using the asynchronous pH replica-exchange titration simulations with the c22 force field for six benchmark proteins, including BBL, hen egg white lysozyme (HEWL), staphylococcal nuclease (SNase), thioredoxin, ribonuclease A (RNaseA), and human muscle creatine kinase (HMCK). The root-mean-square deviation from the experimental pK_a's of Asp, Glu, His, and Cys is 0.76 pH units, and the Pearson's correlation coefficient for the pK_a shifts with respect to model values is 0.80. We demonstrated that a finite-size correction or much enlarged simulation box size can remove a systematic error of the calculated pK_a's and improve agreement with experiment. Importantly, the simulations captured the relevant biology in several challenging cases, e.g., the titration order of the catalytic dyad Glu35/Asp52 in HEWL and the coupled residues Asp19/Asp21 in SNase, the large pK_a upshift of the deeply buried catalytic Asp26 in thioredoxin, and the large pK_a downshift of the deeply buried catalytic Cys283 in HMCK. We anticipate that PME-CpHMD will offer proper pH control to improve the accuracies of MD simulations and enable mechanistic studies of proton-coupled dynamical processes that are ubiquitous in biology but remain poorly understood due to the lack of experimental tools and limitation of current MD simulations.

Constant pH molecular dynamics simulations: Current status and recent applications

Many important protein functions are carried out through proton-coupled conformational dynamics. Thus, the ability to accurately model protonation states dynamically has wide-ranging implications. Over the past two decades, two main types of constant pH methods (discrete and continuous) have been developed to enable proton-coupled molecular dynamics (MD) simulations. In this short review, we discuss the current status of the development and highlight recent applications that have advanced our understanding of protein structure-function relationships. We conclude the review by outlining the remaining challenges in the method development and projecting important areas for future applications.

Photochromic Fentanyl Derivatives for Controlled μ-Opioid Receptor Activation

Photoswitchable ligands as biological tools provide an opportunity to explore the kinetics and dynamics of the clinically relevant μ-opioid receptor. These ligands can potentially activate or deactivate the receptor when desired by using light. Spatial and temporal control of biological activity allows for application in a diverse range of biological investigations. Photoswitchable ligands have been developed in this work, modelled on the known agonist fentanyl, with the aim of expanding the current "toolbox" of fentanyl photoswitchable ligands. In doing so, ligands have been developed that change geometry (isomerize) upon exposure to light, with varying photophysical and biochemical properties. This variation in properties could be valuable in further studying the functional significance of the μ-opioid receptor.

Molecular recognition of morphine and fentanyl by the human μ-opioid receptor

Morphine and fentanyl are among the most used opioid drugs that confer analgesia and unwanted side effects through both G protein and arrestin signaling pathways of μ-opioid receptor (μOR). Here, we report structures of the human μOR-G protein complexes bound to morphine and fentanyl, which uncover key differences in how they bind the receptor. We also report structures of μOR bound to TRV130, PZM21, and SR17018, which reveal preferential interactions of these agonists with TM3 side of the ligand-binding pocket rather than TM6/7 side. In contrast, morphine and fentanyl form dual interactions with both TM3 and TM6/7 regions. Mutations at the TM6/7 interface abolish arrestin recruitment of μOR promoted by morphine and fentanyl. Ligands designed to reduce TM6/7 interactions display preferential G protein signaling. Our results provide crucial insights into fentanyl recognition and signaling of μOR, which may facilitate rational design of next-generation analgesics.

Fentanyl induces autism-like behaviours in mice by hypermethylation of the glutamate receptor gene Grin2b

BACKGROUND

Environmental factors contribute to autism spectrum disorder. Fentanyl, one of the most widely used opioid analgesics in anaesthesia, can induce neurotoxicity, but its role in autism remains unknown. We determined whether fentanyl induced autism-like behaviours in young mice and the underlying mechanisms.

METHODS

Young male and female mice received fentanyl at postnatal days 6, 8, and 10, and performed behavioural tests, including three-chamber social preference, elevated plus maze, grooming behaviour, and open-field test, from postnatal days 30-32. Expression of Grin2b, the gene encoding the GluN2B subunit of the N-methyl-d-aspartate receptor, was assessed in the anterior cingulate cortex of male mice using fluorescence in situ hybridisation histochemistry. We used bisulfite target sequencing to determine Grin2b hypermethylation sites after fentanyl treatment. In the specific activation and rescue experiments, we injected the mu opioid receptor agonist [D-Ala,² N-MePhe,⁴ Gly-ol]-enkephalin (DAMGO) or Grin2b overexpression lentivirus into the anterior cingulate cortex of male mice.

RESULTS

Fentanyl induced autism-like behaviours in both young male and female mice, and downregulated Grin2b expression (0.49-fold [0.08] vs 1.00-fold [0.09]; P<0.01) and GluN2B protein amounts (0.38-fold [0.07] vs 1.00-fold [0.12]; P<0.01) in the anterior cingulate cortex through hypermethylation of Grin2b. The mu-opioid receptor antagonist naloxone and overexpression of Grin2b in anterior cingulate cortex attenuated the fentanyl-induced effects, whereas DAMGO injection into the anterior cingulate cortex induced autism-like behaviours.

CONCLUSIONS

These data suggest that fentanyl induces autism-like behaviours in young mice via an epigenetic mechanism. Further research is required to determine possible clinical relevance to autism risk.

A Guide to the Continuous Constant pH Molecular Dynamics Methods in Amber and CHARMM [Article v1.0]

Like temperature and pressure, solution pH is an important environmental variable in biomolecular simulations. Virtually all proteins depend on pH to maintain their structure and function. In conventional molecular dynamics (MD) simulations of proteins, pH is implicitly accounted for by assigning and fixing protonation states of titratable sidechains. This is a significant limitation, as the assigned protonation states may be wrong and they may change during dynamics. In this tutorial, we guide the reader in learning and using the various continuous constant pH MD methods in Amber and CHARMM packages, which have been applied to predict pK a values and elucidate proton-coupled conformational dynamics of a variety of proteins including enzymes and membrane transporters.

Changes in Protonation States of In-Pathway Residues can Alter Ligand Binding Pathways Obtained From Spontaneous Binding Molecular Dynamics Simulations

Protein-ligand binding processes often involve changes in protonation states that can be key to recognize and orient the ligand in the binding site. The pathways through which (bio)molecules interplay to attain productively bound complexes are intricate and involve a series of interconnected intermediate and transition states. Molecular dynamics (MD) simulations and enhanced sampling techniques are commonly used to characterize the spontaneous binding of a ligand to its receptor. However, the effect of protonation state changes of in-pathway residues in spontaneous binding MD simulations remained mostly unexplored. Here, we used molecular dynamics simulations to reconstruct the trypsin-benzamidine binding pathway considering different protonation states of His57. This residue is part of the trypsin catalytic triad and is located more than 10 Å away from Asp189, which is responsible for benzamidine binding in the trypsin S1 pocket. Our MD simulations showed that the binding pathways that benzamidine follow to target the S1 binding site are critically dependent on the His57 protonation state. Binding of benzamidine frequently occurs when His57 is protonated in the delta nitrogen while the binding process is significantly less frequent when His57 is positively charged. Constant-pH MD simulations retrieved the equilibrium populations of His57 protonation states at trypsin active pH offering a clearer picture of benzamidine recognition and binding. These results indicate that properly accounting for protonation states of distal residues can be important in spontaneous binding MD simulations.

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.

In Vitro, In Vivo and In Silico Characterization of a Novel Kappa-Opioid Receptor Antagonist

Kappa-opioid receptor (KOR) antagonists are promising innovative therapeutics for the treatment of the central nervous system (CNS) disorders. The new scaffold opioid ligand, Compound A, was originally found as a mu-opioid receptor (MOR) antagonist but its binding/selectivity and activation profile at the KOR and delta-opioid receptor (DOR) remain elusive. In this study, we present an in vitro, in vivo and in silico characterization of Compound A by revealing this ligand as a KOR antagonist in vitro and in vivo. In the radioligand competitive binding assay, Compound A bound at the human KOR, albeit with moderate affinity, but with increased affinity than to the human MOR and without specific binding at the human DOR, thus displaying a preferential KOR selectivity profile. Following subcutaneous administration in mice, Compound A effectively reverse the antinociceptive effects of the prototypical KOR agonist, U50,488. In silico investigations were carried out to assess the structural determinants responsible for opioid receptor subtype selectivity of Compound A. Molecular docking, molecular dynamics simulations and dynamic pharmacophore (dynophore) generation revealed differences in the stabilization of the chlorophenyl moiety of Compound A within the opioid receptor binding pockets, rationalizing the experimentally determined binding affinity values. This new chemotype bears the potential for favorable ADMET properties and holds promise for chemical optimization toward the development of potential therapeutics.

Unique Pharmacology, Brain Dysfunction, and Therapeutic Advancements for Fentanyl Misuse and Abuse

Fentanyl is a fully synthetic opioid with analgesic and anesthetic properties. It has become a primary driver of the deadliest opioid crisis in the United States and elsewhere, consequently imposing devastating social, economic, and health burdens worldwide. However, the neural mechanisms that underlie the behavioral effects of fentanyl and its analogs are largely unknown, and approaches to prevent fentanyl abuse and fentanyl-related overdose deaths are scarce. This review presents the abuse potential and unique pharmacology of fentanyl and elucidates its potential mechanisms of action, including neural circuit dysfunction and neuroinflammation. We discuss recent progress in the development of pharmacological interventions, anti-fentanyl vaccines, anti-fentanyl/heroin conjugate vaccines, and monoclonal antibodies to attenuate fentanyl-seeking and prevent fentanyl-induced respiratory depression. However, translational studies and clinical trials are still lacking. Considering the present opioid crisis, the development of effective pharmacological and immunological strategies to prevent fentanyl abuse and overdose are urgently needed.

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.

Computational analyses of the vibrational spectra of fentanyl, carfentanil and remifentanil

The infrared (IR) spectra of fentanyl, carfentanil and remifentanil, and protonated salts, are computed using quantum chemistry methods. New experimental FTIR spectra are also reported and compared to the calculations. The accuracy of two density functional theory methods, B3LYP and M06-2X, are tested against higher level theories (MP2) and the experimental data. Gas phase IR spectra are calculated for both the neutral and protonated molecules in order to compare with the experimental data measured for various salts of fentanyl and its analogues. Key vibrational modes are selected and studied in detail using a vibrational mode locality calculation. The main contributing atomic movements in these vibrational modes are identified.

Fentanyl Structure as a Scaffold for Opioid/Non-Opioid Multitarget Analgesics

One of the strategies in the search for safe and effective analgesic drugs is the design of multitarget analgesics. Such compounds are intended to have high affinity and activity at more than one molecular target involved in pain modulation. In the present contribution we summarize the attempts in which fentanyl or its substructures were used as a μ-opioid receptor pharmacophoric fragment and a scaffold to which fragments related to non-opioid receptors were attached. The non-opioid ‘second’ targets included proteins as diverse as imidazoline I₂ binding sites, CB₁ cannabinoid receptor, NK₁ tachykinin receptor, D₂ dopamine receptor, cyclooxygenases, fatty acid amide hydrolase and monoacylglycerol lipase and σ₁ receptor. Reviewing the individual attempts, we outline the chemistry, the obtained pharmacological properties and structure-activity relationships. Finally, we discuss the possible directions for future work.

Mechanistic Characterization of the Pharmacological Profile of HS-731, a Peripherally Acting Opioid Analgesic, at the µ-, δ-, κ-Opioid and Nociceptin Receptors

Accumulated preclinical and clinical data show that peripheral restricted opioids provide pain relief with reduced side effects. The peripherally acting opioid analgesic HS-731 is a potent dual μ-/δ-opioid receptor (MOR/DOR) full agonist, and a weak, partial agonist at the κ-opioid receptor (KOR). However, its binding mode at the opioid receptors remains elusive. Here, we present a comprehensive in silico evaluation of HS-731 binding at all opioid receptors. We provide insights into dynamic interaction patterns explaining the different binding and activity of HS-731 on the opioid receptors. For this purpose, we conducted docking, performed molecular dynamics (MD) simulations and generated dynamic pharmacophores (dynophores). Our results highlight two residues important for HS-731 recognition at the classical opioid receptors (MOR, DOR and KOR), particular the conserved residue 5.39 (K) and the non-conserved residue 6.58 (MOR: K, DOR: W and KOR: E). Furthermore, we assume a salt bridge between the transmembrane helices (TM) 5 and 6 via K2275.39 and E2976.58 to be responsible for the partial agonism of HS-731 at the KOR. Additionally, we experimentally demonstrated the absence of affinity of HS-731 to the nociceptin/orphanin FQ peptide (NOP) receptor. We consider the morphinan phenol Y1303.33 responsible for this affinity lack. Y1303.33 points deep into the NOP receptor binding pocket preventing HS-731 binding to the orthosteric binding pocket. These findings provide significant structural insights into HS-731 interaction pattern with the opioid receptors that are important for understanding the pharmacology of this peripheral opioid analgesic.

Substance abuse and neurotransmission

The number of people who suffer from a substance abuse disorder has continued to rise over the last decade; particularly, the number of drug-related overdose deaths has sharply increased during the COVID-19 pandemic. Converging lines of clinical observations, supported by imaging and neuropsychological performance testing, have demonstrated that substance abuse-induced dysregulation of neurotransmissions in the brain is critical for development and expression of the addictive properties of abused substances. Recent scientific advances have allowed for better understanding of the neurobiological processes that mediates drugs of abuse and addiction. This chapter presents the past classic concepts and the recent advances in our knowledge about how cocaine, amphetamines, opioids, alcohol, and nicotine alter multiple neurotransmitter systems, which contribute to the behaviors associated with each drug. Additionally, we discuss the interactive effects of HIV-1 or COVID-19 and substance abuse on neurotransmission and neurobiological pathways. Finally, we introduce therapeutic strategies for development of pharmacotherapies for substance abuse disorders.

Kinetics and Mechanism of Fentanyl Dissociation from the μ-Opioid Receptor

Driven by illicit fentanyl, opioid related deaths have reached the highest level in 2020. Currently, an opioid overdose is resuscitated by the use of naloxone, which competitively binds and antagonizes the μ-opioid receptor (mOR). Thus, knowledge of the residence times of opioids at mOR and the unbinding mechanisms is valuable for assessing the effectiveness of naloxone. In the present study, we calculate the fentanyl-mOR dissociation time and elucidate the mechanism by applying an enhanced sampling molecular dynamics (MD) technique. Two sets of metadynamics simulations with different initial structures were performed while accounting for the protonation state of the conserved H297^6.52, which has been suggested to modulate the ligand-mOR affinity and binding mode. Surprisingly, with the Nδ-protonated H297^6.52, fentanyl can descend as much as 10 Å below the level of the conserved D147^3.32 before escaping the receptor and has a calculated residence time τ of 38 s. In contrast, with the Nϵ- and doubly protonated H297^6.52, the calculated τ are 2.6 and 0.9 s, respectively. Analysis suggests that formation of the piperidine-Hid297 hydrogen bond strengthens the hydrophobic contacts with the transmembrane helix (TM) 6, allowing fentanyl to explore a deep pocket. Considering the experimental τ of ∼4 min for fentanyl and the role of TM6 in mOR activation, the deep insertion mechanism may be biologically relevant. The work paves the way for large-scale computational predictions of opioid dissociation rates to inform evaluation of strategies for opioid overdose reversal. The profound role of the histidine protonation state found here may shift the paradigm in computational studies of ligand-receptor kinetics.

Protein pK a Prediction with Machine Learning

Protein pK _a prediction is essential for the investigation of the pH-associated relationship between protein structure and function. In this work, we introduce a deep learning-based protein pK _a predictor DeepKa, which is trained and validated with the pK _a values derived from continuous constant-pH molecular dynamics (CpHMD) simulations of 279 soluble proteins. Here, the CpHMD implemented in the Amber molecular dynamics package has been employed (Huang Y.J. Chem. Inf. Model.2018, 58, 1372-1383). Notably, to avoid discontinuities at the boundary, grid charges are proposed to represent protein electrostatics. We show that the prediction accuracy by DeepKa is close to that by CpHMD benchmarking simulations, validating DeepKa as an efficient protein pK _a predictor. In addition, the training and validation sets created in this study can be applied to the development of machine learning-based protein pK _a predictors in the future. Finally, the grid charge representation is general and applicable to other topics, such as the protein-ligand binding affinity prediction.