Structure versus function—The impact of computational methods on the discovery of specific GPCR–ligands

Structure and dynamics of the somatostatin receptor 3-ligand binding in the presence of lipids examined using computational structural biology methods

In the past two decades, the structural biology studies on G-protein coupled receptors (GPCRs) are on the rise. Understanding the relation between the structure and function of GPCRs is important as they play a huge role in various signaling mechanisms in a eukaryotic cell. Somatostatin receptor 3 (SSTR3), one of the GPCRs, is one such important receptor which oversees different cellular processes including cell-to-cell signaling. However, the information available regarding the structural features of SSTR3 responsible for their bioactivity is scarce. In this study, we report a structural understanding of SSTR3-ligand binding that could be helpful in demystifying the structural complexities related to functioning of the receptor. An integrated protocol consisting of different computational structural biology tools including protein structure prediction via comparative modeling, binding site characterization, three-dimensional quantitative structure-activity relationship based on comparative molecular field analysis and comparative molecular similarity indices analysis, density functional theory, and molecular dynamics simulations were performed. Different understandings from the simulation of SSTR3-ligand complexes, mainly the conditions that are favorable for the formation of lowest bioactive state of SSTR3 ligands are reported. In addition to that, we report the important physicochemical descriptors of SSTR3 ligands that could significantly influence their bioactivity. The results of the study could be helpful in developing novel SSTR3 ligands (both agonists and antagonists) with high potency and receptor selectivity.

Bitter taste in silico: A review on virtual ligand screening and characterization methods for TAS2R-bitterant interactions

The growing pharmaceutical interest in the human bitter taste receptors (hTAS2Rs) has two dimensions; i) evaluation of the bitterness of active pharmaceutical compounds, in order to develop strategies for improving patients' adherence to medication, and ii) application of ligands for extra-cellular hTAS2Rs for potential preventive therapeutic achievements. The result is an increasing demand on robust tools for bitterness assessment and screening the receptor-ligand affinity. In silico tools are useful for aiding experimental-screening, as well as to elucide ligand-receptor interactions. In this review, the ligand-based and structure-based approaches are described as the two main in silico tools for bitter taste analysis. The strengths and weaknesses of each approach are discussed. Both approaches provide key tools for understanding and exploiting bitter taste for human health applications.

Molecular-Level Understanding of the Somatostatin Receptor 1 (SSTR1)-Ligand Binding: A Structural Biology Study Based on Computational Methods

Somatostatin receptor 1 (SSTR1), a subtype of somatostatin receptors, is involved in various signaling mechanisms in different parts of the human body. Like most of the G-protein-coupled receptors (GPCRs), the available information on the structural features of SSTR1 responsible for the biological activity is scarce. In this study, we report a molecular-level understanding of SSTR1-ligand binding, which could be helpful in solving the structural complexities involved in SSTR1 functioning. Based on a three-dimensional quantitative structure-activity relationship (3D-QSAR) study using comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA), we have identified that an electronegative, less-bulkier, and hydrophobic atom substitution can substantially increase the biological activity of SSTR1 ligands. A density functional theory (DFT) study has been followed to study the electron-related properties of the SSTR1 ligands and to validate the results obtained via the 3D-QSAR study. 3D models of SSTR1-ligand systems have been embedded in lipid-lipid bilayer membranes to perform molecular dynamics (MD) simulations. Analysis of the MD trajectories reveals important information about the crucial residues involved in SSTR1-ligand binding and various conformational changes in the protein that occur after ligand binding. Additionally, we have identified the probable ligand-binding site of SSTR1 and validated it using MD. We have also studied the favorable conditions that are essential for forming the most stable and lowest-energy bioactive conformation of the ligands inside the binding site. The results of the study could be useful in constructing more potent and novel SSTR1 antagonists and agonists.

Development of the first in vivo GPR17 ligand through an iterative drug discovery pipeline: A novel disease-modifying strategy for multiple sclerosis

The GPR17 receptor, expressed on oligodendroglial precursors (OPCs, the myelin producing cells), has emerged as an attractive target for a pro-myelinating strategy in multiple sclerosis (MS). However, the proof-of-concept that selective GPR17 ligands actually exert protective activity in vivo is still missing. Here, we exploited an iterative drug discovery pipeline to prioritize novel and selective GPR17 pro-myelinating agents out of more than 1,000,000 compounds. We first performed an in silico high-throughput screening on GPR17 structural model to identify three chemically-diverse ligand families that were then combinatorially exploded and refined. Top-scoring compounds were sequentially tested on reference pharmacological in vitro assays with increasing complexity, ending with myelinating OPC-neuron co-cultures. Successful ligands were filtered through in silico simulations of metabolism and pharmacokinetics, to select the most promising hits, whose dose and ability to target the central nervous system were then determined in vivo. Finally, we show that, when administered according to a preventive protocol, one of them (named by us as galinex) is able to significantly delay the onset of experimental autoimmune encephalomyelitis (EAE), a mouse model of MS. This outcome validates the predictivity of our pipeline to identify novel MS-modifying agents.

Current Advances in Allosteric Modulation of Muscarinic Receptors

Allosteric modulators are ligands that bind to a site on the receptor that is spatially separated from the orthosteric binding site for the endogenous neurotransmitter. Allosteric modulators modulate the binding affinity, potency, and efficacy of orthosteric ligands. Muscarinic acetylcholine receptors are prototypical allosterically-modulated G-protein-coupled receptors. They are a potential therapeutic target for the treatment of psychiatric, neurologic, and internal diseases like schizophrenia, Alzheimer’s disease, Huntington disease, type 2 diabetes, or chronic pulmonary obstruction. Here, we reviewed the progress made during the last decade in our understanding of their mechanisms of binding, allosteric modulation, and in vivo actions in order to understand the translational impact of studying this important class of pharmacological agents. We overviewed newly developed allosteric modulators of muscarinic receptors as well as new spin-off ideas like bitopic ligands combining allosteric and orthosteric moieties and photo-switchable ligands based on bitopic agents.

Strategies for the discovery of biased GPCR ligands

G-protein-coupled receptors (GPCRs) represent important drug targets with complex pharmacological characteristics. Biased signaling represents one important dimension, describing ligand-dependent shifts of naturally imprinted signaling profiles. Because biased GPCR modulators provide potential therapeutic benefits including higher efficiencies and reduced adverse effects, the identification of such ligands as drug candidates is highly desirable. This review aims to provide an overview of the challenges and strategies in the discovery of biased ligands. We show different approaches for biased ligand discovery in the example of G-protein-biased opioid analgesics and discuss possibilities to design biased ligands by targeting extracellular receptor regions.

In Silico and Structural Analyses Demonstrate That Intrinsic Protein Motions Guide T Cell Receptor Complementarity Determining Region Loop Flexibility

T-cell immunity is controlled by T cell receptor (TCR) binding to peptide major histocompatibility complexes (pMHCs). The nature of the interaction between these two proteins has been the subject of many investigations because of its central role in immunity against pathogens, cancer, in autoimmunity, and during organ transplant rejection. Crystal structures comparing unbound and pMHC-bound TCRs have revealed flexibility at the interaction interface, particularly from the perspective of the TCR. However, crystal structures represent only a snapshot of protein conformation that could be influenced through biologically irrelevant crystal lattice contacts and other factors. Here, we solved the structures of three unbound TCRs from multiple crystals. Superposition of identical TCR structures from different crystals revealed some conformation differences of up to 5 Å in individual complementarity determining region (CDR) loops that are similar to those that have previously been attributed to antigen engagement. We then used a combination of rigidity analysis and simulations of protein motion to reveal the theoretical potential of TCR CDR loop flexibility in unbound state. These simulations of protein motion support the notion that crystal structures may only offer an artifactual indication of TCR flexibility, influenced by crystallization conditions and crystal packing that is inconsistent with the theoretical potential of intrinsic TCR motions.

Modeling and Deorphanization of Orphan GPCRs

Despite tremendous efforts, approximately 120 GPCRs remain orphan. Their physiological functions and their potential roles in diseases are poorly understood. Orphan GPCRs are extremely important because they may provide novel therapeutic targets for unmet medical needs. As a complement to experimental approaches, molecular modeling and virtual screening are efficient techniques to discover synthetic surrogate ligands which can help to elucidate the role of oGPCRs. Constitutively activated mutants and recently published active structures of GPCRs provide stimulating opportunities for building active molecular models for oGPCRs and identifying activators using virtual screening of compound libraries. We describe the molecular modeling and virtual screening process we have applied in the discovery of surrogate ligands, and provide examples for CCKA, a simulated oGPCR, and for two oGPCRs, GPR52 and GPR34.

Ligand binding modes from low resolution GPCR models and mutagenesis: chicken bitter taste receptor as a test-case

Bitter taste is one of the basic taste modalities, warning against consuming potential poisons. Bitter compounds activate members of the bitter taste receptor (Tas2r) subfamily of G protein-coupled receptors (GPCRs). The number of functional Tas2rs is species-dependent. Chickens represent an intriguing minimalistic model, because they detect the bitter taste of structurally different molecules with merely three bitter taste receptor subtypes. We investigated the binding modes of several known agonists of a representative chicken bitter taste receptor, ggTas2r1. Because of low sequence similarity between ggTas2r1 and crystallized GPCRs (~10% identity, ~30% similarity at most), the combination of computational approaches with site-directed mutagenesis was used to characterize the agonist-bound conformation of ggTas2r1 binding site between TMs 3, 5, 6 and 7. We found that the ligand interactions with N93 in TM3 and/or N247 in TM5, combined with hydrophobic contacts, are typically involved in agonist recognition. Next, the ggTas2r1 structural model was successfully used to identify three quinine analogues (epiquinidine, ethylhydrocupreine, quinidine) as new ggTas2r1 agonists. The integrated approach validated here may be applicable to additional cases where the sequence identity of the GPCR of interest and the existing experimental structures is low.

Ligand-Specific Restriction of Extracellular Conformational Dynamics Constrains Signaling of the M2 Muscarinic Receptor

G protein-coupled receptors transmit extracellular signals across cell membranes via different G protein classes and β-arrestins. Some pathways may be therapeutically beneficial, whereas others may be detrimental under certain pathophysiological conditions. For many GPCRs, biased agonists are available, which preferentially signal through one pathway or a subset of pathways, and harnessing biased agonism could be a potential novel therapeutic strategy. However, the incomplete mechanistic understanding of biased agonism hampers rational design of biased ligands. Using the muscarinic M₂ receptor as a model system, we have analyzed the relationship between ligand-dependent conformational changes as revealed in all-atom MD simulations and the activation of specific G proteins. We find that the extent of closure of the extracellular, allosteric binding site interferes with the activation of certain G proteins. Our data allow the rational design of G_i-biased agonists at the M₂ receptor and delineate a simple principle which may be translated to other GPRCs.

Exploring the binding properties and structural stability of an opsin in the chytrid Spizellomyces punctatus using comparative and molecular modeling

Background

Opsin proteins are seven transmembrane receptor proteins which detect light. Opsins can be classified into two types and share little sequence identity: type 1, typically found in bacteria, and type 2, primarily characterized in metazoa. The type 2 opsins (Rhodopsins) are a subfamily of G-protein coupled receptors (GPCRs), a large and diverse class of seven transmembrane proteins and are generally restricted to metazoan lineages. Fungi use light receptors including opsins to sense the environment and transduce signals for developmental or metabolic changes. Opsins characterized in the Dikarya (Ascomycetes and Basidiomycetes) are of the type 1 bacteriorhodopsin family but the early diverging fungal lineages have not been as well surveyed. We identified by sequence similarity a rhodopsin-like GPCR in genomes of early diverging chytrids and examined the structural characteristics of this protein to assess its likelihood to be homologous to animal rhodopsins and bind similar chromophores.

Methods

We used template-based structure modeling, automated ligand docking, and molecular modeling to assess the structural and binding properties of an identified opsin-like protein found in Spizellomyces punctatus, a unicellular, flagellated species belonging to Chytridiomycota, one of the earliest diverging fungal lineages. We tested if the sequence and inferred structure were consistent with a solved crystal structure of a type 2 rhodopsin from the squid Todarodes pacificus.

Results

Our results indicate that the Spizellomyces opsin has structural characteristics consistent with functional animal type 2 rhodopsins and is capable of maintaining a stable structure when associated with the retinaldehyde chromophore, specifically the 9-cis-retinal isomer. Together, these results support further the homology of Spizellomyces opsins to animal type 2 rhodopsins.

Discussion

This represents the first test of structure/function relationship of a type 2 rhodopsin identified in early branching fungal lineages, and provides a foundation for future work exploring pathways and components of photoreception in early fungi.

Structural determinants of diphenethylamines for interaction with the κ opioid receptor: Synthesis, pharmacology and molecular modeling studies

The κ opioid (KOP) receptor crystal structure in an inactive state offers nowadays a valuable platform for inquiry into receptor function. We describe the synthesis, pharmacological evaluation and docking calculations of KOP receptor ligands from the class of diphenethylamines using an active-like structure of the KOP receptor attained by molecular dynamics simulations. The structure-activity relationships derived from computational studies was in accordance with pharmacological activities of targeted diphenethylamines at the KOP receptor established by competition binding and G protein activation in vitro assays. Our analysis identified that agonist binding results in breaking of the Arg156-Thr273 hydrogen bond, which stabilizes the inactive receptor conformation, and a crucial hydrogen bond with His291 is formed. Compounds with a phenolic 4-hydroxy group do not form the hydrogen bond with His291, an important residue for KOP affinity and agonist activity. The size of the N-substituent hosted by the hydrophobic pocket formed by Val108, Ile316 and Tyr320 considerably influences binding and selectivity, with the n-alkyl size limit being five carbon atoms, while bulky substituents turn KOP agonists in antagonists. Thus, combination of experimental and molecular modeling strategies provides an initial framework for understanding the structural features of diphenethylamines that are essential to promote binding affinity and selectivity for the KOP receptor, and may be involved in transduction of the ligand binding event into molecular changes, ultimately leading to receptor activation.

Ligand Binding Ensembles Determine Graded Agonist Efficacies at a G Protein-coupled Receptor

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

More than a look into a crystal ball: protein structure elucidation guided by molecular dynamics simulations

The 'form follows function' principle implies that a structural determination of protein structures is indispensable to understand proteins in their biological roles. However, experimental methods still show shortcomings in the description of the dynamic properties of proteins. Therefore, molecular dynamics (MD) simulations represent an essential tool for structural biology to investigate proteins as flexible and dynamic entities. Here, we will give an overview on the impact of MD simulations on structural investigations, including studies that aim at a prediction of protein-folding pathways, protein-assembly processes and the sampling of conformational space by computational means.

In silico Exploration of the Conformational Universe of GPCRs

The structural plasticity of G protein coupled receptors (GPCRs) leads to a conformational universe going from inactive to active receptor states with several intermediate states. Many of them have not been captured yet and their role for GPCR activation is not well understood. The study of this conformational space and the transition dynamics between different receptor populations is a major challenge in molecular biophysics. The rational design of effector molecules that target such receptor populations allows fine-tuning receptor signalling with higher specificity to produce drugs with safer therapeutic profiles. In this minireview, we outline highly conserved receptor regions which are considered determinant for the establishment of distinct receptor states. We then discuss in-silico approaches such as dimensionality reduction methods and Markov State Models to explore the GPCR conformational universe and exploit the obtained conformations through structure-based drug design.

Beyond Membrane Protein Structure: Drug Discovery, Dynamics and Difficulties

Most of the previous content of this book has focused on obtaining the structures of membrane proteins. In this chapter we explore how those structures can be further used in two key ways. The first is their use in structure based drug design (SBDD) and the second is how they can be used to extend our understanding of their functional activity via the use of molecular dynamics. Both aspects now heavily rely on computations. This area is vast, and alas, too large to consider in depth in a single book chapter. Thus where appropriate we have referred the reader to recent reviews for deeper assessment of the field. We discuss progress via the use of examples from two main drug target areas; G-protein coupled receptors (GPCRs) and ion channels. We end with a discussion of some of the main challenges in the area.