Aminergic GPCR-Ligand Interactions: A Chemical and Structural Map of Receptor Mutation Data

AlphaFold accelerated discovery of psychotropic agonists targeting the trace amine-associated receptor 1

Artificial intelligence is revolutionizing protein structure prediction, providing unprecedented opportunities for drug design. To assess the potential impact on ligand discovery, we compared virtual screens using protein structures generated by the AlphaFold machine learning method and traditional homology modeling. More than 16 million compounds were docked to models of the trace amine-associated receptor 1 (TAAR1), a G protein-coupled receptor of unknown structure and target for treating neuropsychiatric disorders. Sets of 30 and 32 highly ranked compounds from the AlphaFold and homology model screens, respectively, were experimentally evaluated. Of these, 25 were TAAR1 agonists with potencies ranging from 12 to 0.03 μM. The AlphaFold screen yielded a more than twofold higher hit rate (60%) than the homology model and discovered the most potent agonists. A TAAR1 agonist with a promising selectivity profile and drug-like properties showed physiological and antipsychotic-like effects in wild-type but not in TAAR1 knockout mice. These results demonstrate that AlphaFold structures can accelerate drug discovery.

A high-affinity, cis-on photoswitchable beta blocker to optically control β2-adrenergic receptors in vitro and in vivo

This study introduces (S)-Opto-prop-2, a second-generation photoswitchable ligand designed for precise modulation of β2-adrenoceptor (β2AR). Synthesised by incorporating an azobenzene moiety with propranolol, (S)-Opto-prop-2 exhibited a high PSScis (photostationary state for cis isomer) percentage (∼90 %) and a favourable half-life (>10 days), facilitating diverse bioassay measurements. In vitro, the cis-isomer displayed substantially higher β2AR binding affinity than the trans-isomer (1000-fold), making (S)-Opto-prop-2 one of the best photoswitchable GPCR (G protein-coupled receptor) ligands reported so far. Molecular docking of (S)-Opto-prop-2 in the X-ray structure of propranolol-bound β2AR followed by site-directed mutagenesis studies, identified D1133.32, N3127.39 and F2896.51 as crucial residues that contribute to ligand-receptor interactions at the molecular level. In vivo efficacy was assessed using a rabbit ocular hypertension model, revealing that the cis isomer mimicked propranolol's effects in reducing intraocular pressure, while the trans isomer was inactive. Dynamic optical modulation of β2AR by (S)-Opto-prop-2 was demonstrated in two different cAMP bioassays and using live-cell confocal imaging, indicating reversible and dynamic control of β2AR activity using the new photopharmacology tool. In conclusion, (S)-Opto-prop-2 emerges as a promising photoswitchable ligand for precise and reversible β2AR modulation with light. The new tool shows superior cis-on binding affinity, one of the largest reported differences in affinity (1000-fold) between its two configurations, in vivo efficacy, and dynamic modulation. This study contributes valuable insights into the evolving field of photopharmacology, offering a potential avenue for targeted therapy in β2AR-associated pathologies.

PI3 kinase inhibitor PI828 uncouples aminergic GPCRs and Ca2+ mobilization irrespectively of its primary target

The phosphoinositide 3-kinase (PI3K) is involved in regulation of multiple intracellular processes. Although the inhibitory analysis is generally employed for validating a physiological role of PI3K, increasing body of evidence suggests that PI3K inhibitors can exhibit PI3K-unrelated activity as well. Here we studied Ca2+ signaling initiated by aminergic agonists in a variety of different cells and analyzed effects of the PI3K inhibitor PI828 on cell responsiveness. It turned out that PI828 inhibited Ca2+ transients elicited by acetylcholine (ACh), histamine, and serotonin, but did not affect Ca2+ responses to norepinephrine and ATP. Another PI3K inhibitor wortmannin negligibly affected Ca2+ signaling initiated by any one of the tested agonists. Using the genetically encoded PIP3 sensor PH(Akt)-Venus, we confirmed that both PI828 and wortmannin effectively inhibited PI3K and ascertained that this kinase negligibly contributed to ACh transduction. These findings suggested that PI828 inhibited Ca2+ responses to aminergic agonists tested, involving an unknown cellular mechanism unrelated to the PI3K inhibition. Complementary physiological experiments provided evidence that PI828 could inhibit Ca2+ signals induced by certain agonists, by acting extracellularly, presumably, through their surface receptors. For the muscarinic M3 receptor, this possibility was verified with molecular docking and molecular dynamics. As demonstrated with these tools, wortmannin could be bound in the extracellular vestibule at the muscarinic M3 receptor but this did not preclude binding of ACh to the M3 receptor followed by its activation. In contrast, PI828 could sterically block the passage of ACh into the allosteric site, preventing activation of the muscarinic M3 receptor.

Structural basis of psychedelic LSD recognition at dopamine D1 receptor

Understanding the kinetics of LSD in receptors and subsequent induced signaling is crucial for comprehending both the psychoactive and therapeutic effects of LSD. Despite extensive research on LSD's interactions with serotonin 2A and 2B receptors, its behavior on other targets, including dopamine receptors, has remained elusive. Here, we present cryo-EM structures of LSD/PF6142-bound dopamine D1 receptor (DRD1)-legobody complexes, accompanied by a β-arrestin-mimicking nanobody, NBA3, shedding light on the determinants of G protein coupling versus β-arrestin coupling. Structural analysis unveils a distinctive binding mode of LSD in DRD1, particularly with the ergoline moiety oriented toward TM4. Kinetic investigations uncover an exceptionally rapid dissociation rate of LSD in DRD1, attributed to the flexibility of extracellular loop 2 (ECL2). Moreover, G protein can stabilize ECL2 conformation, leading to a significant slowdown in ligand's dissociation rate. These findings establish a solid foundation for further exploration of G protein-coupled receptor (GPCR) dynamics and their relevance to signal transduction.

Occurrence of "Natural Selection" in Successful Small Molecule Drug Discovery

Published compounds from ChEMBL version 32 are used to seek evidence for the occurrence of "natural selection" in drug discovery. Three measures of natural product (NP) character were applied, to compare time- and target-matched compounds reaching the clinic (clinical compounds in phase 1-3 development and approved drugs) with background compounds (reference compounds). Pseudo-NPs (PNPs), containing NP fragments combined in ways inaccessible by nature, are increasing over time, reaching 67% of clinical compounds first disclosed since 2010. PNPs are 54% more likely to be found in post-2008 clinical versus reference compounds. The majority of target classes show increased clinical compound NP character versus their reference compounds. Only 176 NP fragments appear in >1000 clinical compounds published since 2008, yet these make up on average 63% of the clinical compound's core scaffolds. There is untapped potential awaiting exploitation, by applying nature's building blocks─"natural intelligence"─to drug design.

Computational insights into diverse binding modes of the allosteric modulator and their regulation on dopamine D1 receptor

Allosteric drugs hold the promise of addressing many challenges in the current drug development of GPCRs. However, the molecular mechanism underlying their allosteric modulations remain largely elusive. The dopamine D1 receptor (DRD1), a member of Class A GPCRs, is critical for treating psychiatric disorders, and LY3154207 serves as its promising positive allosteric modulator (PAM). In the work, we utilized extensive Gaussian-accelerated molecular dynamics simulations (a total of 41μs) for the first time probe the diverse binding modes of the allosteric modulator and their regulation effects, based on the DRD1 and LY3154207 as representative. Our simulations identify four binding modes of LY3154207 (one boat mode, two metastable vertical modes and a novel cleft-anchored mode), in which the boat mode is the most stable while there three modes are similar in the stability. However, it is interesting to observed that the most stable boat mode inversely exhibits the weakest positive allosteric effect on influencing the orthosteric ligand binding and maintaining the activity of the transducer binding site. It should result from its induced weaker correlation between the allosteric site and the orthosteric site, and between the orthosteric site and the transducer binding site than the other three binding modes, as well as its weakened interaction between a crucial activation-related residue (S2025.46) and the orthosteric ligand (dopamine). Overall, the work offers atomic-level information to advance our understanding of the complex allosteric regulation on GPCRs, which is beneficial to the allosteric modulator design and development.

Analysing the effect caused by increasing the molecular volume in M1-AChR receptor agonists and antagonists: a structural and computational study

M1 muscarinic acetylcholine receptor (M1-AChR), a member of the G protein-coupled receptors (GPCR) family, plays a crucial role in learning and memory, making it an important drug target for Alzheimer's disease (AD) and schizophrenia. M1-AChR activation and deactivation have shown modifying effects in AD and PD preclinical models, respectively. However, understanding the pharmacology associated with M1-AChR activation or deactivation is complex, because of the low selectivity among muscarinic subtypes, hampering their therapeutic applications. In this regard, we constructed two quantitative structure-activity relationship (QSAR) models, one for M1-AChR agonists (total and partial), and the other for the antagonists. The binding mode of 59 structurally different compounds, including agonists and antagonists with experimental binding affinity values (pKi), were analyzed employing computational molecular docking over different structures of M1-AChR. Furthermore, we considered the interaction energy (Einter), the number of rotatable bonds (NRB), and lipophilicity (ilogP) for the construction of the QSAR model for agonists (R2 = 89.64, QLMO2 = 78, and Qext2 = 79.1). For the QSAR model of antagonists (R2 = 88.44, QLMO2 = 82, and Qext2 = 78.1) we considered the Einter, the fraction of sp3 carbons fCsp3, and lipophilicity (MlogP). Our results suggest that the ligand volume is a determinant to establish its biological activity (agonist or antagonist), causing changes in binding energy, and determining the affinity for M1-AChR.

Identification of a Novel Subtype-Selective α1B-Adrenoceptor Antagonist

α1A-, α1B-, and α1D-adrenoceptors (α1-ARs) are members of the adrenoceptor G protein-coupled receptor family that are activated by adrenaline (epinephrine) and noradrenaline. α1-ARs are clinically targeted using antagonists that have minimal subtype selectivity, such as prazosin and tamsulosin, to treat hypertension and benign prostatic hyperplasia, respectively. Abundant expression of α1-ARs in the heart and central nervous system (CNS) makes these receptors potential targets for the treatment of cardiovascular and CNS disorders, such as heart failure, epilepsy, and Alzheimer's disease. Our understanding of the precise physiological roles of α1-ARs, however, and their involvement in disease has been hindered by the lack of sufficiently subtype-selective tool compounds, especially for α1B-AR. Here, we report the discovery of 4-[(2-hydroxyethyl)amino]-6-methyl-2H-chromen-2-one (Cpd1), as an α1B-AR antagonist that has 10-15-fold selectivity over α1A-AR and α1D-AR. Through computational and site-directed mutagenesis studies, we have identified the binding site of Cpd1 in α1B-AR and propose the molecular basis of α1B-AR selectivity, where the nonconserved V19745.52 residue plays a major role, with contributions from L3146.55 within the α1B-AR pocket. By exploring the structure-activity relationships of Cpd1 at α1B-AR, we have also identified 3-[(cyclohexylamino)methyl]-6-methylquinolin-2(1H)-one (Cpd24), which has a stronger binding affinity than Cpd1, albeit with reduced selectivity for α1B-AR. Cpd1 and Cpd24 represent potential leads for α1B-AR-selective drug discovery and novel tool molecules to further study the physiology of α1-ARs.

Development of 2-Aminotetralin-Type Serotonin 5-HT1 Agonists: Molecular Determinants for Selective Binding and Signaling at 5-HT1A, 5-HT1B, 5-HT1D, and 5-HT1F Receptors

The serotonin (5-hydroxytryptamine, 5-HT) 5-HT1 G-protein coupled receptor subtypes (5-HT1A/1B/1D/1E/1F) share a high sequence homology, confounding development of subtype-specific ligands. This study used a 5-HT1 structure-based ligand design approach to develop subtype-selective ligands using a 5-substituted-2-aminotetralin (5-SAT) chemotype, leveraging results from pharmacological, molecular modeling, and mutagenesis studies to delineate molecular determinants for 5-SAT binding and function at 5-HT1 subtypes. 5-SATs demonstrated high affinity (Ki ≤ 25 nM) and at least 50-fold stereoselective preference ([2S] > [2R]) at 5-HT1A, 5-HT1B, and 5-HT1D receptors but essentially nil affinity (Ki > 1 μM) at 5-HT1F receptors. The 5-SATs tested were agonists with varying degrees of potency and efficacy, depending on chemotype substitution and 5-HT1 receptor subtype. Models were built from the 5-HT1A (cryo-EM), 5-HT1B (crystal), and 5-HT1D (cryo-EM) structures, and 5-SATs underwent docking studies with up to 1 μs molecular dynamics simulations. 5-SAT interactions observed at positions 3.33, 5.38, 5.42, 5.43, and 7.39 of 5-HT1 subtypes were confirmed with point mutation experiments. Additional 5-SATs were designed and synthesized to exploit experimental and computational results, yielding a new full efficacy 5-HT1A agonist with 100-fold selectivity over 5-HT1B/1D receptors. The results presented lay the foundation for the development of additional 5-HT1 subtype selective ligands for drug discovery purposes.

Structure-Based Discovery of Mouse Trace Amine-Associated Receptor 5 Antagonists

Trace amine-associated receptors (TAARs) were discovered in 2001 as new members of class A G protein-coupled receptors (GPCRs). With the only exception of TAAR1, TAAR members (TAAR2-9, also known as noncanonical olfactory receptors) were originally described exclusively in the olfactory epithelium and believed to mediate the innate perception of volatile amines. However, most noncanonical olfactory receptors are still orphan receptors. Given its recently discovered nonolfactory expression and therapeutic potential, TAAR5 has been the focus of deorphanization campaigns that led to the discovery of a few druglike antagonists. Here, we report four novel TAAR5 antagonists identified through high-throughput screening, which, along with the four ligands published in the literature, constituted our starting point to design a computational strategy for the identification of TAAR5 ligands. We developed a structure-based virtual screening protocol that allowed us to identify three new TAAR5 antagonists with a hit rate of 10%. Despite lacking an experimental structure, we accurately modeled the TAAR5 binding site by integrating comparative sequence- and structure-based analyses of serotonin receptors with homology modeling and side-chain optimization. In summary, we have identified seven new TAAR5 antagonists that could serve as lead candidates for the development of new treatments for depression, anxiety, and neurodegenerative diseases.

Ligand binding affinity prediction with fusion of graph neural networks and 3D structure-based complex graph

Accurate prediction of protein-ligand binding affinity is pivotal for drug design and discovery. Here, we proposed a novel deep fusion graph neural networks framework named FGNN to learn the protein-ligand interactions from the 3D structures of protein-ligand complexes. Unlike 1D sequences for proteins or 2D graphs for ligands, the 3D graph of protein-ligand complex enables the more accurate representations of the protein-ligand interactions. Benchmark studies have shown that our fusion models FGNN can achieve more accurate prediction of binding affinity than any individual algorithm. The advantages of fusion strategies have been demonstrated in terms of expressive power of data, learning efficiency and model interpretability. Our fusion models show satisfactory performances on diverse data sets, demonstrating their generalization ability. Given the good performances in both binding affinity prediction and virtual screening, our fusion models are expected to be practically applied for drug screening and design. Our work highlights the potential of the fusion graph neural network algorithm in solving complex prediction problems in computational biology and chemistry. The fusion graph neural networks (FGNN) model is freely available in https://github.com/LinaDongXMU/FGNN.

Structure-Activity Assessment and In-Depth Analysis of Biased Agonism in a Set of Phenylalkylamine 5-HT2A Receptor Agonists

Serotonergic psychedelics are described to have activation of the serotonin 2A receptor (5-HT2A) as their main pharmacological action. Despite their relevance, the molecular mechanisms underlying the psychedelic effects induced by certain 5-HT2A agonists remain elusive. One of the proposed hypotheses is the occurrence of biased agonism, defined as the preferential activation of certain signaling pathways over others. This study comparatively monitored the efficiency of a diverse panel of 4-position-substituted (and N-benzyl-derived) phenylalkylamines to induce recruitment of β-arrestin2 (βarr2) or miniGαq to the 5-HT2A, allowing us to assess structure-activity relationships and biased agonism. All test compounds exhibited agonist properties with a relatively large range of both EC50 and Emax values. Interestingly, the lipophilicity of the 2C-X phenethylamines was correlated with their efficacy in both assays but yielded a stronger correlation in the miniGαq- than in the βarr2-assay. Molecular docking suggested that accommodation of the 4-substituent of the 2C-X analogues in a hydrophobic pocket between transmembrane helices 4 and 5 of 5-HT2A may contribute to this differential effect. Aside from previously used standard conditions (lysergic acid diethylamide (LSD) as a reference agonist and a 2 h activation profile to assess a compound's activity), serotonin was included as a second reference agonist, and the compounds' activities were also assessed using the first 30 min of the activation profile. Under all assessed circumstances, the qualitative structure-activity relationships remained unchanged. Furthermore, the use of two reference agonists allowed for the estimation of both "benchmark bias" (relative to LSD) and "physiology bias" (relative to serotonin).

Structural basis of α1A-adrenergic receptor activation and recognition by an extracellular nanobody

The α_1A-adrenergic receptor (α_1AAR) belongs to the family of G protein-coupled receptors that respond to adrenaline and noradrenaline. α_1AAR is involved in smooth muscle contraction and cognitive function. Here, we present three cryo-electron microscopy structures of human α_1AAR bound to the endogenous agonist noradrenaline, its selective agonist oxymetazoline, and the antagonist tamsulosin, with resolutions range from 2.9 Å to 3.5 Å. Our active and inactive α_1AAR structures reveal the activation mechanism and distinct ligand binding modes for noradrenaline compared with other adrenergic receptor subtypes. In addition, we identified a nanobody that preferentially binds to the extracellular vestibule of α_1AAR when bound to the selective agonist oxymetazoline. These results should facilitate the design of more selective therapeutic drugs targeting both orthosteric and allosteric sites in this receptor family.

Insight into the Mode of Action of 8-Hydroxyquinoline-Based Blockers on the Histamine Receptor 2

Histamine receptor 2 (HR_H2) blockers are used to treat peptic ulcers and gastric reflux. Chlorquinaldol and chloroxine, which contain an 8-hydroxyquinoline (8HQ) core, have recently been identified as blocking HR_H2. To gain insight into the mode of action of 8HQ-based blockers, here, we leverage an HR_H2-based sensor in yeast to evaluate the role of key residues in the HR_H2 active site on histamine and 8HQ-based blocker binding. We find that the HR_H2 mutations D98A, F254A, Y182A, and Y250A render the receptor inactive in the presence of histamine, while HR_H2:D186A and HR_H2:T190A retain residual activity. Based on molecular docking studies, this outcome correlates with the ability of the pharmacologically relevant histamine tautomers to interact with D98 via the charged amine. Docking studies also suggest that, unlike established HR_H2 blockers that interact with both ends of the HR_H2 binding site, 8HQ-based blockers interact with only one end, either the end framed by D98/Y250 or T190/D186. Experimentally, we find that chlorquinaldol and chloroxine still inactivate HR_H2:D186A by shifting their engagement from D98 to Y250 in the case of chlorquinaldol and D186 to Y182 in the case of chloroxine. Importantly, the tyrosine interactions are supported by the intramolecular hydrogen bonding of the 8HQ-based blockers. The insight gained in this work will aid in the development of improved HR_H2 therapeutics. More generally, this work demonstrates that Gprotein-coupled receptor (GPCR)-based sensors in yeast can help elucidate the mode of action of novel ligands for GPCRs, a family of receptors that bind 30% of FDA therapeutics.

Dopamine Receptor Ligand Selectivity-An In Silico/In Vitro Insight

Different dopamine receptor (DR) subtypes are involved in pathophysiological conditions such as Parkinson's Disease (PD), schizophrenia and depression. While many DR-targeting drugs have been approved by the U.S. Food and Drug Administration (FDA), only a very small number are truly selective for one of the DR subtypes. Additionally, most of them show promiscuous activity at related G-protein coupled receptors, thus suffering from diverse side-effect profiles. Multiple studies have shown that combined in silico/in vitro approaches are a valuable contribution to drug discovery processes. They can also be applied to divulge the mechanisms behind ligand selectivity. In this study, novel DR ligands were investigated in vitro to assess binding affinities at different DR subtypes. Thus, nine D2R/D3R-selective ligands (micro- to nanomolar binding affinities, D3R-selective profile) were successfully identified. The most promising ligand exerted nanomolar D3R activity (Ki = 2.3 nM) with 263.7-fold D2R/D3R selectivity. Subsequently, ligand selectivity was rationalized in silico based on ligand interaction with a secondary binding pocket, supporting the selectivity data determined in vitro. The developed workflow and identified ligands could aid in the further understanding of the structural motifs responsible for DR subtype selectivity, thus benefitting drug development in D2R/D3R-associated pathologies such as PD.

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.

Structural basis for recognition of antihistamine drug by human histamine receptor

The histamine receptors belong to the G protein-coupled receptor (GPCR) superfamily, and play important roles in the regulation of histamine and other neurotransmitters in the central nervous system, as potential targets for the treatment of neurologic and psychiatric disorders. Here we report the crystal structure of human histamine receptor H₃R bound to an antagonist PF-03654746 at 2.6 Å resolution. Combined with the computational and functional assays, our structure reveals binding modes of the antagonist and allosteric cholesterol. Molecular dynamic simulations and molecular docking of different antihistamines further elucidate the conserved ligand-binding modes. These findings are therefore expected to facilitate the structure-based design of novel antihistamines.

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 105 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.

Evolutionary association of receptor-wide amino acids with G protein-coupling selectivity in aminergic GPCRs

G protein-coupled receptors (GPCRs) induce signal transduction pathways through coupling to four main subtypes of G proteins (Gs, Gi, Gq, and G12/13), selectively. However, G protein selective activation mechanisms and residual determinants in GPCRs have remained obscure. Herein, we performed extensive phylogenetic analysis and identified specifically conserved residues for the aminergic receptors having similar coupling profiles. By integrating our methodology of differential evolutionary conservation of G protein-specific amino acids with structural analyses, we identified specific activation networks for Gs, Gi1, Go, and Gq To validate that these networks could determine coupling selectivity we further analyzed Gs-specific activation network and its association with Gs selectivity. Through molecular dynamics simulations, we showed that previously uncharacterized Glycine at position 7x41 plays an important role in receptor activation and it may determine Gs coupling selectivity by facilitating a larger TM6 movement. Finally, we gathered our results into a comprehensive model of G protein selectivity called "sequential switches of activation" describing three main molecular switches controlling GPCR activation: ligand binding, G protein selective activation mechanisms, and G protein contact.

Trends in small molecule drug properties: A developability molecule assessment perspective

Developability molecule assessment is a key interfacial capability across the biopharmaceutical industry, screening and staging molecules discovered by medicinal chemists for successful chemistry manufacturing controls (CMC) development and launch. The breadth of responsibility and expertise such teams possess puts them in a unique position to understand the impact of the physicochemical properties of a drug during its initial discovery and subsequent development. However, most of the publications describing trends in physicochemical properties are written from a medicinal chemistry perspective with the aim to identify molecules with better ADMET profiles that are either lead-like or drug-like, failing to describe the impact these properties have on CMC development. To systematically uncover knowledge obtained from recent trends in physicochemical properties and the corresponding impact on CMC development, a comprehensive analysis was conducted on molecules in the drug repurposing hub dataset. The only physicochemical property that seems to have been preserved in FDA-approved oral molecules over the decades (1900-2020) is a constant H-bond donor count, highlighting the importance this property has on cell permeability and lattice energy. Pharmaceutical attrition analysis suggests that partition-distribution coefficient, H-bond acceptors, polar surface area and the fraction of sp³ carbons are properties that are associated with compound attrition. Looking at pharmaceutical attrition asynchronously with the temporal analysis of FDA-approved oral molecules highlights the opposing trends, risks and diminishing effects some of these physiochemical properties (cLogP, cLogD and Fsp3) have on describing compound attrition during the past decade. Trellising the dataset by target class suggests that certain formulation and drug delivery strategies can be anticipated or put into place based on target class of a molecule. For example, molecules binding to nuclear hormone receptors are amenable to lipid-based drug delivery systems with proven commercial success. Although the poor solubility of kinase inhibitors is a combination of hydrophobicity (due to aromaticity) required to bind to its target and high lattice energy (melting point), they are a challenging target class to formulate. The influence of drug targets on physicochemical properties and the temporal nature of these properties is highlighted when comparing molecules in the drug repurposing dataset to those developed at Amgen. An improved understanding of the impact of molecular properties on performance attributes can accelerate decisions and facilitate risk assessments during candidate selection and development.

Discovery of 8-Hydroxyquinoline as a Histamine Receptor 2 Blocker Scaffold

Histamine receptor 2 (HR_H2) activation in the stomach results in gastric acid secretion, and HR_H2 blockers are used for the treatment of peptidic ulcers and acid reflux. Over-the-counter HR_H2 blockers carry a five-membered aromatic heterocycle, with two of them additionally carrying a tertiary amine that decomposes to N-nitrosodimethylamine, a human carcinogen. To discover a novel HR_H2 blocker scaffold to serve in the development of next-generation HR_H2 blockers, we developed an HR_H2-based sensor in yeast by linking human HR_H2 activation to cell luminescence. We used the HR_H2-based sensor to screen a 403-member anti-infection chemical library and identified three HR_H2 blockers, chlorquinaldol, chloroxine, and broxyquinoline, all sharing an 8-hydroxyquinoline scaffold, which is not found among known HR_H2 antagonists. Critically, we validate their HR_H2-blocking ability in mammalian cells. Molecular docking suggests that the HR_H2 blockers bind the histamine binding pocket and structure–activity data point toward these blockers acting as competitive antagonists. Chloroxine and broxyquinoline are antimicrobials that can be found in the gastrointestinal tract at concentrations that would block HR_H2, thus likely modulating gastric acid secretion. Taken together, this work demonstrates the utility of GPCR-based sensors for rapid drug discovery applications, identifies a novel HR_H2 blocker scaffold, and provides further evidence that antimicrobials not only target the human microbiota but also the human host.

GPCRs steer Gi and Gs selectivity via TM5-TM6 switches as revealed by structures of serotonin receptors

Serotonin (or 5-hydroxytryptamine, 5-HT) is an important neurotransmitter that activates 12 different G protein-coupled receptors (GPCRs) through selective coupling of G_s, G_i, or G_q proteins. The structural basis for G protein subtype selectivity by these GPCRs remains elusive. Here, we report the structures of the serotonin receptors 5-HT₄, 5-HT₆, and 5-HT₇ with G_s, and 5-HT₄ with G_i1. The structures reveal that transmembrane helices TM5 and TM6 alternate lengths as a macro-switch to determine receptor's selectivity for G_s and G_i, respectively. We find that the macro-switch by the TM5-TM6 length is shared by class A GPCR-G protein structures. Furthermore, we discover specific residues within TM5 and TM6 that function as micro-switches to form specific interactions with G_s or G_i. Together, these results present a common mechanism of G_s versus G_i protein coupling selectivity or promiscuity by class A GPCRs and extend the basis of ligand recognition at serotonin receptors.

"Selective" serotonin 5-HT2A receptor antagonists

Blockade of the serotonin 5-HT2A G protein-coupled receptor (5-HT2AR) is a fundamental pharmacological characteristic of numerous antipsychotic medications, which are FDA-approved to treat schizophrenia, bipolar disorder, and as adjunctive therapies in major depressive disorder. Meanwhile, activation of the 5-HT2AR by serotonergic psychedelics may be useful in treating neuropsychiatric indications, including major depressive and substance use disorders. Serotonergic psychedelics and other 5-HT2AR agonists, however, often bind other receptors, and standard 5-HT2AR antagonists lack sufficient selectivity to make well-founded mechanistic conclusions about the 5-HT2AR-dependent effects of these compounds and the general neurobiological function of 5-HT2ARs. This review discusses the limitations and strengths of currently available "selective" 5-HT2AR antagonists, the molecular determinants of antagonist selectivity at 5-HT2ARs, and the utility of molecular pharmacology and computational methods in guiding the discovery of novel unambiguously selective 5-HT2AR antagonists.

Serotonin G Protein-Coupled Receptor-Based Biosensing Modalities in Yeast

Serotonin is a key neurotransmitter involved in numerous physiological processes and serves as an important precursor for manufacturing bioactive indoleamines and alkaloids used in the treatment of human pathologies. In humans, serotonin sensing and signaling can occur by 12 G protein-coupled receptors (GPCRs) coupled to Gα proteins. In yeast, human serotonin GPCRs coupled to Gα proteins have previously been shown to function as whole-cell biosensors of serotonin. However, systematic characterization of serotonin biosensing modalities between variant serotonin GPCRs and application thereof for high-resolution serotonin quantification is still awaiting. To systematically assess GPCR signaling in response to serotonin, we characterized reporter gene expression at two different pHs of a 144-sized library encoding all 12 human serotonin GPCRs in combination with 12 different Gα proteins engineered in yeast. From this screen, we observed changes in the biosensor sensitivities of >4 orders of magnitude. Furthermore, adopting optimal biosensing designs and pH conditions enabled high-resolution high-performance liquid chromatography-validated sensing of serotonin produced in yeast. Lastly, we used the yeast platform to characterize 19 serotonin GPCR polymorphisms found in human populations. While major differences in signaling were observed among the individual polymorphisms when studied in yeast, a cross-comparison of selected variants in mammalian cells showed both similar and disparate results. Taken together, our study highlights serotonin biosensing modalities of relevance to both biotechnological and potential human health applications.

Using GPCRs as Molecular Beacons to Target Ovarian Cancer with Nanomedicines

The five-year survival rate for women with ovarian cancer is very poor despite radical cytoreductive surgery and chemotherapy. Although most patients initially respond to platinum-based chemotherapy, the majority experience recurrence and ultimately develop chemoresistance, resulting in fatal outcomes. The current administration of cytotoxic compounds is hampered by dose-limiting severe adverse effects. There is an unmet clinical need for targeted drug delivery systems that transport chemotherapeutics selectively to tumor cells while minimizing off-target toxicity. G protein-coupled receptors (GPCRs) are the largest family of membrane receptors, and many are overexpressed in solid tumors, including ovarian cancer. This review summarizes the progress in engineered nanoparticle research for drug delivery for ovarian cancer and discusses the potential use of GPCRs as molecular entry points to deliver anti-cancer compounds into ovarian cancer cells. A newly emerging treatment paradigm could be the personalized design of nanomedicines on a case-by-case basis.

Binding and Activation of Serotonergic G-Protein Coupled Receptors by the Multimodal Antidepressant Vortioxetine

G-protein coupled receptors (GPCRs) are important pharmacological targets. Despite substantial progress, important questions still remain concerning the details of activation: how can a ligand act as an agonist in one receptor but as an antagonist in a homologous receptor, and how can agonists activate a receptor despite lacking polar functional groups able to interact with helix 5 as is the case for the related adrenergic receptors? Studying vortioxetine (VXT), an important multimodal antidepressant drug, may elucidate both questions. Herein, we present a thorough in silico analysis of VXT binding to 5-HT_1A, 5-HT_1B, and 5-HT₇ receptors and compare it with available experimental data. We are able to rationalize the differential mode of action of VXT at different receptors, but also, in the case of the 5-HT_1A receptor, we observe the initial steps of activation that inform about an activation mechanism that does not involve polar interaction with helix 5. The results extend our current understanding of agonist and antagonist action at aminergic GPCRs.

Ligand-directed bias of G protein signaling at the dopamine D2 receptor

G-protein-coupled receptors (GPCRs) represent the largest family of drug targets. Upon activation, GPCRs signal primarily via a diverse set of heterotrimeric G proteins. Most GPCRs can couple to several different G protein subtypes. However, how drugs act at GPCRs contributing to the selectivity of G protein recognition is poorly understood. Here, we examined the G protein selectivity profile of the dopamine D₂ receptor (D₂), a GPCR targeted by antipsychotic drugs. We show that D₂ discriminates between six individual members of the G_i/o family, and its profile of functional selectivity is remarkably different across its ligands, which all engaged D₂ with a distinct G protein coupling pattern. Using structural modeling, receptor mutagenesis, and pharmacological evaluation, we identified residues in the D₂ binding pocket that shape these ligand-directed biases. We further provide pharmacogenomic evidence that natural variants in D₂ differentially affect its G protein biases in response to different ligands.

Two entry tunnels in mouse TAAR9 suggest the possibility of multi-entry tunnels in olfactory receptors

Orthosteric binding sites of olfactory receptors have been well understood for ligand-receptor interactions. However, a lack of explanation for subtle differences in ligand profile of olfactory receptors even with similar orthosteric binding sites promotes more exploration into the entry tunnels of the receptors. An important question regarding entry tunnels is the number of entry tunnels, which was previously believed to be one. Here, we used TAAR9 that recognizes important biogenic amines such as cadaverine, spermine, and spermidine as a model for entry tunnel study. We identified two entry tunnels in TAAR9 and described the residues that form the tunnels. In addition, we found two vestibular binding pockets, each located in one tunnel. We further confirmed the function of two tunnels through site-directed mutagenesis. Our study challenged the existing views regarding the number of entry tunnels in the subfamily of olfactory receptors and demonstrated the possible mechanism how the entry tunnels function in odorant recognition.

Crystal structure of the α1B-adrenergic receptor reveals molecular determinants of selective ligand recognition

α-adrenergic receptors (αARs) are G protein-coupled receptors that regulate vital functions of the cardiovascular and nervous systems. The therapeutic potential of αARs, however, is largely unexploited and hampered by the scarcity of subtype-selective ligands. Moreover, several aminergic drugs either show off-target binding to αARs or fail to interact with the desired subtype. Here, we report the crystal structure of human α_1BAR bound to the inverse agonist (+)-cyclazosin, enabled by the fusion to a DARPin crystallization chaperone. The α_1BAR structure allows the identification of two unique secondary binding pockets. By structural comparison of α_1BAR with α₂ARs, and by constructing α_1BAR-α_2CAR chimeras, we identify residues 3.29 and 6.55 as key determinants of ligand selectivity. Our findings provide a basis for discovery of α_1BAR-selective ligands and may guide the optimization of aminergic drugs to prevent off-target binding to αARs, or to elicit a selective interaction with the desired subtype.

This study reports the X-ray structure of the α_1B-adrenergic G protein-coupled receptor bound to an inverse agonist, and unveils key determinants of subtype-selective ligand binding that may help the design of aminergic drugs with fewer side-effects.

Structure-based design of a novel third-generation antipsychotic drug lead with potential antidepressant properties

Partial agonist activity at the dopamine D₂ receptor (DRD2) is a key feature of third-generation antipsychotics (TGAs). However, TGAs also act as antagonists or weak partial agonists to the serotonin (5-hydroxytryptamine; 5-HT) 2A receptor (5-HT_2AR). Here we present the crystal structures of aripiprazole- and cariprazine-bound human 5-HT_2AR. Both TGAs adopt an unexpected 'upside-down' pose in the 5-HT_2AR binding pocket, with secondary pharmacophores inserted in a similar way to a 'bolt'. This insight into the binding modes of TGAs offered a structural mechanism underlying their varied partial efficacies at 5-HT_2AR and DRD2. These structures enabled the design of a partial agonist at DRD2/3 and 5-HT_1AR with negligible 5-HT_2AR binding that displayed potent antipsychotic-like activity without motor side effects in mice. This TGA lead also had antidepressant-like effects and improved cognitive performance in mouse models via 5-HT_1AR. This work indicates that 5-HT_2AR affinity is a dispensable contributor to the therapeutic actions of TGAs.

A new class of serotonin 5-HT2A /5-HT2C receptor inverse agonists: Synthesis, molecular modeling, in vitro and in vivo pharmacology of novel 2-aminotetralins

BACKGROUND AND PURPOSE

The 5-HT receptor (5-HTR) subtypes 5-HT_2A and 5-HT_2C are important neurotherapeutic targets, though, obtaining selectivity over 5-HT_2B and closely related histamine H₁ Rs is challenging. Here, we delineated molecular determinants of selective binding to 5-HT_2A and 5-HT_2C Rs for novel 4-phenyl-2-dimethylaminotetralins (4-PATs).

EXPERIMENTAL APPROACH

We synthesized 42 novel 4-PATs with halogen or aryl moieties at the C(4)-phenyl meta position. Affinity, function, molecular modeling, and 5-HT_2A R mutagenesis studies were undertaken to understand structure-activity relationships at 5-HT₂ -type and H₁ Rs. Lead 4-PAT-type selective 5-HT_2A /5-HT_2C R inverse agonists were compared to pimavanserin, a selective 5-HT_2A /5-HT_2C R inverse agonist approved to treat psychoses, in the mouse head twitch response, and locomotor activity assays, as models relevant to antipsychotic drug development.

KEY RESULTS

Most 4-PAT diastereomers in the (2S,4R)-configuration bound non-selectively to 5-HT_2A , 5-HT_2C, and H₁ Rs, with >100-fold selectivity over 5-HT_2B Rs, whereas, diastereomers in the (2R,4R)-configuration bound preferentially to 5-HT_2A over 5-HT_2C Rs and had >100-fold selectivity over 5-HT_2B and H₁ Rs. Results suggest that G238^5.42 and V235^5.39 in 5-HT_2A Rs (conserved in 5-HT_2C Rs) are important for high affinity binding, whereas, interactions with T194^5.42 and W158^4.56 determine H₁ R affinity. The 4-PAT (2S,4R)-2k, a potent and selective 5-HT_2A /5-HT_2C R inverse agonist, had activity like pimavanserin in the mouse head-twitch response assay, but was distinct in not suppressing locomotor activity.

CONCLUSIONS AND IMPLICATIONS

We provide evidence that the novel 4-PAT chemotype can yield selective 5-HT_2A /5-HT_2C R inverse agonists for antipsychotic drug development by optimizing ligand-receptor interactions in transmembrane domain 5. We also show that chirality can be exploited to attain selectivity over H₁ Rs which may circumvent sedative effects.

From structure to clinic: Design of a muscarinic M1 receptor agonist with potential to treatment of Alzheimer's disease

Current therapies for Alzheimer's disease seek to correct for defective cholinergic transmission by preventing the breakdown of acetylcholine through inhibition of acetylcholinesterase, these however have limited clinical efficacy. An alternative approach is to directly activate cholinergic receptors responsible for learning and memory. The M1-muscarinic acetylcholine (M1) receptor is the target of choice but has been hampered by adverse effects. Here we aimed to design the drug properties needed for a well-tolerated M1-agonist with the potential to alleviate cognitive loss by taking a stepwise translational approach from atomic structure, cell/tissue-based assays, evaluation in preclinical species, clinical safety testing, and finally establishing activity in memory centers in humans. Through this approach, we rationally designed the optimal properties, including selectivity and partial agonism, into HTL9936-a potential candidate for the treatment of memory loss in Alzheimer's disease. More broadly, this demonstrates a strategy for targeting difficult GPCR targets from structure to clinic.

Data-Driven Analysis of Fluorination of Ligands of Aminergic G Protein Coupled Receptors

Currently, G protein-coupled receptors are the targets with the highest number of drugs in many therapeutic areas. Fluorination has become a common strategy in designing highly active biological compounds, as evidenced by the steadily increasing number of newly approved fluorine-containing drugs. Herein, we identified in the ChEMBL database and analysed 1554 target-based FSAR sets (non-fluorinated compounds and their fluorinated analogues) comprising 966 unique non-fluorinated and 2457 unique fluorinated compounds active against 33 different aminergic GPCRs. Although a relatively small number of activity cliffs (defined as a pair of structurally similar compounds showing significant differences of activity -ΔpPot > 1.7) was found in FSAR sets, it is clear that appropriately introduced fluorine can increase ligand potency more than 50-fold. The analysis of matched molecular pairs (MMPs) networks indicated that the fluorination of the aromatic ring showed no clear trend towards a positive or negative effect on affinity; however, a favourable site for a positive potency effect of fluorination was the ortho position. Fluorination of aliphatic fragments more often led to a decrease in biological activity. The results may constitute the rules of thumb for fluorination of aminergic receptor ligands and provide insights into the role of fluorine substitutions in medicinal chemistry.

Convergent olfactory trace amine-associated receptors detect biogenic polyamines with distinct motifs via a conserved binding site

Biogenic amines activate G-protein-coupled receptors (GPCRs) in the central nervous system in vertebrate animals. Several biogenic amines, when excreted, stimulate trace amine-associated receptors (TAARs), a group of GPCRs in the main olfactory epithelium, and elicit innate behaviors. How TAARs recognize amines with varying numbers of amino groups is largely unknown. We reasoned that a comparison between lamprey and mammalian olfactory TAARs, which are thought to have evolved independently but show convergent responses to polyamines, may reveal structural determinants of amine recognition. Here, we demonstrate that sea lamprey TAAR365 (sTAAR365) responds strongly to biogenic polyamines cadaverine, putrescine, and spermine, and shares a similar response profile as a mammalian TAAR (mTAAR9). Docking and site-directed mutagenesis analyses show that both sTAAR365 and mTAAR9 recognize the two amino groups of cadaverine with the conserved Asp^3.32 and Tyr^6.51 residues. sTAAR365, which has remarkable sensitivity for cadaverine (EC₅₀ = 4 nM), uses an extra residue, Thr^7.42, to stabilize ligand binding. These cadaverine recognition sites also interact with amines with four and three amino groups (spermine and spermidine, respectively). Glu^7.36 of sTAAR365 cooperates with Asp^3.32 and Thr^7.42 to recognize spermine, whereas mTAAR9 recognizes spermidine through an additional aromatic residue, Tyr^7.43. These results suggest a conserved mechanism whereby independently evolved TAAR receptors recognize amines with two, three, or four amino groups using the same recognition sites, at which sTAAR365 and mTAAR9 evolved distinct motifs. These motifs interact directly with the amino groups of the polyamines, a class of potent and ecologically important odorants, mediating olfactory signaling.

Structure-Guided Design of G-Protein-Coupled Receptor Polypharmacology

Many diseases are polygenic and can only be treated efficiently with drugs that modulate multiple targets. However, rational design of compounds with multi-target profiles is rarely pursued because it is considered too difficult, in particular if the drug must enter the central nervous system. Here, a structure-based strategy to identify dual-target ligands of G-protein-coupled receptors is presented. We use this approach to design compounds that both antagonize the A_2A adenosine receptor and activate the D₂ dopamine receptor, which have excellent potential as antiparkinson drugs. Atomic resolution models of the receptors guided generation of a chemical library with compounds designed to occupy orthosteric and secondary binding pockets in both targets. Structure-based virtual screens identified ten compounds, of which three had affinity for both targets. One of these scaffolds was optimized to nanomolar dual-target activity and showed the predicted pharmacodynamic effect in a rat model of Parkinsonism.

Guanidine Derivatives: How Simple Structural Modification of Histamine H3R Antagonists Has Led to the Discovery of Potent Muscarinic M2R/M4R Antagonists

This article describes the discovery of novel potent muscarinic receptor antagonists identified during a search for more active histamine H₃ receptor (H₃R) ligands. The idea was to replace the flexible seven methylene linker with a semirigid 1,4-cyclohexylene or p-phenylene substituted group of the previously described histamine H₃R antagonists ADS1017 and ADS1020. These simple structural modifications of the histamine H₃R antagonist led to the emergence of additional pharmacological effects, some of which unexpectedly showed strong antagonist potency at muscarinic receptors. This paper reports the routes of synthesis and pharmacological characterization of guanidine derivatives, a novel chemotype of muscarinic receptor antagonists binding to the human muscarinic M₂ and M₄ receptors (hM₂R and hM₄R, respectively) in nanomolar concentration ranges. The affinities of the newly synthesized ADS10227 (1-{4-{4-{[4-(phenoxymethyl)cyclohexyl]methyl}piperazin-1-yl}but-1-yl}-1-(benzyl)guanidine) at hM₂R and hM₄R were 2.8 nM and 5.1 nM, respectively.

Abolishing Dopamine D2long/D3 Receptor Affinity of Subtype-Selective Carbamoylguanidine-Type Histamine H2 Receptor Agonists

3-(2-Amino-4-methylthiazol-5-yl)propyl-substituted carbamoylguanidines are potent, subtype-selective histamine H₂ receptor (H₂R) agonists, but their applicability as pharmacological tools to elucidate the largely unknown H₂R functions in the central nervous system (CNS) is compromised by their concomitant high affinity toward dopamine D₂-like receptors (especially to the D₃R). To improve the selectivity, a series of novel carbamoylguanidine-type ligands containing various heterocycles, spacers, and side residues were rationally designed, synthesized, and tested in binding and/or functional assays at H_1-4 and D_2long/3 receptors. This study revealed a couple of selective candidates (among others 31 and 47), and the most promising ones were screened at several off-target receptors, showing good selectivities. Docking studies suggest that the amino acid residues (3.28, 3.32, E2.49, E2.51, 5.42, and 7.35) are responsible for the different affinities at the H₂- and D_2long/3-receptors. These results provide a solid base for the exploration of the H₂R functions in the brain in further studies.

Controlling the selectivity of aminergic GPCR ligands from the extracellular vestibule

In addition to the orthosteric binding pocket (OBP) of GPCRs, recent structural studies have revealed that there are several allosteric sites available for pharmacological intervention. The secondary binding pocket (SBP) of aminergic GPCRs is located in the extracellular vestibule of these receptors, and it has been suggested to be a potential selectivity pocket for bitopic ligands. Here, we applied a virtual screening protocol based on fragment docking to the SBP of the orthosteric ligand-receptor complex. This strategy was employed for a number of aminergic receptors. First, we designed dopamine D₃ preferring bitopic compounds from a D₂ selective orthosteric ligand. Next, we designed 5-HT_2B selective bitopic compounds starting from the 5-HT_1B preferring ergoline core of LSD. Comparing the serotonergic profiles of the new derivatives to that of LSD, we found that these derivatives became significantly biased towards the desired 5-HT_2B receptor target. Finally, addressing the known limitations of H₁ antihistamines, our protocol was successfully used to eliminate the well-known side effects related to the muscarinic M₁ activity of amitriptyline while preserving H₁ potency in some of the designed bitopic compounds. These applications highlight the usefulness of our new virtual screening protocol and offer a powerful strategy towards bitopic GPCR ligands with designed receptor profiles.

Structural insights into the human D1 and D2 dopamine receptor signaling complexes

The D1- and D2-dopamine receptors (D1R and D2R), which signal through Gs and Gi, respectively, represent the principal stimulatory and inhibitory dopamine receptors in the central nervous system. D1R and D2R also represent the main therapeutic targets for Parkinson's disease, schizophrenia, and many other neuropsychiatric disorders, and insight into their signaling is essential for understanding both therapeutic and side effects of dopaminergic drugs. Here, we report four cryoelectron microscopy (cryo-EM) structures of D1R-Gs and D2R-Gi signaling complexes with selective and non-selective dopamine agonists, including two currently used anti-Parkinson's disease drugs, apomorphine and bromocriptine. These structures, together with mutagenesis studies, reveal the conserved binding mode of dopamine agonists, the unique pocket topology underlying ligand selectivity, the conformational changes in receptor activation, and potential structural determinants for G protein-coupling selectivity. These results provide both a molecular understanding of dopamine signaling and multiple structural templates for drug design targeting the dopaminergic system.

Current and Future Challenges in Modern Drug Discovery

Drug discovery is an expensive, time-consuming, and risky business. To avoid late-stage failure, learnings from past projects and the development of new approaches are crucial. New modalities and emerging new target spaces allow the exploration of unprecedented indications or to address so far undrugable targets. Late-stage attrition is usually attributed to the lack of efficacy or to compound-related safety issues. Efficacy has been shown to be related to a strong genetic link to human disease, a better understanding of the target biology, and the availability of biomarkers to bridge from animals to humans. Compound safety can be improved by ligand optimization, which is becoming increasingly demanding for difficult targets. Therefore, new strategies include the design of allosteric ligands, covalent binders, and other modalities. Design methods currently heavily rely on artificial intelligence and advanced computational methods such as free energy calculations and quantum chemistry. Especially for quantum chemical methods, a more detailed overview is given in this chapter.

Current computer-aided drug design methodologies in discovery of novel drug candidates for neuropsychiatric and inflammatory diseases

Drug discovery and development is a very challenging, expensive and time-consuming process. Impressive technological advances in computer sciences and molecular biology have made it possible to use computer-aided drug design (CADD) methods in various stages of the drug discovery and development pipeline. Nowadays, CADD presents an efficacious and indispensable tool, widely used in medicinal chemistry, to lead rational drug design and synthesis of novel compounds. In this article, an overview of commonly used CADD approaches from hit identification to lead optimization was presented. Moreover, different aspects of design of multitarget ligands for neuropsychiatric and anti-inflammatory diseases were summarized. Apparently, designing multi-target directed ligands for treatment of various complex diseases may offer better efficacy, and fewer side effects. Antipsychotics that act through aminergic G protein-coupled receptors (GPCRs), especially Dopamine D2 and serotonin 5-HT2A receptors, are the best option for treatment of various symptoms associated with neuropsychiatric disorders. Furthermore, multi-target directed cyclooxygenase-2 (COX-2) and 5-lipoxygenase (5-LOX) inhibitors are also a successful approach to aid the discovery of new anti-inflammatory drugs with fewer side effects. Overall, employing CADD approaches in the process of rational drug design provides a great opportunity for future development, allowing rapid identification of compounds with the optimal polypharmacological profile.

Biased Ligands Differentially Shape the Conformation of the Extracellular Loop Region in 5-HT2B Receptors

G protein-coupled receptors are linked to various intracellular transducers, each pathway associated with different physiological effects. Biased ligands, capable of activating one pathway over another, are gaining attention for their therapeutic potential, as they could selectively activate beneficial pathways whilst avoiding those responsible for adverse effects. We performed molecular dynamics simulations with known β-arrestin-biased ligands like lysergic acid diethylamide and ergotamine in complex with the 5-HT_2B receptor and discovered that the extent of ligand bias is directly connected with the degree of closure of the extracellular loop region. Given a loose allosteric coupling of extracellular and intracellular receptor regions, we delineate a concept for biased signaling at serotonin receptors, by which conformational interference with binding pocket closure restricts the signaling repertoire of the receptor. Molecular docking studies of biased ligands gathered from the BiasDB demonstrate that larger ligands only show plausible docking poses in the ergotamine-bound structure, highlighting the conformational constraints associated with bias. This emphasizes the importance of selecting the appropriate receptor conformation on which to base virtual screening workflows in structure-based drug design of biased ligands. As this mechanism of ligand bias has also been observed for muscarinic receptors, our studies provide a general mechanism of signaling bias transferable between aminergic receptors.

The Role of Orthosteric Building Blocks of Bitopic Ligands for Muscarinic M1 Receptors

The muscarinic M₁ acetylcholine receptor is an important drug target for the treatment of various neurological disorders. Designing M₁ receptor-selective drugs has proven challenging, mainly due to the high conservation of the acetylcholine binding site among muscarinic receptor subtypes. Therefore, less conserved and topographically distinct allosteric binding sites have been explored to increase M₁ receptor selectivity. In this line, bitopic ligands, which target orthosteric and allosteric binding sites simultaneously, may provide a promising strategy. Here, we explore the allosteric, M₁-selective BQCAd scaffold derived from BQCA as a starting point for the design, synthesis, and pharmacological evaluation of a series of novel bitopic ligands in which the orthosteric moieties and linker lengths are systematically varied. Since β-arrestin recruitment seems to be favorable to therapeutic implication, all the compounds were investigated by G protein and β-arrestin assays. Some bitopic ligands are partial to full agonists for G protein activation, some activate β-arrestin recruitment, and the degree of β-arrestin recruitment varies according to the respective modification. The allosteric BQCAd scaffold controls the positioning of the orthosteric ammonium group of all ligands, suggesting that this interaction is essential for stimulating G protein activation. However, β-arrestin recruitment is not affected. The novel set of bitopic ligands may constitute a toolbox to study the requirements of β-arrestin recruitment during ligand design for therapeutic usage.

Calcium signaling mediated by aminergic GPCRs is impaired by the PI3K inhibitor LY294002 and its analog LY303511 in a PI3K-independent manner

The phosphoinositide 3-kinase (PI3K) inhibitor LY294002 (LY294) and its much less active analog LY303511 (LY303) constitute the paired probe that is commonly used to demonstrate the involvement of PI3K in intracellular signaling. We studied effects of LY294 and LY303 on Ca²⁺ signaling initiated by certain GPCR agonists in cells of several lines, including CHO cells expressing the recombinant serotonin receptor 5-HT2C and mesenchymal stromal cells derived from the human adipose tissue (AD-MSCs) and umbilical cord (UD-MSCs). The LY294/LY303 pair exerted apparently specific effects on responsiveness of AD-MSCs to ATP, suggesting the involvement of PI3K in ATP transduction. Surprisingly, LY303 inhibited Ca²⁺ transients elicited by histamine in the same cells, while LY294 was ineffective. This observation and other findings implicated a PI3K-unrelated mechanism in mediating effects of the LY compound on AD-MSC responsiveness to histamine. With LY303 in the bath, the dose dependence of histamine responses was shifted positively at the invariable number of responsive cells, as would be the case with a competitive antagonist of histamine receptors. Moreover, LY303 and LY294 inhibited Ca²⁺ transients elicited by acetylcholine and serotonin in UD-MSCs and CHO/5-HT2C cells, respectively. Our overall results argued for the possibility that LY294 and LY303 could directly affect activity of aminergic GPCRs. Thus, LY303511 and LY294002 should be used cautiously in studies of PI3K as a factor of GPCR signaling.

Design and Synthesis of Bitopic 2-Phenylcyclopropylmethylamine (PCPMA) Derivatives as Selective Dopamine D3 Receptor Ligands

2-Phenylcyclopropylmethylamine (PCPMA) analogues have been reported as selective serotonin 2C agonists. On the basis of the same scaffold, we designed and synthesized a series of bitopic derivatives as dopamine D3R ligands. A number of these new compounds show a high binding affinity for D3R with excellent selectivity. Compound (1R,2R)-22e and its enantiomer (1S,2S)-22e show a comparable binding affinity for the D3R, but the former is a potent D3R agonist, while the latter acts as an antagonist. Molecular docking studies revealed different binding poses of the PCPMA moiety within the orthosteric binding pocket of the D3R, which might explain the different functional profiles of the enantiomers. Compound (1R,2R)-30q shows a high binding affinity for the D3R (K_i = 2.2 nM) along with good selectivity, as well as good bioavailability and brain penetration properties in mice. These results reveal that the PCPMA scaffold may serve as a privileged scaffold for the design of aminergic GPCR ligands.

Binding Interactions of Ergotamine and Dihydroergotamine to 5-Hydroxytryptamine Receptor 1B (5-HT1b) Using Molecular Dynamics Simulations and Dynamic Network Analysis

Ergotamine (ERG) and dihydroergotamine (DHE), common migraine drugs, have small structural differences but lead to clinically important distinctions in their pharmacological profiles. For example, DHE is less potent than ERG by about 10-fold at the 5-hydroxytrptamine receptor 1B (5-HT_1B). Although the high-resolution crystal structures of the 5-HT_1B receptor with both ligands have been solved, the high similarity between these two complex structures does not sufficiently explain their activity differences and the activation mechanism of the receptor. Hence, an examination of the dynamic motion of both drugs with the receptor is required. In this study, we ran a total of 6.0 μs molecular dynamics simulations on each system. Our simulation data show the subtle variations between the two systems in terms of the ligand-receptor interactions and receptor secondary structures. More importantly, the ligand and protein root-mean-square fluctuations (RMSFs) for the two systems were distinct, with ERG having a trend of lower RMSF values, indicating it to be bound tighter to 5-HT_1B with less fluctuations. The molecular mechanism-general born surface area (MM-GBSA) binding energies illustrate this further, proving ERG has an overall stronger MM-GBSA binding energy. Analysis of several different microswitches has shown that the 5-HT_1B-ERG complex is in a more active conformation state than 5-HT_1B-DHE, which is further supported by the dynamic network model, with reference to mutagenesis data with the critical nodes and the first three low-energy modes from the normal mode analysis. We also identify Trp327^6.48 and Phe331^6.52 as key residues involved in the active state 5-HT_1B for both ligands. Using the detailed dynamic information from our analysis, we made predictions for possible modifications to DHE and ERG that yielded five derivatives that might have more favorable binding energies and reduced structural fluctuations.