No abstract

3D pharmacophore models are three‐dimensional ensembles of chemically defined interactions of a ligand in its bioactive conformation. They represent an elegant way to decipher chemically encoded ligand information and have therefore become a valuable tool in drug design. In this review, we provide an overview on the basic concept of this method and summarize key studies for applying 3D pharmacophore models in virtual screening and mechanistic studies for protein functionality. Moreover, we discuss recent developments in the field. The combination of 3D pharmacophore models with molecular dynamics simulations could be a quantum leap forward since these approaches consider macromolecule–ligand interactions as dynamic and therefore show a physiologically relevant interaction pattern. Other trends include the efficient usage of 3D pharmacophore information in machine learning and artificial intelligence applications or freely accessible web servers for 3D pharmacophore modeling. The recent developments show that 3D pharmacophore modeling is a vibrant field with various applications in drug discovery and beyond. This article is categorized under: Computer and Information Science > Chemoinformatics Computer and Information Science > Computer Algorithms and Programming Molecular and Statistical Mechanics > Molecular Interactions

GPCRs modulate a plethora of physiological processes and mediate the effects of one-third of FDA-approved drugs. Depending on which ligand activates a receptor, it can engage different intracellular transducers. This 'biased signalling' paradigm requires that we now characterize physiological signalling not just by receptors but by ligand-receptor pairs. Ligands eliciting biased signalling may constitute better drugs with higher efficacy and fewer adverse effects. However, ligand bias is very complex, making reproducibility and description challenging. Here, we provide guidelines and terminology for any scientists to design and report ligand bias experiments.The guidelines will aid consistency and clarity, as the basic receptor research and drug discovery communities continue to advance our understanding and exploitation of ligand bias. Scientific insight, biosensors, and analytical methods are still evolving and should benefit from and contribute to the implementation of the guidelines, together improving translation from in vitro to disease-relevant in vivo models. | INTRODUCTIONThe $800 human GPCRs transduce sensory inputs and systemic signals into appropriate cellular responses in numerous physiological processes. They recognize a vast diversity of signals ranging from photons, tastants and odours to ions, neurotransmitters, hormones, and cytokines (Harding et al., 2021;Wacker, Stevens, & Roth, 2017).Even though GPCRs represent the primary target of 34% of FDAapproved drugs, more than 220 non-olfactory GPCRs have disease associations which are as yet untapped in clinical research (Hauser, Attwood, Rask-Andersen, Schioth, & Gloriam, 2017;Sriram & Insel, 2018). Despite the diversity of extracellular ligands and physiological roles of GPCRs, these cell surface receptors share a conserved molecular fold and intracellular transducers. Agonist binding stabilizes active conformations of the receptor, facilitating the binding of one or more cytosolic transducer proteins. These include the heterotrimeric G proteins consisting of α, β and γ subunits that dissociate to α and βγ upon activation by the receptor. G proteins comprise 16 distinct α subunits and are divided into four families based on homology and associated downstream signalling pathways: G s (G s and

G protein-coupled receptors constitute the largest family of membrane receptors and modulate almost every physiological process in humans. Binding of agonists to G protein-coupled receptors induces a shift from inactive to active receptor conformations. Biophysical studies of the dynamic equilibrium of receptors suggest that a portion of receptors can remain in inactive states even in the presence of saturating concentrations of agonist and G protein mimetic. However, the molecular details of agonist-bound inactive receptors are poorly understood.Here we use the model of bitopic orthosteric/allosteric (i.e. dualsteric) agonists for muscarinic M 2 receptors to demonstrate the existence and function of such inactive agonist⅐receptor complexes on a molecular level. Using all-atom molecular dynamics simulations, dynophores (i.e. a combination of static three-dimensional pharmacophores and molecular dynamics-based conformational sampling), ligand design, and receptor mutagenesis, we show that inactive agonist⅐receptor complexes can result from agonist binding to the allosteric vestibule alone, whereas the dualsteric binding mode produces active receptors. Each agonist forms a distinct ligand binding ensemble, and different agonist efficacies depend on the fraction of purely allosteric (i.e. inactive) versus dualsteric (i.e. active) binding modes. We propose that this concept may explain why agonist⅐receptor complexes can be inactive and that adopting multiple binding modes may be generalized also to small agonists where binding modes will be only subtly different and confined to only one binding site.Specific and coordinated cell-to-cell communication regulates the flow of information between cells, and proper information processing ensures physiological functions of biological systems. G protein-coupled receptors (GPCRs), 8 constituting the largest class of membrane proteins in mammals, are essential mediators of chemically and light-encoded information (1-4). GPCRs sense a great variety of extracellular stimuli, e.g. neurotransmitters and hormones, and subsequently translate this information into an intracellular response via G proteins, ␤-arrestins, and possibly GPCR-interacting proteins (2-5). Because of their abundance and relevance in regulating the majority of (patho-)physiological processes in humans, GPCRs have for a long time represented the most important drug targets being addressed by at least a third of all currently marketed drugs (6, 7).Agonist binding leads to receptor activation, which is followed by intracellular G protein recruitment and subsequent cell signaling. Breakthroughs in GPCR crystallography have led to inactive and active crystal structures of the same receptor protein. Among these are rhodopsin (8 -10) and more recently the ␤ 2 -adrenergic (11-14), M 2 muscarinic (15, 16), and -opioid receptors (17, 18). These structures most likely represent energetically favored, relatively stable inactive and active receptor⅐ligand complexes. Despite the diversity of crystallized receptors, a common...

Small molecules interfering with Rac1 activation are considered as potential drugs and are already studied in animal models. A widely used inhibitor without reported attenuation of RhoA activity is NSC23766 [(. We found that NSC23766 inhibits the M 2 muscarinic acetylcholine receptor (M 2 mAChR)-induced Rac1 activation in neonatal rat cardiac myocytes. Surprisingly, NSC27366 concomitantly suppressed the carbachol-induced RhoA activation and a M 2 mAChR-induced inotropic response in isolated neonatal rat hearts requiring the activation of Rhodependent kinases. We therefore aimed to identify the mechanisms by which NSC23766 interferes with the differentially mediated, M 2 mAChR-induced responses. Interestingly, NSC23766 caused a rightward shift of the carbachol concentration response curve for the positive inotropic response without modifying carbachol efficacy. To analyze the specificity of NSC23766, we compared the carbachol and the similarly G i bg-mediated, adenosine-induced activation of G i protein-regulated potassium channel (GIRK) channels in human atrial myocytes. Application of NSC23766 blocked the carbachol-induced K 1 current but had no effect on the adenosine-induced GIRK current. Similarly, an adenosine A 1 receptor-induced positive inotropic response in neonatal rat hearts was not attenuated by NSC23766. To investigate its specificity toward the different mAChR types, we studied the carbachol-induced elevation of intracellular Ca 21 concentrations in human embryonic kidney 293 (HEK-293) cells expressing M 1 , M 2 , or M 3 mAChRs. NSC23766 caused a concentration-dependent rightward shift of the carbachol concentration response curves at all mAChRs. Thus, NSC23766 is not only an inhibitor of Rac1 activation, but it is within the same concentration range a competitive antagonist at mAChRs. Molecular docking analysis at M 2 and M 3 mAChR crystal structures confirmed this interpretation.