Comparison of Orexin 1 and Orexin 2 Ligand Binding Modes Using X-ray Crystallography and Computational Analysis

Structure-guided discovery of orexin receptor-binding PET ligands

Molecular imaging using positron emission tomography (PET) can serve as a promising tool for visualizing biological targets in the brain. Insights into the expression pattern and the in vivo imaging of the G protein-coupled orexin receptors OX1R and OX2R will further our understanding of the orexin system and its role in various physiological and pathophysiological processes. Guided by crystal structures of our lead compound JH112 and the approved hypnotic drug suvorexant bound to OX1R and OX2R, respectively, we herein describe the design and synthesis of two novel radioligands, [18F]KD23 and [18F]KD10. Key to the success of our structural modifications was a bioisosteric replacement of the triazole moiety with a fluorophenyl group. The 19F-substituted analog KD23 showed high affinity for the OX1R and selectivity over OX2R, while the high affinity ligand KD10 displayed similar Ki values for both subtypes. Radiolabeling starting from the respective pinacol ester precursors resulted in excellent radiochemical yields of 93% and 88% for [18F]KD23 and [18F]KD10, respectively, within 20 min. The new compounds will be useful in PET studies aimed at subtype-selective imaging of orexin receptors in brain tissue.

Distinct binding hotspots for natural and synthetic agonists of FFA4 from in silico approaches

FFA4 has gained interest in recent years since its deorphanization in 2005 and the characterization of the Free Fatty Acids receptors family for their therapeutic potential in metabolic disorders. The expression of FFA4 (also known as GPR120) in numerous organs throughout the human body makes this receptor a highly potent target, particularly in fat sensing and diet preference. This offers an attractive approach to tackle obesity and related metabolic diseases. Recent cryo-EM structures of the receptor have provided valuable information for a potential active state although the previous studies of FFA4 presented diverging information. We performed molecular docking and molecular dynamics simulations of four agonist ligands, TUG-891, Linoleic acid, α-Linolenic acid, and Oleic acid, based on a homology model. Our simulations, which accumulated a total of 2 μs of simulation, highlighted two binding hotspots at Arg992.64 and Lys293 (ECL3). The results indicate that the residues are located in separate areas of the binding pocket and interact with various types of ligands, implying different potential active states of FFA4 and a highly adaptable binding intra-receptor pocket. This article proposes additional structural characteristics and mechanisms for agonist binding that complement the experimental structures.

G Protein-Coupled Receptor-Ligand Pose and Functional Class Prediction

G protein-coupled receptor (GPCR) transmembrane protein family members play essential roles in physiology. Numerous pharmaceuticals target GPCRs, and many drug discovery programs utilize virtual screening (VS) against GPCR targets. Improvements in the accuracy of predicting new molecules that bind to and either activate or inhibit GPCR function would accelerate such drug discovery programs. This work addresses two significant research questions. First, do ligand interaction fingerprints provide a substantial advantage over automated methods of binding site selection for classical docking? Second, can the functional status of prospective screening candidates be predicted from ligand interaction fingerprints using a random forest classifier? Ligand interaction fingerprints were found to offer modest advantages in sampling accurate poses, but no substantial advantage in the final set of top-ranked poses after scoring, and, thus, were not used in the generation of the ligand-receptor complexes used to train and test the random forest classifier. A binary classifier which treated agonists, antagonists, and inverse agonists as active and all other ligands as inactive proved highly effective in ligand function prediction in an external test set of GPR31 and TAAR2 candidate ligands with a hit rate of 82.6% actual actives within the set of predicted actives.

Machine learning models to predict ligand binding affinity for the orexin 1 receptor

The orexin 1 receptor (OX1R) is a G-protein coupled receptor that regulates a variety of physiological processes through interactions with the neuropeptides orexin A and B. Selective OX1R antagonists exhibit therapeutic effects in preclinical models of several behavioral disorders, including drug seeking and overeating. However, currently there are no selective OX1R antagonists approved for clinical use, fueling demand for novel compounds that act at this target. In this study, we meticulously curated a dataset comprising over 1300 OX1R ligands using a stringent filter and criteria cascade. Subsequently, we developed highly predictive quantitative structure-activity relationship (QSAR) models employing the optimized hyper-parameters for the random forest machine learning algorithm and twelve 2D molecular descriptors selected by recursive feature elimination with a 5-fold cross-validation process. The predictive capacity of the QSAR model was further assessed using an external test set and enrichment study, confirming its high predictivity. The practical applicability of our final QSAR model was demonstrated through virtual screening of the DrugBank database. This revealed two FDA-approved drugs (isavuconazole and cabozantinib) as potential OX1R ligands, confirmed by radiolabeled OX1R binding assays. To our best knowledge, this study represents the first report of highly predictive QSAR models on a large comprehensive dataset of diverse OX1R ligands, which should prove useful for the discovery and design of new compounds targeting this receptor.

Interaction modes of human orexin 2 receptor with selective and nonselective antagonists studied by NMR spectroscopy

Orexin neuropeptides have many physiological roles in the sleep-wake cycle, feeding behavior, reward demands, and stress responses by activating cognitive receptors, the orexin receptors (OX1R and OX2R), distributed in the brain. There are only subtle differences between OX1R and OX2R in the orthosteric site, which has hindered the rational development of subtype-selective antagonists. In this study, we utilized solution-state NMR to capture the structural plasticity of OX2R labeled with 13CH3-ε-methionine in complex with antagonists. Mutations in the orthosteric site allosterically affected the intracellular tip of TM6. Ligand exchange experiments with the subtype-selective EMPA and the nonselective suvorexant identified three methionine residues that were substantially perturbed. The NMR spectra suggested that the suvorexant-bound state exhibited more structural plasticity than the EMPA-bound state, which has not been foreseen from the close similarity of their crystal structures, providing insights into dynamic features to be considered in understanding the ligand recognition mode.

Discovery and first-time disclosure of CVN766, an exquisitely selective orexin 1 receptor antagonist

Modulators of orexin receptors are being developed for neurological illnesses such as sleep disorders, addictive behaviours and other psychiatric diseases. We herein describe the discovery of CVN766, a potent orexin 1 receptor antagonist that has greater than 1000-fold selectivity for the orexin 1 receptor over the orexin 2 receptor and demonstrates low off target hits in a diversity screen. In agreement with its in vitro ADME data, CVN766 demonstrated moderate in vivo clearance in rodents and displayed good brain permeability and target occupancy. This drug candidate is currently being investigated in clinical trials for schizophrenia and related psychiatric conditions.

New drug approvals for 2022: Synthesis and clinical applications

The U.S. Food and Drug Administration has approved a total of 37 new drugs in 2022, which are composed of 20 chemical entities and 17 biologics. In particular, 20 chemical entities, including 17 small molecule drugs, 1 radiotherapy, and 2 diagnostic agents, provide privileged scaffolds, breakthrough clinical benefits, and a new mechanism of action for the discovery of more potent clinical candidates. The structure-based drug development with clear targets and fragment-based drug development with privileged scaffolds have always been the important modules in the field of drug discovery, which could easily bypass the patent protection and bring about improved biological activity. Therefore, we summarized the relevant valuable information about clinical application, mechanism of action, and chemical synthesis of 17 newly approved small molecule drugs in 2022. We hope this timely and comprehensive review could bring about creative and elegant inspiration on the synthetic methodologies and mechanism of action for the discovery of new drugs with novel chemical scaffolds and extended clinical indications.

Targeting orexin receptors: Recent advances in the development of subtype selective or dual ligands for the treatment of neuropsychiatric disorders

Orexin-A and orexin-B, also named hypocretin-1 and hypocretin-2, are two hypothalamic neuropeptides highly conserved across mammalian species. Their effects are mediated by two distinct G protein-coupled receptors, namely orexin receptor type 1 (OX1-R) and type 2 (OX2-R), which share 64% amino acid identity. Given the wide expression of OX-Rs in different central nervous system and peripheral areas and the several pathophysiological functions in which they are involved, including sleep-wake cycle regulation (mainly mediated by OX2-R), emotion, panic-like behaviors, anxiety/stress, food intake, and energy homeostasis (mainly mediated by OX1-R), both subtypes represent targets of interest for many structure-activity relationship (SAR) campaigns carried out by pharmaceutical companies and academies. However, before 2017 the research was predominantly directed towards dual-orexin ligands, and limited chemotypes were investigated. Analytical characterizations, including resolved structures for both OX1-R and OX2-R in complex with agonists and antagonists, have improved the understanding of the molecular basis of receptor recognition and are assets for medicinal chemists in the design of subtype-selective ligands. This review is focused on the medicinal chemistry aspects of small molecules acting as dual or subtype selective OX1-R/OX2-R agonists and antagonists belonging to different chemotypes and developed in the last years, including radiolabeled OX-R ligands for molecular imaging. Moreover, the pharmacological effects of the most studied ligands in different neuropsychiatric diseases, such as sleep, mood, substance use, and eating disorders, as well as pain, have been discussed. Poly-pharmacology applications and multitarget ligands have also been considered.

The Role of Allylic Strain for Conformational Control in Medicinal Chemistry

It is axiomatic in medicinal chemistry that optimization of the potency of a small molecule at a macromolecular target requires complementarity between the ligand and target. In order to minimize the conformational penalty on binding, both enthalpically and entropically, it is therefore preferred to have the ligand preorganized in the bound conformation. In this Perspective, we highlight the role of allylic strain in controlling conformational preferences. Allylic strain was originally described for carbon-based allylic systems, but the same principles apply to other types of structure with sp² or pseudo-sp² arrangements. These systems include benzylic (including heteroaryl methyl) positions, amides, N-aryl groups, aryl ethers, and nucleotides. We have derived torsion profiles from small molecule X-ray structures for these systems. Through multiple examples, we show how these effects have been applied in drug discovery and how they can be used prospectively to influence conformation in the design process.

The abuse potential of lemborexant, a dual orexin receptor antagonist, according to the 8 factors of the Controlled Substances Act

RATIONALE

Lemborexant (LEM) is a dual orexin receptor antagonist (DORA) approved in multiple countries including the USA, Japan, Canada, Australia, and several Asian countries for the treatment of insomnia in adults. As a compound with central nervous system activity, it is important to understand the abuse potential of LEM with respect to public health.

OBJECTIVES

This review discusses data for LEM relevant to each of the 8 factors of the United States Controlled Substances Act.

RESULTS

LEM did not demonstrate abuse potential in nonclinical testing and was associated with a low incidence of abuse-related adverse events in clinical study participants with insomnia disorder. Similar to other DORAs that have been evaluated (eg., almorexant, suvorexant (SUV), and daridorexant), LEM and the positive controls (zolpidem and SUV) also showed drug liking in a phase 1 abuse potential study that enrolled subjects who used sedatives recreationally. However, internet surveillance of SUV and the FDA Adverse Events Reporting System suggests that drugs in the DORA class display very low abuse-related risks in the community. Additionally, as described in FDA-approved labeling, it does not carry physical dependence and withdrawal risks.

CONCLUSIONS

LEM, similar to most other prescription insomnia medications, was placed into Schedule IV. However, LEM and other drugs in the DORA class may have a lower potential for abuse as suggested by real-world postmarketing data from federal surveys and internet surveillance, and thus may have lower risks to public health than Schedule IV benzodiazepines and nonbenzodiazepine hypnotics that potentiate GABA signaling.

Design and synthesis of novel orexin 2 receptor agonists with a 1,3,5‑trioxazatriquinane skeleton

A novel series of 1,3,5‑trioxazatriquinane with multiple effective residues (TriMER) derivatives with amino-methylene side chains was designed and synthesized based on the docking-simulation results between orexin receptors (OXRs) and TriMER-type OXR antagonists. In vitro screening against orexin receptors identified six TriMER derivatives with a cis side-chain configuration, and, among these, 20d and 28d showed full agonist activity against OX₂R at a concentration of 10 µM. To determine the absolute stereochemistry of these hit compounds, we also conducted the first asymmetric synthesis of a 1,3,5‑trioxazatriquinane skeleton using a Katsuki-Sharpless asymmetric epoxidation as the key reaction and obtained a set of the individual stereoisomers. After evaluating their activity, (+)-20d (EC₅₀ = 3.87 μM for OX₂R) and (+)-28d (EC₅₀ = 1.62 μM for OX₂R) were determined as eutomers for OX₂R agonist activity. Our results provide a new class of skeleton consisting of an (R)-1,3,5‑trioxazatriquinane core with flexible methylene linkers and hydrophobic substituents at the terminals of the side chains via carbamates/sulfonamides as OX₂R agonists.

Protein Design Strategies for the Structural–Functional Studies of G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are an important family of membrane proteins responsible for many physiological functions in human body. High resolution GPCR structures are required to understand their molecular mechanisms and perform rational drug design, as GPCRs play a crucial role in a variety of diseases. That is difficult to obtain for the wild-type proteins because of their low stability. In this review, we discuss how this problem can be solved by using protein design strategies developed to obtain homogeneous stabilized GPCR samples for crystallization and cryoelectron microscopy.

Antagonization of OX1 Receptor Potentiates CB2 Receptor Function in Microglia from APPSw/Ind Mice Model

Microdialysis assays demonstrated a possible role of orexin in the regulation of amyloid beta peptide (Aß) levels in the hippocampal interstitial fluid in the APP transgenic model. CB₂R is overexpressed in activated microglia, showing a neuroprotective effect. These two receptors may interact, forming CB₂-OX₁-Hets and becoming a new target to combat Alzheimer's disease. Aims: Demonstrate the potential role of CB₂-OX₁-Hets expression and function in microglia from animal models of Alzheimer's disease. Receptor heteromer expression was detected by immunocytochemistry, bioluminescence resonance energy transfer (BRET) and proximity ligation assay (PLA) in transfected HEK-293T cells and microglia primary cultures. Quantitation of signal transduction events in a heterologous system and in microglia cells was performed using the AlphaScreen^® SureFire^® kit, western blot, the GCaMP6 calcium sensor and the Lance Ultra cAMP kit (PerkinElmer). The formation of CB₂-OX₁ receptor complexes in transfected HEK-293T cells has been demonstrated. The tetrameric complex is constituted by one CB₂R homodimer, one OX₁R homodimer and two G proteins, a G_i and a G_q. The use of TAT interfering peptides showed that the CB₂-OX₁ receptor complex interface is TM4-TM5. At the functional level it has been observed that the OX₁R antagonist, SB334867, potentiates the action induced by CB₂R agonist JWH133. This effect is observed in transfected HEK-293T cells and microglia, and it is stronger in the Alzheimer's disease (AD) animal model APPSw/Ind where the expression of the complex assessed by the proximity ligation assay indicates an increase in the number of complexes compared to resting microglia. The CB₂-OX₁ receptor complex is overexpressed in microglia from AD animal models where OX₁R antagonists potentiate the neuroprotective actions of CB₂R activation. Taken together, these results point to OX₁R antagonists as drugs with therapeutic potential to combat AD. Data access statement: Raw data will be provided by the corresponding author upon reasonable requirement.

Augmented Hill-Climb increases reinforcement learning efficiency for language-based de novo molecule generation

A plethora of AI-based techniques now exists to conduct de novo molecule generation that can devise molecules conditioned towards a particular endpoint in the context of drug design. One popular approach is using reinforcement learning to update a recurrent neural network or language-based de novo molecule generator. However, reinforcement learning can be inefficient, sometimes requiring up to 10⁵ molecules to be sampled to optimize more complex objectives, which poses a limitation when using computationally expensive scoring functions like docking or computer-aided synthesis planning models. In this work, we propose a reinforcement learning strategy called Augmented Hill-Climb based on a simple, hypothesis-driven hybrid between REINVENT and Hill-Climb that improves sample-efficiency by addressing the limitations of both currently used strategies. We compare its ability to optimize several docking tasks with REINVENT and benchmark this strategy against other commonly used reinforcement learning strategies including REINFORCE, REINVENT (version 1 and 2), Hill-Climb and best agent reminder. We find that optimization ability is improved ~ 1.5-fold and sample-efficiency is improved ~ 45-fold compared to REINVENT while still delivering appealing chemistry as output. Diversity filters were used, and their parameters were tuned to overcome observed failure modes that take advantage of certain diversity filter configurations. We find that Augmented Hill-Climb outperforms the other reinforcement learning strategies used on six tasks, especially in the early stages of training or for more difficult objectives. Lastly, we show improved performance not only on recurrent neural networks but also on a reinforcement learning stabilized transformer architecture. Overall, we show that Augmented Hill-Climb improves sample-efficiency for language-based de novo molecule generation conditioning via reinforcement learning, compared to the current state-of-the-art. This makes more computationally expensive scoring functions, such as docking, more accessible on a relevant timescale.

Interaction between the Human OX2 Orexin Receptor and Suvorexant and Some of Its Analogues: SAPT (DFT) Interaction Energy Decomposition Analysis

In this study, the interaction energy (E_int) of suvorexant (as an orexin receptor antagonist) and some of its analogues with the important residues of the human OX2 orexin receptor, determined by molecular docking, is calculated using the symmetry-adapted perturbation theory-density functional theory (SAPT (DFT)) method. Also, the important residues with the dominant interaction with each ligand are determined based on the obtained SAPT (DFT) interaction energies. To analyze the interaction of the receptor with each ligand, the decomposition of E_int to its constituent components including electrostatic (E_ele), exchange (E_ex), induction (E_ind), and exchange-induction (E_ex-ind), dispersion (E_disp), and exchange-dispersion (E_ex-disp) is performed. The change of interaction energy components with the replacement of the benzoxazole part of suvorexant by pyrimidine containing different functional groups, thieno pyrimidine, and furo pyrimidine is also investigated, separately. It is found that the change in E_int, due to these replacements, is controlled more by the variation of the electrostatic interaction energy component of E_int than by the other interaction energy components. A linear correlation (R² = 0.91) is found for the variation of E_int versus experimental ligand-binding affinities. Also, the existence of the linear correlation for the variation of the interaction energy components with experimental ligand-binding affinities is investigated. The variation of the electrostatic component versus experimental ligand-binding affinities shows a more linear correlation compared to the other interaction energy components.

Applications of Isosteres of Piperazine in the Design of Biologically Active Compounds: Part 2

Applications of piperazine and homopiperazine in drug design are well-established, and these heterocycles have found use as both scaffolding and terminal elements and also as a means of introducing a water-solubilizing element into a molecule. In the accompanying review (10.1021/acs.jafc.2c00726), we summarized applications of piperazine and homopiperazine and their fused ring homologues in bioactive compound design along with illustrations of the use of 4-substituted piperidines and a sulfoximine-based mimetic. In this review, we discuss applications of pyrrolidine- and fused-pyrrolidine-based mimetics of piperazine and homopiperazine and illustrate derivatives of azetidine that include stretched and spirocyclic motifs, along with applications of a series of diaminocycloalkanes.

Applications of Isosteres of Piperazine in the Design of Biologically Active Compounds: Part 1

Piperazine and homopiperazine are well-studied heterocycles in drug design that have found gainful application as scaffolds and terminal elements and for enhancing the aqueous solubility of a molecule. The optimization of drug candidates that incorporate these heterocycles in an effort to refine potency, selectivity, and developability properties has stimulated the design and evaluation of a wide range of bioisosteres that can offer advantage. In this review, we summarize the design and application of bioisosteres of piperazine and homopiperazine that have almost exclusively been in the drug design arena. While there are ∼100 approved drugs that incorporate a piperazine ring, only a single marketed agricultural product is built on this heterocycle. As part of the review, we discuss some of the potential reasons underlying the relatively low level of importance of this heterocycle to the design of agrochemicals and highlight the potential opportunities for their use in contemporary research programs.

Daridorexant for the treatment of insomnia disorder: findings and implications

PURPOSE

The involvement of the orexin system in the physiopathology of insomnia has been rapidly increasing in understanding. In this sense, daridorexant was the third orexin receptor antagonist approved by the FDA in January 2022. This review aims to summarize the chemistry, pharmacodynamics, pharmacokinetics, efficacy, safety, and tolerability profile of daridorexant for the treatment of insomnia disorder.

METHODS

We performed a review of daridorexant for the treatment of insomnia disorder. The search was carried out in Medline via PubMed, Embase, and clinical trials, up to March 2022.

RESULTS

Daridorexant 25 and 50 mg had more significant improvement for the wake after sleep onset (WASO), latency to persistent sleep (LPS), and subjective total sleep time (sTST) than placebo. In addition, daridorexant 50 mg was better for Insomnia Daytime Symptoms and Impacts Questionnaire (IDSIQ) than placebo. The most common adverse events were nasopharyngitis and headache.

CONCLUSION

Daridorexant was efficacious and safe. Studies that evaluate the long-term safety and compare daridorexant with benzodiazepines, benzodiazepine receptor agonists, sedative antidepressants, and other orexin receptor antagonists are required.

Molecular mechanism of the wake-promoting agent TAK-925

The OX2 orexin receptor (OX2R) is a highly expressed G protein-coupled receptor (GPCR) in the brain that regulates wakefulness and circadian rhythms in humans. Antagonism of OX2R is a proven therapeutic strategy for insomnia drugs, and agonism of OX2R is a potentially powerful approach for narcolepsy type 1, which is characterized by the death of orexinergic neurons. Until recently, agonism of OX2R had been considered 'undruggable.' We harness cryo-electron microscopy of OX2R-G protein complexes to determine how the first clinically tested OX2R agonist TAK-925 can activate OX2R in a highly selective manner. Two structures of TAK-925-bound OX2R with either a Gq mimetic or Gi reveal that TAK-925 binds at the same site occupied by antagonists, yet interacts with the transmembrane helices to trigger activating microswitches. Our structural and mutagenesis data show that TAK-925's selectivity is mediated by subtle differences between OX1 and OX2 receptor subtypes at the orthosteric pocket. Finally, differences in the polarity of interactions at the G protein binding interfaces help to rationalize OX2R's coupling selectivity for Gq signaling. The mechanisms of TAK-925's binding, activation, and selectivity presented herein will aid in understanding the efficacy of small molecule OX2R agonists for narcolepsy and other circadian disorders.

Investigation of Tumor Cells and Receptor-Ligand Simulation Models for the Development of PET Imaging Probes Targeting PSMA and GRPR and a Possible Crosstalk between the Two Receptors

Prostate-specific membrane antigen (PSMA) and gastrin-releasing peptide receptor (GRPR) have both been used in nuclear medicine as targets for molecular imaging and therapy of prostate (PCa) and breast cancer (BCa). Three bioconjugate probes, the PSMA specific: [⁶⁸Ga]Ga-1, ((HBED-CC)-Ahx-Lys-NH-CO-NH Glu or PSMA-11), the GRPR specific: [⁶⁸Ga]Ga-2, ((HBED-CC)-4-amino-1-carboxymethyl piperidine-[D-Phe⁶, Sta¹³]BN(6-14), a bombesin (BN) analogue), and 3 (the BN analogue: 4-amino-1-carboxymethyl piperidine-[(R)-Phe⁶, Sta¹³]BN(6-14) connected with the fluorescent dye, BDP-FL), were synthesized and tested in vitro with PCa and BCa cell lines, more specifically, with PCa cells, PC-3 and LNCaP, with BCa cells, T47D, MDA-MB-231, and with the in-house created PSMA-overexpressing PC-3^(PSMA), T47D^(PSMA), and MDA-MB-231^(PSMA). In addition, biomolecular simulations were conducted on the association of 1 and 2 with PSMA and GRPR. The PSMA overexpression resulted in an increase of cell-bound radioligand [⁶⁸Ga]Ga-1 (PSMA) for PCa and BCa cells and also of [⁶⁸Ga]Ga-2 (GRPR), especially in those cell lines already expressing GRPR. The results were confirmed by fluorescence-activated cell sorting with a PE-labeled PSMA-specific antibody and the fluorescence tracer 3. The docking calculations and molecular dynamics simulations showed how 1 enters the PSMA funnel region and how pharmacophore Glu-urea-Lys interacts with the arginine patch, the S1', and S1 subpockets by forming hydrogen and van der Waals bonds. The chelating moiety of 1, that is, HBED-CC, forms additional stabilizing hydrogen bonding and van der Waals interactions in the arene-binding site. Ligand 2 is diving into the GRPR transmembrane (TM) helical cavity, thereby forming hydrogen bonds through its amidated end, water-mediated hydrogen bonds, and π-π interactions. Our results provide valuable information regarding the molecular mechanisms involved in the interactions of 1 and 2 with PSMA and GRPR, which might be useful for the diagnostic imaging and therapy of PCa and BCa.

Modulating Conformational Preferences by Allylic Strain toward Improved Physical Properties and Binding Interactions

The preference of the axial over the equatorial orientation of 2-substitutent for both phenyl-1-piperidines and N-acylpiperidines is studied at the M06-2X level of theory. For phenyl-1-piperidines, the axial 2-substituent is modestly favored over the equatorial one. In contrast, the pseudoallylic strain in N-acylpiperidines dictates the axial orientation of 2-substituent with a ΔG up to -3.2 kcal/mol. The calculations agree well with the statistics from both the Cambridge Structural Database of small-molecule organic crystal structures and the Protein Data Bank. The equilibrium between the twist-boat and chair conformations for N-acylpiperidines with a 2-substituent was further investigated. The twist-boat conformation is found to be around 1.5 kcal/mol less favorable. Finally, the three-dimensionality in shape resulting from minimization of the pseudoallylic strain is characterized, and its implication in protein-ligand interactions is briefly reviewed.

Structural basis of neuropeptide Y signaling through Y1 receptor

Neuropeptide Y (NPY) is highly abundant in the brain and involved in various physiological processes related to food intake and anxiety, as well as human diseases such as obesity and cancer. However, the molecular details of the interactions between NPY and its receptors are poorly understood. Here, we report a cryo-electron microscopy structure of the NPY-bound neuropeptide Y1 receptor (Y1R) in complex with Gi1 protein. The NPY C-terminal segment forming the extended conformation binds deep into the Y1R transmembrane core, where the amidated C-terminal residue Y36 of NPY is located at the base of the ligand-binding pocket. Furthermore, the helical region and two N-terminal residues of NPY interact with Y1R extracellular loops, contributing to the high affinity of NPY for Y1R. The structural analysis of NPY-bound Y1R and mutagenesis studies provide molecular insights into the activation mechanism of Y1R upon NPY binding.

Clinical, Genetic and Functional Characterization of a Novel AVPR2 Missense Mutation in a Woman with X-Linked Recessive Nephrogenic Diabetes Insipidus

Nephrogenic diabetes insipidus (NDI) is a rare disorder characterized by renal unresponsiveness to the hormone vasopressin, leading to excretion of large volumes of diluted urine. Mutations in the arginine vasopressin receptor-2 (AVPR2) gene cause congenital NDI and have an X-linked recessive inheritance. The disorder affects almost exclusively male family members, but female carriers occasionally present partial phenotypes due to skewed inactivation of the X-chromosome. Here, we report a rare case of a woman affected with X-linked recessive NDI, presenting an average urinary output of 12 L/day. Clinical and biochemical studies showed incomplete responses to water deprivation and vasopressin stimulation tests. Genetic analyses revealed a novel heterozygous missense mutation (c.493G > C, p.Ala165Pro) in the AVPR2 gene. Using a combination of in-silico protein modeling with human cellular models and molecular phenotyping, we provide functional evidence for phenotypic effects. The mutation destabilizes the helical structure of the AVPR2 transmembrane domains and disrupts its plasma membrane localization and downstream intracellular signaling pathways upon activation with its agonist vasopressin. These defects lead to deficient aquaporin 2 (AQP2) membrane translocation, explaining the inability to concentrate urine in this patient.

Identification of Highly Selective Orexin 1 Receptor Antagonists Driven by Structure-Based Design

OX1 receptor antagonists are of interest to treat, for example, substance abuse disorders, personality disorders, eating disorders, or anxiety-related disorders. However, known dual OX1/OX2 receptor antagonists are not suitable due to their sleep-inducing effects; therefore, we were interested in identifying a highly OX1 selective antagonist with a sufficient window to OX2-mediated effects. Herein, we describe the design of highly selective OX1 receptor antagonists driven by the X-ray structure of OX1 with suvorexant, a dual OX1/OX2 receptor antagonist. Moderately selective OX1 antagonists comprising a [2.2.1]-bicyclic scaffold served as our starting point. Based on our binding mode hypothesis, we postulated which part of the scaffold points toward one of the regions where the two binding pockets differ. Structural changes in this part resulted in a modified core with higher inherent selectivity compared to the [2.2.1]-bicyclic template. The structure-based design, synthesis, and hit-to-lead evaluation of this novel OX1 receptor-selective scaffold are discussed herein.

Orexin receptors in GtoPdb v.2021.3

Orexin receptors (nomenclature as agreed by the NC-IUPHAR Subcommittee on Orexin receptors [42]) are activated by the endogenous polypeptides orexin-A and orexin-B (also known as hypocretin-1 and -2; 33 and 28 aa) derived from a common precursor, preproorexin or orexin precursor, by proteolytic cleavage and some typical peptide modifications [109]. Currently the only orexin receptor ligands in clinical use are suvorexant and lemborexant, which are used as hypnotics. Orexin receptor crystal structures have been solved [134, 133, 54, 117, 46].

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.

The Novel Application of Geometric Morphometrics with Principal Component Analysis to Existing G Protein-Coupled Receptor (GPCR) Structures

The G protein-coupled receptor (GPCR) superfamily is a large group of membrane proteins which, because of their vast involvement in cell signalling pathways, are implicated in a plethora of disease states and are therefore considered to be key drug targets. Despite advances in techniques to study these receptors, current prophylaxis is often limited due to the challenging nature of their dynamic, complex structures. Greater knowledge and understanding of their intricate structural rearrangements will therefore undoubtedly aid structure-based drug design against GPCRs. Disciplines such as anthropology and palaeontology often use geometric morphometrics to measure variation between shapes and we have therefore applied this technique to analyse GPCR structures in a three-dimensional manner, using principal component analysis. Our aim was to create a novel system able to discriminate between GPCR structures and discover variation between them, correlated with a variety of receptor characteristics. This was conducted by assessing shape changes at the extra- and intracellular faces of the transmembrane helix bundle, analysing the XYZ coordinates of the amino acids at those positions. We have demonstrated that GPCR structures can be classified based on characteristics such as activation state, bound ligands and fusion proteins, with the most significant results focussed at the intracellular face. Conversely, our analyses provide evidence that thermostabilising mutations do not cause significant differences when compared to non-mutated GPCRs. We believe that this is the first time geometric morphometrics has been applied to membrane proteins on this scale, and believe it can be used as a future tool in sense-checking newly resolved structures and planning experimental design.

Molecular Dynamics Simulations for the Determination of the Characteristic Structural Differences between Inactive and Active States of Wild Type and Mutants of the Orexin2 Receptor

The orexin2 receptor (OX2R), which is classified as a class A G protein-coupled receptor (GPCR), is the target of our study. We performed over 20 several-microsecond-scale molecular dynamics simulations of the wild type and mutants of OX2R to extract the characteristics of the structural changes taking place in the active state. We introduced mutations that exhibited the stable inactive state and the constitutively active state in class A GPCRs. In these simulations, significant characteristic structural changes were observed in the V309^6.40Y mutant, which corresponded to a constitutively active mutant. These conformational changes include the outward movement of the transmembrane helix 6 (TM6) and the inward movement of TM7, which are common structural changes in the activation of GPCRs. In addition, we extracted a suitable index for the quantitative evaluation of the active and inactive states of GPCRs, namely, the inter-atomic distance of Cα atoms between x(3.46) and Y(7.53). The structures of the inactive and active states solved by X-ray crystallography and cryo-electron microscopy can be classified using the inter-atomic distance. Furthermore, we clarified that the inward movement of TM7 requires the swapping of M305^6.36 on TM6 and L367^7.56 on TM7. Finally, we discussed the structural advantages of TM7 inward movement for GPCR activation.

Put a ring on it: application of small aliphatic rings in medicinal chemistry

Aliphatic three- and four-membered rings including cyclopropanes, cyclobutanes, oxetanes, azetidines and bicyclo[1.1.1]pentanes have been increasingly exploited in medicinal chemistry for their beneficial physicochemical properties and applications as functional group bioisosteres. This review provides a historical perspective and comparative up to date overview of commonly applied small rings, exemplifying key principles with recent literature examples. In addition to describing the merits and advantages of each ring system, potential hazards and liabilities are also illustrated and explained, including any significant chemical or metabolic stability and toxicity risks.

Structure-Activity Relationship Studies on Oxazolo[3,4-a]pyrazine Derivatives Leading to the Discovery of a Novel Neuropeptide S Receptor Antagonist with Potent In Vivo Activity

Neuropeptide S modulates important neurobiological functions including locomotion, anxiety, and drug abuse through interaction with its G protein-coupled receptor known as neuropeptide S receptor (NPSR). NPSR antagonists are potentially useful for the treatment of substance abuse disorders against which there is an urgent need for new effective therapeutic approaches. Potent NPSR antagonists in vitro have been discovered which, however, require further optimization of their in vivo pharmacological profile. This work describes a new series of NPSR antagonists of the oxazolo[3,4-a]pyrazine class. The guanidine derivative 16 exhibited nanomolar activity in vitro and 5-fold improved potency in vivo compared to SHA-68, a reference pharmacological tool in this field. Compound 16 can be considered a new tool for research studies on the translational potential of the NPSergic system. An in-depth molecular modeling investigation was also performed to gain new insights into the observed structure-activity relationships and provide an updated model of ligand/NPSR interactions.

Structures of active-state orexin receptor 2 rationalize peptide and small-molecule agonist recognition and receptor activation

Narcolepsy type 1 (NT1) is a chronic neurological disorder that impairs the brain’s ability to control sleep-wake cycles. Current therapies are limited to the management of symptoms with modest effectiveness and substantial adverse effects. Agonists of the orexin receptor 2 (OX₂R) have shown promise as novel therapeutics that directly target the pathophysiology of the disease. However, identification of drug-like OX₂R agonists has proven difficult. Here we report cryo-electron microscopy structures of active-state OX₂R bound to an endogenous peptide agonist and a small-molecule agonist. The extended carboxy-terminal segment of the peptide reaches into the core of OX₂R to stabilize an active conformation, while the small-molecule agonist binds deep inside the orthosteric pocket, making similar key interactions. Comparison with antagonist-bound OX₂R suggests a molecular mechanism that rationalizes both receptor activation and inhibition. Our results enable structure-based discovery of therapeutic orexin agonists for the treatment of NT1 and other hypersomnia disorders.

Agonists of the orexin receptor 2 (OX₂R) show promise in the treatment of narcolepsy. Cryo-EM structures of active-state OX₂R bound to an endogenous peptide agonist and a small-molecule agonist suggest a molecular mechanism that rationalizes both receptor activation and inhibition.

Homology Modeling of Class A G-Protein-Coupled Receptors in the Age of the Structure Boom

With 700 members, G protein-coupled receptors (GPCRs) of the rhodopsin family (class A) form the largest membrane receptor family in humans and are the target of about 30% of presently available pharmaceutical drugs. The recent boom in GPCR structures led to the structural resolution of 57 unique receptors in different states (39 receptors in inactive state only, 2 receptors in active state only and 16 receptors in different activation states). In spite of these tremendous advances, most computational studies on GPCRs, including molecular dynamics simulations, virtual screening and drug design, rely on GPCR models obtained by homology modeling. In this protocol, we detail the different steps of homology modeling with the MODELLER software, from template selection to model evaluation. The present structure boom provides closely related templates for most receptors. If, in these templates, some of the loops are not resolved, in most cases, the numerous available structures enable to find loop templates with similar length for equivalent loops. However, simultaneously, the large number of putative templates leads to model ambiguities that may require additional information based on multiple sequence alignments or molecular dynamics simulations to be resolved. Using the modeling of the human bradykinin receptor B1 as a case study, we show how several templates are managed by MODELLER, and how the choice of template(s) and of template fragments can improve the quality of the models. We also give examples of how additional information and tools help the user to resolve ambiguities in GPCR modeling.

Frameshift Variant in Novel Adenosine-A1-Receptor Homolog Associated With Bovine Spastic Syndrome/Late-Onset Bovine Spastic Paresis in Holstein Sires

Since their first description almost 100 years ago, bovine spastic paresis (BSP) and bovine spastic syndrome (BSS) are assumed to be inherited neuronal-progressive diseases in cattle. Affected animals are characterized by (frequent) spasms primarily located in the hind limbs, accompanied by severe pain symptoms and reduced vigor, thus initiating premature slaughter or euthanasia. Due to the late onset of BSP and BSS and the massively decreased lifespan of modern cattle, the importance of these diseases is underestimated. In the present study, BSP/BSS-affected German Holstein breeding sires from artificial insemination centers were collected and pedigree analysis, genome-wide association studies, whole genome resequencing, protein-protein interaction network analysis, and protein-homology modeling were performed to elucidate the genetic background. The analysis of 46 affected and 213 control cattle revealed four significantly associated positions on chromosome 15 (BTA15), i.e., AC_000172.1:g.83465449A>G (-log₁₀P = 19.17), AC_000172.1:g.81871849C>T (-log₁₀P = 8.31), AC_000172.1:g.81872621A>T (-log₁₀P = 6.81), and AC_000172.1:g.81872661G>C (-log₁₀P = 6.42). Two additional loci were significantly associated located on BTA8 and BTA19, i.e., AC_000165.1:g.71177788T>C and AC_000176.1:g.30140977T>G, respectively. Whole genome resequencing of five affected individuals and six unaffected relatives (two fathers, two mothers, a half sibling, and a full sibling) belonging to three different not directly related families was performed. After filtering, a homozygous loss of function variant was identified in the affected cattle, causing a frameshift in the so far unknown gene locus LOC100848076 encoding an adenosine-A1-receptor homolog. An allele frequency of the variant of 0.74 was determined in 3,093 samples of the 1000 Bull Genomes Project.

Membrane protein crystallography in the era of modern structural biology

The aim of structural biology has been always the study of biological macromolecules structures and their mechanistic behaviour at molecular level. To achieve its goal, multiple biophysical methods and approaches have become part of the structural biology toolbox. Considered as one of the pillars of structural biology, X-ray crystallography has been the most successful method for solving three-dimensional protein structures at atomic level to date. It is however limited by the success in obtaining well-ordered protein crystals that diffract at high resolution. This is especially true for challenging targets such as membrane proteins (MPs). Understanding structure-function relationships of MPs at the biochemical level is vital for medicine and drug discovery as they play critical roles in many cellular processes. Though difficult, structure determination of MPs by X-ray crystallography has significantly improved in the last two decades, mainly due to many relevant technological and methodological developments. Today, numerous MP crystal structures have been solved, revealing many of their mechanisms of action. Yet the field of structural biology has also been through significant technological breakthroughs in recent years, particularly in the fields of single particle electron microscopy (cryo-EM) and X-ray free electron lasers (XFELs). Here we summarise the most important advancements in the field of MP crystallography and the significance of these developments in the present era of modern structural biology.

The Quest for the Best Dual Orexin Receptor Antagonist (Daridorexant) for the Treatment of Insomnia Disorders

Since its discovery in 1998, the orexin system has been of interest to the research community as a potential therapeutic target for the treatment of sleep/wake disorders, stress and anxiety disorders, addiction or eating disorders. It consists of two G protein-coupled receptors, the orexin 1 and orexin 2 receptors, and two neuropeptides with agonistic effects, the orexin A and orexin B peptides. Herein we describe our efforts leading to the identification of a promising set of dual orexin receptor antagonists (DORAs) which subsequently went through physiology-based pharmacokinetic and pharmacodynamic modelling>^[1] and finally led to the selection of daridorexant, currently in phase 3 clinical trials for the treatment of insomnia disorders.

Substituted Azabicyclo[2.2.1]heptanes as Selective Orexin-1 Antagonists: Discovery of JNJ-54717793

The orexin system consists of two neuropeptides (orexin-A and orexin-B) that exert their mode of action on two receptors (orexin-1 and orexin-2). While the role of the orexin-2 receptor is established as an important modulator of sleep wake states, the role of the orexin-1 receptor is believed to play a role in addiction, panic, or anxiety. In this manuscript, we describe the optimization of a nonselective substituted azabicyclo[2.2.1]heptane dual orexin receptor antagonist (DORA) into orally bioavailable, brain penetrating, selective orexin-1 receptor (OX1R) antagonists. This resulted in the discovery of our first candidate for clinical development, JNJ-54717793.

Cyclic tailor-made amino acids in the design of modern pharmaceuticals

Tailor-made AAs are indispensable components of modern medicinal chemistry and are becoming increasingly prominent in new drugs. In fact, about 30% of small-molecule pharmaceuticals contain residues of tailor-made AAs or structurally related diamines and amino-alcohols. Cyclic tailor-made AAs present a particular value to rational structural design by virtue of their local conformational constraints and are widely used in lead optimization programs. The present review article highlights 34 compounds, all of which are derived from cyclic AAs, representing recently-approved, small-molecule pharmaceuticals as well as promising drug candidates currently in various phases of clinical study. For each compound, the discussion includes the discovery, therapeutic profile and optimized synthesis, with a focus on the preparation of cyclic tailor-made AA as the principal structural feature. The present review article is intended to serve as a reference source for organic, medicinal and process chemists along with other professionals working in the fields of drug design and pharmaceutical discovery.

Synthesis and biological activity evaluation of azacycloheptane sulfonamide derivatives as potential orexin receptor antagonists

As the orexin signaling system is crucial for the regulation of the sleep/wake cycle, inhibitors of orexin-1 and orexin-2 receptors are of significant interest in the treatment of insomnia. Herein, a series of novel azacycloheptane sulfonamide derivatives were designed and synthesized, and all the compounds were evaluated as potential orexin receptor inhibitors by FLIPR Tetra calcium assay. A majority of the tested azacycloheptane sulfonamide derivatives showed OX1R and OX2R inhibitory activity. Chloro-substituted derivatives functionalized at the C5 or C6 position of the benzoxazole group exhibited better inhibitory activity for OX1R and OX2R than unsubstituted derivatives functionalized at C5 or C6. In addition, phenyl group modification had positive effects on the inhibitory activities, and an electron-withdrawing fluorine group at the ortho or meta position of the phenyl ring improved the OX2R inhibitory activity of the derivatives. This suggests that azacycloheptane sulfonamide derivatives are promising scaffolds for the development of OX1R and OX2R antagonists.

Azacycloheptane sulfonamide derivatives are promising scaffolds for the development of OX1R and OX2R antagonists.

Structure-based development of a subtype-selective orexin 1 receptor antagonist

Orexins are neuropeptides that activate the rhodopsin-like G protein-coupled receptors OX1R and OX2R. The orexin system plays an important role in the regulation of the sleep-wake cycle and the regulation of feeding and emotions. The nonselective orexin receptor antagonist suvorexant has been the first drug on the market targeting the orexin system and is prescribed for the treatment of insomnia. Subtype-selective OX1R antagonists are valuable tools to further investigate the functions and physiological role of the OX1R in vivo and promising lead compounds for the treatment of drug addiction, anxiety, pain or obesity. Starting from the OX1R and OX2R crystal structures bound to suvorexant, we exploited a single amino acid difference in the orthosteric binding site by using molecular docking and structure-based drug design to optimize ligand interactions with the OX1R while introducing repulsive interactions with the OX2R. A newly established enantiospecific synthesis provided ligands showing up to 75-fold selectivity for the OX1R over the OX2R subtype. The structure of a new OX1R antagonist with subnanomolar affinity (JH112) was determined by crystallography in complex with the OX1R and corresponded closely to the docking-predicted geometry. JH112 exhibits high selectivity over a panel of different GPCRs, is able to cross the blood-brain barrier and acts as slowly diffusing and insurmountable antagonist for G_q protein activation and in particular β-arrestin-2 recruitment at OX1R. This study demonstrates the potential of structure-based drug design to develop more subtype-selective GPCR ligands with potentially reduced side effects and provides an attractive probe molecule and lead compound.

Accurate Prediction of GPCR Ligand Binding Affinity with Free Energy Perturbation

The computational prediction of relative binding free energies is a crucial goal for drug discovery, and G protein-coupled receptors (GPCRs) are arguably the most important drug target class. However, they present increased complexity to model compared to soluble globular proteins. Despite breakthroughs, experimental X-ray crystal and cryo-EM structures are challenging to attain, meaning computational models of the receptor and ligand binding mode are sometimes necessary. This leads to uncertainty in understanding ligand-protein binding induced changes such as, water positioning and displacement, side chain positioning, hydrogen bond networks, and the overall structure of the hydration shell around the ligand and protein. In other words, the very elements that define structure activity relationships (SARs) and are crucial for accurate binding free energy calculations are typically more uncertain for GPCRs. In this work we use free energy perturbation (FEP) to predict the relative binding free energies for ligands of two different GPCRs. We pinpoint the key aspects for success such as the important role of key water molecules, amino acid ionization states, and the benefit of equilibration with specific ligands. Initial calculations following typical FEP setup and execution protocols delivered no correlation with experiment, but we show how results are improved in a logical and systematic way. This approach gave, in the best cases, a coefficient of determination (R²) compared with experiment in the range of 0.6-0.9 and mean unsigned errors compared to experiment of 0.6-0.7 kcal/mol. We anticipate that our findings will be applicable to other difficult-to-model protein ligand data sets and be of wide interest for the community to continue improving FE binding energy predictions.

Orexin Antagonists: A Wake-Up Call for Medicinal Chemists?

Orexin receptors are G protein coupled receptors that may be useful targets for sleep disorders, eating disorders, or addictive behavior. Recent work shows that binding of antagonists to these receptors is complex, with strong dependence on hydrophobic hot spots and networks of water-mediated hydrogen bonds. Despite the minimal structural differences between receptor types, selectivity can be achieved in a number of different ways.