Selectivity Challenges in Docking Screens for GPCR Targets and Antitargets

DrugGym: A testbed for the economics of autonomous drug discovery

Drug discovery is stochastic. The effectiveness of candidate compounds in satisfying design objectives is unknown ahead of time, and the tools used for prioritization-predictive models and assays-are inaccurate and noisy. In a typical discovery campaign, thousands of compounds may be synthesized and tested before design objectives are achieved, with many others ideated but deprioritized. These challenges are well-documented, but assessing potential remedies has been difficult. We introduce DrugGym, a framework for modeling the stochastic process of drug discovery. Emulating biochemical assays with realistic surrogate models, we simulate the progression from weak hits to sub-micromolar leads with viable ADME. We use this testbed to examine how different ideation, scoring, and decision-making strategies impact statistical measures of utility, such as the probability of program success within predefined budgets and the expected costs to achieve target candidate profile (TCP) goals. We also assess the influence of affinity model inaccuracy, chemical creativity, batch size, and multi-step reasoning. Our findings suggest that reducing affinity model inaccuracy from 2 to 0.5 pIC50 units improves budget-constrained success rates tenfold. DrugGym represents a realistic testbed for machine learning methods applied to the hit-to-lead phase. Source code is available at www.drug-gym.org.

Discovery of Selective Dopamine Receptor Ligands Derived from (-)-Stepholidine via C-3 Alkoxylation and C-3/C-9 Dialkoxylation

We evaluated C-3 alkoxylated and C-3/C-9 dialkoxylated (-)-stepholidine analogues to probe the tolerance at the C-3 and C-9 positions of the tetrahydroprotoberberine (THPB) template toward affinity for dopamine receptors. A C-9 ethoxyl substituent appears optimal for D1R affinity since high D1R affinities were observed for compounds that contain an ethyl group at C-9, with larger C-9 substituents tending to decrease D1R affinity. A number of novel ligands were identified, such as compounds 12a and 12b, with nanomolar affinities for D1R and no affinity for either D2R or D3R, with compound 12a being identified as a D1R antagonist for both G-protein- and β-arrestin-based signaling. Compound 23b was identified as the most potent and selective D3R ligand containing a THPB template to date and functions as an antagonist for both G-protein- and β-arrestin-based signaling. Molecular docking and molecular dynamics studies validated the D1R and D3R affinity and selectivity of 12a, 12b, and 23b.

Structure-based virtual screening discovers potent and selective adenosine A1 receptor antagonists

Development of subtype-selective leads is essential in drug discovery campaigns targeting G protein-coupled receptors (GPCRs). Herein, a structure-based virtual screening approach to rationally design subtype-selective ligands was applied to the A₁ and A_2A adenosine receptors (A₁R and A_2AR). Crystal structures of these closely related subtypes revealed a non-conserved subpocket in the binding sites that could be exploited to identify A₁R selective ligands. A library of 4.6 million compounds was screened computationally against both receptors using molecular docking and 20 A₁R selective ligands were predicted. Of these, seven antagonized the A₁R with micromolar activities and several compounds displayed slight selectivity for this subtype. Twenty-seven analogs of two discovered scaffolds were designed, resulting in antagonists with nanomolar potency and up to 76-fold A₁R-selectivity. Our results show the potential of structure-based virtual screening to guide discovery and optimization of subtype-selective ligands, which could facilitate the development of safer drugs.

Enhanced compound-protein binding affinity prediction by representing protein multimodal information via a coevolutionary strategy

Due to the lack of a method to efficiently represent the multimodal information of a protein, including its structure and sequence information, predicting compound-protein binding affinity (CPA) still suffers from low accuracy when applying machine-learning methods. To overcome this limitation, in a novel end-to-end architecture (named FeatNN), we develop a coevolutionary strategy to jointly represent the structure and sequence features of proteins and ultimately optimize the mathematical models for predicting CPA. Furthermore, from the perspective of data-driven approach, we proposed a rational method that can utilize both high- and low-quality databases to optimize the accuracy and generalization ability of FeatNN in CPA prediction tasks. Notably, we visually interpret the feature interaction process between sequence and structure in the rationally designed architecture. As a result, FeatNN considerably outperforms the state-of-the-art (SOTA) baseline in virtual drug evaluation tasks, indicating the feasibility of this approach for practical use. FeatNN provides an outstanding method for higher CPA prediction accuracy and better generalization ability by efficiently representing multimodal information of proteins via a coevolutionary strategy.

Computational Chemistry for the Identification of Lead Compounds for Radiotracer Development

The use of computer-aided drug design (CADD) for the identification of lead compounds in radiotracer development is steadily increasing. Traditional CADD methods, such as structure-based and ligand-based virtual screening and optimization, have been successfully utilized in many drug discovery programs and are highlighted throughout this review. First, we discuss the use of virtual screening for hit identification at the beginning of drug discovery programs. This is followed by an analysis of how the hits derived from virtual screening can be filtered and culled to highly probable candidates to test in in vitro assays. We then illustrate how CADD can be used to optimize the potency of experimentally validated hit compounds from virtual screening for use in positron emission tomography (PET). Finally, we conclude with a survey of the newest techniques in CADD employing machine learning (ML).

In silico selectivity modeling of pyridine and pyrimidine based CYP11B1 and CYP11B2 inhibitors: A case study

DESIGN

of selective drug candidates for highly structural similar targets is a challenging task for researchers. The main objective of this study was to explore the selectivity modeling of pyridine and pyrimidine scaffold towards the highly homologous targets CYP11B1 and CYP11B2 enzymes by in silico (Molecular docking and QSAR) approaches. In this regard, a big dataset (n = 228) of CYP11B1 and CYP11B2 inhibitors were gathered and classified based on heterocyclic ring and the exhaustive analysis was carried out for pyridine and pyrimidinescaffolds. The LibDock algorithm was used to explore the binding pattern, screening, and identify the structural feature responsible for the selectivity of the ligands towards the studied targets. Finally, QSAR analysis was done to explore the correlation between various binding parameters and structural features responsible for the inhibitory activity and selectivity of the ligands in a quantitative way. The docking and QSAR analysis clearly revealed and distinguished the importance of structural features, functional groups attached for CYP11B2 and CYP11B1 selectivity for pyridine and pyrimidine analogs. Additionally, the docking analysis highlighted the differentiating amino acids residues for selectivity for ligands for each of the enzymes. The results obtained from this research work will be helpful in designing the selective CYP11B1/CYP11B2 inhibitors.

Docking-Based Virtual Screening Enables Prioritizing Protein Kinase Inhibitors With In Vitro Phenotypic Activity Against Schistosoma mansoni

Schistosomiasis is a parasitic neglected disease with praziquantel (PZQ) utilized as the main drug for treatment, despite its low effectiveness against early stages of the worm. To aid in the search for new drugs to tackle schistosomiasis, computer-aided drug design has been proved a helpful tool to enhance the search and initial identification of schistosomicidal compounds, allowing fast and cost-efficient progress in drug discovery. The combination of high-throughput in silico data followed by in vitro phenotypic screening assays allows the assessment of a vast library of compounds with the potential to inhibit a single or even several biological targets in a more time- and cost-saving manner. Here, we describe the molecular docking for in silico screening of predicted homology models of five protein kinases (JNK, p38, ERK1, ERK2, and FES) of Schistosoma mansoni against approximately 85,000 molecules from the Managed Chemical Compounds Collection (MCCC) of the University of Nottingham (UK). We selected 169 molecules predicted to bind to SmERK1, SmERK2, SmFES, SmJNK, and/or Smp38 for in vitro screening assays using schistosomula and adult worms. In total, 89 (52.6%) molecules were considered active in at least one of the assays. This approach shows a much higher efficiency when compared to using only traditional high-throughput in vitro screening assays, where initial positive hits are retrieved from testing thousands of molecules. Additionally, when we focused on compound promiscuity over selectivity, we were able to efficiently detect active compounds that are predicted to target all kinases at the same time. This approach reinforces the concept of polypharmacology aiming for “one drug-multiple targets”. Moreover, at least 17 active compounds presented satisfactory drug-like properties score when compared to PZQ, which allows for optimization before further in vivo screening assays. In conclusion, our data support the use of computer-aided drug design methodologies in conjunction with high-throughput screening approach.

Recent Advances in Structure, Function, and Pharmacology of Class A Lipid GPCRs: Opportunities and Challenges for Drug Discovery

Great progress has been made over the past decade in understanding the structural, functional, and pharmacological diversity of lipid GPCRs. From the first determination of the crystal structure of bovine rhodopsin in 2000, much progress has been made in the field of GPCR structural biology. The extraordinary progress in structural biology and pharmacology of GPCRs, coupled with rapid advances in computational approaches to study receptor dynamics and receptor-ligand interactions, has broadened our comprehension of the structural and functional facets of the receptor family members and has helped usher in a modern age of structure-based drug design and development. First, we provide a primer on lipid mediators and lipid GPCRs and their role in physiology and diseases as well as their value as drug targets. Second, we summarize the current advancements in the understanding of structural features of lipid GPCRs, such as the structural variation of their extracellular domains, diversity of their orthosteric and allosteric ligand binding sites, and molecular mechanisms of ligand binding. Third, we close by collating the emerging paradigms and opportunities in targeting lipid GPCRs, including a brief discussion on current strategies, challenges, and the future outlook.

Structure-Based Virtual Screening for Ligands of G Protein-Coupled Receptors: What Can Molecular Docking Do for You?

G protein-coupled receptors (GPCRs) constitute the largest family of membrane proteins in the human genome and are important therapeutic targets. During the last decade, the number of atomic-resolution structures of GPCRs has increased rapidly, providing insights into drug binding at the molecular level. These breakthroughs have created excitement regarding the potential of using structural information in ligand design and initiated a new era of rational drug discovery for GPCRs. The molecular docking method is now widely applied to model the three-dimensional structures of GPCR-ligand complexes and screen for chemical probes in large compound libraries. In this review article, we first summarize the current structural coverage of the GPCR superfamily and the understanding of receptor-ligand interactions at atomic resolution. We then present the general workflow of structure-based virtual screening and strategies to discover GPCR ligands in chemical libraries. We assess the state of the art of this research field by summarizing prospective applications of virtual screening based on experimental structures. Strategies to identify compounds with specific efficacy and selectivity profiles are discussed, illustrating the opportunities and limitations of the molecular docking method. Our overview shows that structure-based virtual screening can discover novel leads and will be essential in pursuing the next generation of GPCR drugs. SIGNIFICANCE STATEMENT: Extraordinary advances in the structural biology of G protein-coupled receptors have revealed the molecular details of ligand recognition by this large family of therapeutic targets, providing novel avenues for rational drug design. Structure-based docking is an efficient computational approach to identify novel chemical probes from large compound libraries, which has the potential to accelerate the development of drug candidates.

Nonantimicrobial Actions of Macrolides: Overview and Perspectives for Future Development

Macrolides are among the most widely prescribed broad spectrum antibacterials, particularly for respiratory infections. It is now recognized that these drugs, in particular azithromycin, also exert time-dependent immunomodulatory actions that contribute to their therapeutic benefit in both infectious and other chronic inflammatory diseases. Their increased chronic use in airway inflammation and, more recently, of azithromycin in COVID-19, however, has led to a rise in bacterial resistance. An additional crucial aspect of chronic airway inflammation, such as chronic obstructive pulmonary disease, as well as other inflammatory disorders, is the loss of epithelial barrier protection against pathogens and pollutants. In recent years, azithromycin has been shown with time to enhance the barrier properties of airway epithelial cells, an action that makes an important contribution to its therapeutic efficacy. In this article, we review the background and evidence for various immunomodulatory and time-dependent actions of macrolides on inflammatory processes and on the epithelium and highlight novel nonantibacterial macrolides that are being studied for immunomodulatory and barrier-strengthening properties to circumvent the risk of bacterial resistance that occurs with macrolide antibacterials. We also briefly review the clinical effects of macrolides in respiratory and other inflammatory diseases associated with epithelial injury and propose that the beneficial epithelial effects of nonantibacterial azithromycin derivatives in chronic inflammation, even given prophylactically, are likely to gain increasing attention in the future. SIGNIFICANCE STATEMENT: Based on its immunomodulatory properties and ability to enhance the protective role of the lung epithelium against pathogens, azithromycin has proven superior to other macrolides in treating chronic respiratory inflammation. A nonantibiotic azithromycin derivative is likely to offer prophylactic benefits against inflammation and epithelial damage of differing causes while preserving the use of macrolides as antibiotics.

Subtype-selective mechanisms of negative allosteric modulators binding to group I metabotropic glutamate receptors

Group I metabotropic glutamate receptors (mGlu₁ and mGlu₅) are promising targets for multiple psychiatric and neurodegenerative disorders. Understanding the subtype selectivity of mGlu₁ and mGlu₅ allosteric sites is essential for the rational design of novel modulators with single- or dual-target mechanism of action. In this study, starting from the deposited mGlu₁ and mGlu₅ crystal structures, we utilized computational modeling approaches integrating docking, molecular dynamics simulation, and efficient post-trajectory analysis to reveal the subtype-selective mechanism of mGlu₁ and mGlu₅ to 10 diverse drug scaffolds representing known negative allosteric modulators (NAMs) in the literature. The results of modeling identified six pairs of non-conserved residues and four pairs of conserved ones as critical features to distinguish the selective NAMs binding to the corresponding receptors. In addition, nine pairs of residues are beneficial to the development of novel dual-target NAMs of group I metabotropic glutamate receptors. Furthermore, the binding modes of a reported dual-target NAM (VU0467558) in mGlu₁ and mGlu₅ were predicted to verify the identified residues that play key roles in the receptor selectivity and the dual-target binding. The results of this study can guide rational structure-based design of novel NAMs, and the approach can be generally applicable to characterize the features of selectivity for other G-protein-coupled receptors.

Preliminary modelling as the first stage of targeted organic synthesis

The review aims to present a classification and applicability analysis of methods for preliminary molecular modelling for targeted organic, catalytic and biocatalytic synthesis. The following three main approaches are considered as a primary classification of the methods: modelling of the target – ligand coordination without structural information on both the target and the resulting complex; calculations based on experimentally obtained structural information about the target; and dynamic simulation of the target – ligand complex and the reaction mechanism with calculation of the free energy of the reaction. The review is meant for synthetic chemists to be used as a guide for building an algorithm for preliminary modelling and synthesis of structures with specified properties.

The bibliography includes 353 references.

Structure-based in silico approaches for drug discovery against Mycobacterium tuberculosis

Mycobacterium tuberculosis is the causative agent of TB and was estimated to cause 1.4 million death in 2019, alongside 10 million new infections. Drug resistance is a growing issue, with multi-drug resistant infections representing 3.3% of all new infections, hence novel antimycobacterial drugs are urgently required to combat this growing health emergency. Alongside this, increased knowledge of gene essentiality in the pathogenic organism and larger compound databases can aid in the discovery of new drug compounds. The number of protein structures, X-ray based and modelled, is increasing and now accounts for greater than > 80% of all predicted M. tuberculosis proteins; allowing novel targets to be investigated. This review will focus on structure-based in silico approaches for drug discovery, covering a range of complexities and computational demands, with associated antimycobacterial examples. This includes molecular docking, molecular dynamic simulations, ensemble docking and free energy calculations. Applications of machine learning onto each of these approaches will be discussed. The need for experimental validation of computational hits is an essential component, which is unfortunately missing from many current studies. The future outlooks of these approaches will also be discussed.

Structural Insights Accelerate the Discovery of Opioid Alternatives

Opioids such as morphine and oxycodone are analgesics frequently prescribed for the treatment of moderate or severe pain. Unfortunately, these medications are associated with exceptionally high abuse potentials and often cause fatal side effects, mainly through the μ-opioid receptor (MOR). Efforts to discover novel, safer, and more efficacious analgesics targeting MOR have encountered challenges. In this review, we summarize alternative strategies and targets that could be used to develop safer nonopioid analgesics. A molecular understanding of G protein-coupled receptor activation and signaling has illuminated not only the complexities of receptor pharmacology but also the potential for pathway-selective agonists and allosteric modulators as safer medications. The availability of structures of pain-related receptors, in combination with high-throughput computational tools, has accelerated the discovery of multitarget ligands with promising pharmacological profiles. Emerging clinical evidence also supports the notion that drugs targeting peripheral opioid receptors have potential as improved analgesic agents.

Ligand- and Structure-Based Analysis of Deep Learning-Generated Potential α2a Adrenoceptor Agonists

The α2a adrenoceptor is a medically relevant subtype of the G protein-coupled receptor family. Unfortunately, high-throughput techniques aimed at producing novel drug leads for this receptor have been largely unsuccessful because of the complex pharmacology of adrenergic receptors. As such, cutting-edge in silico ligand- and structure-based assessment and de novo deep learning methods are well positioned to provide new insights into protein-ligand interactions and potential active compounds. In this work, we (i) collect a dataset of α2a adrenoceptor agonists and provide it as a resource for the drug design community; (ii) use the dataset as a basis to generate candidate-active structures via deep learning; and (iii) apply computational ligand- and structure-based analysis techniques to gain new insights into α2a adrenoceptor agonists and assess the quality of the computer-generated compounds. We further describe how such assessment techniques can be applied to putative chemical probes with a case study involving proposed medetomidine-based probes.

An open source chemical structure curation pipeline using RDKit

BACKGROUND

The ChEMBL database is one of a number of public databases that contain bioactivity data on small molecule compounds curated from diverse sources. Incoming compounds are typically not standardised according to consistent rules. In order to maintain the quality of the final database and to easily compare and integrate data on the same compound from different sources it is necessary for the chemical structures in the database to be appropriately standardised.

RESULTS

A chemical curation pipeline has been developed using the open source toolkit RDKit. It comprises three components: a Checker to test the validity of chemical structures and flag any serious errors; a Standardizer which formats compounds according to defined rules and conventions and a GetParent component that removes any salts and solvents from the compound to create its parent. This pipeline has been applied to the latest version of the ChEMBL database as well as uncurated datasets from other sources to test the robustness of the process and to identify common issues in database molecular structures.

CONCLUSION

All the components of the structure pipeline have been made freely available for other researchers to use and adapt for their own use. The code is available in a GitHub repository and it can also be accessed via the ChEMBL Beaker webservices. It has been used successfully to standardise the nearly 2 million compounds in the ChEMBL database and the compound validity checker has been used to identify compounds with the most serious issues so that they can be prioritised for manual curation.

The Advent of Generative Chemistry

Generative adversarial networks (GANs), first published in 2014, are among the most important concepts in modern artificial intelligence (AI). Bridging deep learning and game theory, GANs are used to generate or "imagine" new objects with desired properties. Since 2016, multiple GANs with reinforcement learning (RL) have been successfully applied in pharmacology for de novo molecular design. Those techniques aim at a more efficient use of the data and a better exploration of the chemical space. We review recent advances for the generation of novel molecules with desired properties with a focus on the applications of GANs, RL, and related techniques. We also discuss the current limitations and challenges in the new growing field of generative chemistry.

Nicotine Bound to Its Receptors: New Structures for a Vexing Pathopharmacological Problem

In this issue of Neuron, Gharpure et al. (2019) nearly complete atomic-level descriptions for binding, permeation, and subunit interactions at the two major heteropentameric nicotinic receptors-those governing nicotine's rewarding and aversive effects. Can we finally design highly selective and useful nicotinic drugs?

Structure-based discovery of potent and selective melatonin receptor agonists

Melatonin receptors MT1 and MT2 are involved in synchronizing circadian rhythms and are important targets for treating sleep and mood disorders, type-2 diabetes and cancer. Here, we performed large scale structure-based virtual screening for new ligand chemotypes using recently solved high-resolution 3D crystal structures of agonist-bound MT receptors. Experimental testing of 62 screening candidates yielded the discovery of 10 new agonist chemotypes with sub-micromolar potency at MT receptors, with compound 21 reaching EC50 of 0.36 nM. Six of these molecules displayed selectivity for MT2 over MT1. Moreover, two most potent agonists, including 21 and a close derivative of melatonin, 28, had dramatically reduced arrestin recruitment at MT2, while compound 37 was devoid of Gi signaling at MT1, implying biased signaling. This study validates the suitability of the agonist-bound orthosteric pocket in the MT receptor structures for the structure-based discovery of selective agonists.

Performance of virtual screening against GPCR homology models: Impact of template selection and treatment of binding site plasticity

Rational drug design for G protein-coupled receptors (GPCRs) is limited by the small number of available atomic resolution structures. We assessed the use of homology modeling to predict the structures of two therapeutically relevant GPCRs and strategies to improve the performance of virtual screening against modeled binding sites. Homology models of the D2 dopamine (D2R) and serotonin 5-HT2A receptors (5-HT2AR) were generated based on crystal structures of 16 different GPCRs. Comparison of the homology models to D2R and 5-HT2AR crystal structures showed that accurate predictions could be obtained, but not necessarily using the most closely related template. Assessment of virtual screening performance was based on molecular docking of ligands and decoys. The results demonstrated that several templates and multiple models based on each of these must be evaluated to identify the optimal binding site structure. Models based on aminergic GPCRs showed substantial ligand enrichment and there was a trend toward improved virtual screening performance with increasing binding site accuracy. The best models even yielded ligand enrichment comparable to or better than that of the D2R and 5-HT2AR crystal structures. Methods to consider binding site plasticity were explored to further improve predictions. Molecular docking to ensembles of structures did not outperform the best individual binding site models, but could increase the diversity of hits from virtual screens and be advantageous for GPCR targets with few known ligands. Molecular dynamics refinement resulted in moderate improvements of structural accuracy and the virtual screening performance of snapshots was either comparable to or worse than that of the raw homology models. These results provide guidelines for successful application of structure-based ligand discovery using GPCR homology models.

Residues and residue pairs of evolutionary importance differentially direct signaling bias of D2 dopamine receptors

The D2 dopamine receptor and the serotonin 5-hydroxytryptamine 2A receptor (5-HT2A) are closely-related G-protein-coupled receptors (GPCRs) from the class A bioamine subfamily. Despite structural similarity, they respond to distinct ligands through distinct downstream pathways, whose dysregulation is linked to depression, bipolar disorder, addiction, and psychosis. They are important drug targets, and it is important to understand how their bias toward G-protein versus β-arrestin signaling pathways is regulated. Previously, evolution-based computational approaches, difference Evolutionary Trace and Evolutionary Trace-Mutual information (ET-Mip), revealed residues and residue pairs that, when switched in the D2 receptor to the corresponding residues from 5-HT2A, altered ligand potency and G-protein activation efficiency. We have tested these residue swaps for their ability to trigger recruitment of β-arrestin2 in response to dopamine or serotonin. The results reveal that the selected residues modulate agonist potency, maximal efficacy, and constitutive activity of β-arrestin2 recruitment. Whereas dopamine potency for most variants was similar to that for WT and lower than for G-protein activation, potency in β-arrestin2 recruitment for N124H^3.42 was more than 5-fold higher. T205M^5.54 displayed high constitutive activity, enhanced dopamine potency, and enhanced efficacy in β-arrestin2 recruitment relative to WT, and L379F^6.41 was virtually inactive. These striking differences from WT activity were largely reversed by a compensating mutation (T205M^5.54/L379F^6.41) at residues previously identified by ET-Mip as functionally coupled. The observation that the signs and relative magnitudes of the effects of mutations in several cases are at odds with their effects on G-protein activation suggests that they also modulate signaling bias.

Molecular pharmacology of metabotropic receptors targeted by neuropsychiatric drugs

Metabotropic receptors are responsible for so-called 'slow synaptic transmission' and mediate the effects of hundreds of peptide and non-peptide neurotransmitters and neuromodulators. Over the past decade or so, a revolution in membrane-protein structural determination has clarified the molecular determinants responsible for the actions of these receptors. This Review focuses on the G protein-coupled receptors (GPCRs) that are targets of neuropsychiatric drugs and shows how insights into the structure and function of these important synaptic proteins are accelerating understanding of their actions. Notably, elucidating the structure and function of GPCRs should enhance the structure-guided discovery of novel chemical tools with which to manipulate and understand these synaptic proteins.

Electrostatic potential and non-covalent interactions analysis for the design of selective 5-HT7 ligands

Despite the significant improvement of methodology in the field of the in silico drug discovery, the search for selective drugs is still far from trivial. This is especially relevant in the case of designing new medicaments for treatment of central nervous system disorders. In this work, we present a new approach based on the molecular docking and the following electronic properties analysis of ligands' binding poses (electrostatic potential distribution analysis and quantitative topological analysis of the electron density distribution). The proposed protocol significantly increases the success rate of the selective 5-HT₇R ligands against 5-HT_1AR selection from the prepared databases (the rise from 33.3% to 77.8% and from 22.7% to 62.5% for training and testing sets, respectively). The presented approach can be applied as a supportive method in the virtual screening of ligands databases.

Bifunctional opioid receptor ligands as novel analgesics

Prolonged treatment of chronic severe pain with opioid analgesics is frought with problematic adverse effects including tolerance, dependence, and life-threatening respiratory depression. Though these effects are mediated predominately through preferential activation of μ opioid peptide (μOP) receptors, there is an emerging appreciation that actions at κOP and δOP receptors contribute to the observed pharmacologic and behavioral profile of μOP receptor agonists and may be targeted simultaneously to afford improved analgesic effects. Recent developments have also identified the related nociceptin opioid peptide (NOP) receptor as a key modulator of the effects of μOP receptor signaling. We review here the available literature describing OP neurotransmitter systems and highlight recent drug and probe design strategies.