An ontology for pharmaceutical ligands and its application for in silico screening and library design

Annotation efforts in biosciences have focused in past years mainly on the annotation of genomic sequences. Only very limited effort has been put into annotation schemes for pharmaceutical ligands. Here we propose annotation schemes for the ligands of four major target classes, enzymes, G protein-coupled receptors (GPCRs), nuclear receptors (NRs), and ligand-gated ion channels (LGICs), and outline their usage for in silico screening and combinatorial library design. The proposed schemes cover ligand functionality and hierarchical levels of target classification. The classification schemes are based on those established by the EC, GPCRDB, NuclearDB, and LGICDB. The ligands of the MDL Drug Data Report (MDDR) database serve as a reference data set of known pharmacologically active compounds. All ligands were annotated according to the schemes when attribution was possible based on the activity classification provided by the reference database. The purpose of the ligand-target classification schemes is to allow annotation-based searching of the ligand database. In addition, the biological sequence information of the target is directly linkable to the ligand, hereby allowing sequence similarity-based identification of ligands of next homologous receptors. Ligands of specified levels can easily be retrieved to serve as comprehensive reference sets for cheminformatics-based similarity searches and for design of target class focused compound libraries. Retrospective in silico screening experiments within the MDDR01.1 database, searching for structures binding to dopamine D2, all dopamine receptors and all amine-binding class A GPCRs using known dopamine D2 binding compounds as a reference set, have shown that such reference sets are in particular useful for the identification of ligands binding to receptors closely related to the reference system. The potential for ligand identification drops with increasing phylogenetic distance. The analysis of the focus of a tertiary amine based combinatorial library compared to known amine binding class A GPCRs, peptide binding class A GPCRs, and LGIC ligands constitutes a second application scenario which illustrates how the focus of a combinatorial library can be treated quantitatively. The provided annotation schemes, which bridge chem- and bioinformatics by linking ligands to sequences, are expected to be of key utility for further systematic chemogenomics exploration of previously well explored target families.     

Virtual screen for ligands of orphan G protein-coupled receptors

This paper describes a virtual screening methodology that generates a ranked list of high-binding small molecule ligands for orphan G protein-coupled receptors (oGPCRs), circumventing the requirement for receptor three-dimensional structure determination. Features representing the receptor are based only on physicochemical properties of primary amino acid sequence, and ligand features use the two-dimensional atomic connection topology and atomic properties. An experimental screen comprised nearly 2 million hypothetical oGPCR-ligand complexes, from which it was observed that the top 1.96% predicted affinity scores corresponded to "highly active" ligands against orphan receptors. Results representing predicted high-scoring novel ligands for many oGPCRs are presented here. Validation of the method was carried out in several ways: (1) A random permutation of the structure-activity relationship of the training data was carried out; by comparing test statistic values of the randomized and nonshuffled data, we conclude that the value obtained with nonshuffled data is unlikely to have been encountered by chance. (2) Biological activities linked to the compounds with high cross-target binding affinity were analyzed using computed log-odds from a structure-based program. This information was correlated with literature citations where GPCR-related pathways or processes were linked to the bioactivity in question. (3) Anecdotal, out-of-sample predictions for nicotinic targets and known ligands were performed, with good accuracy in the low-to-high "active" binding range. (4) An out-of-sample consistency check using the commercial antipsychotic drug olanzapine produced "active" to "highly-active" predicted affinities for all oGPCRs in our study, an observation that is consistent with documented findings of cross-target affinity of this compound for many different GPCRs. It is suggested that this virtual screening approach may be used in support of the functional characterization of oGPCRs by identifying potential cognate ligands. Ultimately, this approach may have implications for pharmaceutical therapies to modulate the activity of faulty or disease-related cellular signaling pathways. In addition to application to cell surface receptors, this approach is a generalized strategy for discovery of small molecules that may bind intracellular enzymes and involve protein-protein interactions.     

Structural basis for thrombin activation of a protease-activated receptor: inhibition of intramolecular liganding

Protease-activated G protein-coupled receptors (PAR1-4) are tethered-ligand receptors that are activated by proteolytic cleavage of the extracellular domain (exodomain) of the receptor. PAR1, the prototypic member of the PAR family, is the high-affinity thrombin receptor of platelets and vascular endothelium and plays a critical role in blood coagulation, thrombosis, and inflammation. Here, we describe the solution structure of the thrombin-cleaved exodomain of PAR1. The side chains of a hydrophobic hirudin-like (Hir) sequence and adjacent anionic motif project into solution. Docking of the exodomain Hir sequence to exosite I of thrombin reveals that the tethered ligand in the cleaved exodomain bends away from thrombin, leaving its active site available to another large macromolecular substrate. The N-terminal ligand is longer than anticipated and forms an intramolecular complex with a region located in the C terminus of the exodomain. Mutational analysis confirmed that this C-terminal region is a ligand binding site for both intra- and intermolecular ligands. A lipidated-ligand binding site peptide was found to be an effective inhibitor of thrombin-induced platelet aggregation.     

Modeling and active site refinement for G protein-coupled receptors: application to the β-2 adrenergic receptor

It is well known that G protein-coupled receptors are prime targets for drug discovery. At the present time there is only one protein from this class that has an X-ray crystal structure, bovine rhodopsin. Crystal structures of rhodopsin have become invaluable templates for the modeling of class-A G protein-coupled receptors as they likely represent the overall topology of this family of proteins. However, because of low sequence homology within the class and the inherent mobility of integral membrane proteins, it is unlikely that this single structural template reflects the ensemble of conformations accessible for any given receptor. We have devised a procedure based upon comparative modeling that uses induced fit modeling coupled with binding site expansion. The modeling protocol enables an ensemble approach to binding mode prediction. The utility of models for β-2 adrenergic receptor will be discussed.     

Chemogenomic discovery of allosteric antagonists at the GPRC6A receptor

GPRC6A is a Family C G protein-coupled receptor recently discovered and deorphanized by our group. This study integrates chemogenomic ligand inference, homology modeling, compound synthesis, and pharmacological mechanism-of-action studies to disclose two noticeable results of methodological and pharmacological character: (1) chemogenomic lead identification through the first, to our knowledge, ligand inference between two different GPCR families, Families A and C; and (2) the discovery of the most selective GPRC6A allosteric antagonists discovered to date. The unprecedented inference of?pharmacological activity across GPCR families provides proof-of-concept for in?silico approaches against Family C targets based on Family A templates, greatly expanding the prospects of successful drug design and discovery. The antagonists were tested against a panel of seven Family A and C G protein-coupled receptors containing the chemogenomic binding sequence motif where some of the identified GPRC6A antagonists showed some activity. However, three compounds with at least ～3-fold selectivity for GPRC6A were discovered, which present a significant step forward compared with the previously published GPRC6A antagonists, calindol and NPS 2143, which both display ～30-fold selectivity for the calcium-sensing receptor compared to GPRC6A. The antagonists constitute novel research tools toward investigating the signaling mechanism of the GPRC6A receptor at the cellular level and serve as initial ligands for further optimization of potency and selectivity enabling future ex?vivo/in?vivo pharmacological studies.     

Biophysical approaches to G protein-coupled receptors: Structure, function and dynamics

G protein-coupled receptors (GPCR) represent a large family of drug targets for which there is no high-resolution structural information. In order to understand the mechanisms of ligand recognition and receptor activation, there is a strong need for novel biophysical methods. In this Perspective we provide an overview of recent experimental approaches used to explore the molecular architecture and dynamics of GPCR and their interactions with ligands and G proteins using biophysical, non-crystallographic, methods.     

The Role of Lipids in Allosteric Modulation of Dopamine D2 Receptor—In Silico Study

The dopamine D₂ receptor, belonging to the class A G protein-coupled receptors (GPCRs), is an important drug target for several diseases, including schizophrenia and Parkinson’s disease. The D₂ receptor can be activated by the natural neurotransmitter dopamine or by synthetic ligands, which in both cases leads to the receptor coupling with a G protein. In addition to receptor modulation by orthosteric or allosteric ligands, it has been shown that lipids may affect the behaviour of membrane proteins. We constructed a model of a D₂ receptor with a long intracellular loop (ICL3) coupled with G_iα1 or G_iα2 proteins, embedded in a complex asymmetric membrane, and simulated it in complex with positive, negative or neutral allosteric ligands. In this study, we focused on the influence of ligand binding and G protein coupling on the membrane–receptor interactions. We show that there is a noticeable interplay between the cell membrane, G proteins, D₂ receptor and its modulators.     

New approach to molecular docking and its application to virtual screening of chemical databases

This paper describes the validation of a molecular docking method and its application to virtual database screening. The code flexibly docks ligand molecules into rigid receptor structures using a tabu search methodology driven by an empirically derived function for estimating the binding affinity of a protein-ligand complex. The docking method has been tested on 70 ligand-receptor complexes for which the experimental binding affinity and binding geometry are known. The lowest energy geometry produced by the docking protocol is within 2.0 A root mean square of the experimental binding mode for 79% of the complexes. The method has been applied to the problem of virtual database screening to identify known ligands for thrombin, factor Xa, and the estrogen receptor. A database of 10,000 randomly chosen "druglike" molecules has been docked into the three receptor structures. In each case known receptor ligands were included in the study. The results showed good separation between the predicted binding affinities of the known ligand set and the database subset.     

Structural Insights into Ligand—Receptor Interactions Involved in Biased Agonism of G-Protein Coupled Receptors

G protein-coupled receptors (GPCRs) are versatile signaling proteins that mediate complex cellular responses to hormones and neurotransmitters. Ligand directed signaling is observed when agonists, upon binding to the same receptor, trigger significantly different configuration of intracellular events. The current work reviews the structurally defined ligand – receptor interactions that can be related to specific molecular mechanisms of ligand directed signaling across different receptors belonging to class A of GPCRs. Recent advances in GPCR structural biology allow for mapping receptors’ binding sites with residues particularly important in recognition of ligands’ structural features that are responsible for biased signaling. Various studies show particular role of specific residues lining the extended ligand binding domains, biased agonists may alternatively affect their interhelical interactions and flexibility what can be translated into intracellular loop rearrangements. Studies on opioid and angiotensin receptors indicate importance of residues located deeper within the binding cavity and direct interactions with receptor residues linking the ortosteric ligand binding site with the intracellular transducer binding domain. Collection of results across different receptors may suggest elements of common molecular mechanisms which are responsible for passing alternative signals from biased agonists.     

Influencing receptor-ligand binding mechanisms with multivalent ligand architecture

Multivalent ligands can function as inhibitors or effectors of biological processes. Potent inhibitory activity can arise from the high functional affinities of multivalent ligand-receptor interactions. Effector functions, however, are influenced not only by apparent affinities but also by alternate factors, including the ability of a ligand to cluster receptors. Little is known about the molecular features of a multivalent ligand that determine whether it will function as an inhibitor or effector. We envisioned that, by altering multivalent ligand architecture, ligands with preferences for different binding mechanisms would be generated. To this end, a series of 28 ligands possessing structural diversity was synthesized. This series provides the means to explore the effects of ligand architecture on the inhibition and clustering of a model protein, the lectin concanavalin A (Con A). The structural parameters that were varied include scaffold shape, size, valency, and density of binding elements. We found that ligands with certain architectures are effective inhibitors, but others mediate receptor clustering. Specifically, high molecular weight, polydisperse polyvalent ligands are effective inhibitors of Con A binding, whereas linear oligomeric ligands generated by the ring-opening metathesis polymerization have structural properties that favor clustering. The shape of a multivalent ligand also influences specific aspects of receptor clustering. These include the rate at which the receptor is clustered, the number of receptors in the clusters, and the average interreceptor distance. Our results indicate that the architecture of a multivalent ligand is a key parameter in determining its activity as an inhibitor or effector. Diversity-oriented syntheses of multivalent ligands coupled with effective assays that can be used to compare the contributions of different binding parameters may afford ligands that function by specific mechanisms.     

Molecular modeling of A1 and A2A adenosine receptors: comparison of rhodopsin- and beta2-adrenergic-based homology models through the docking studies

Adenosine receptors (ARs) are members of the superfamily of G protein-coupled receptors. The homology models of adenosine A1 and A2A receptors were constructed. The high-resolution X-ray structure of bovine rhodopsin and crystal structure of beta2-adrenergic receptor were used as templates. The binding sites of the A1 and A2A ARs were constructed by using data obtained from mutagenesis experiments as well as docking simulations of the respective AR antagonsists DPCPX and XAC. To compare rhodopsin- and beta2-adrenergic-based models, the binding mode of A1 (KW-3902, LUF-5437) and A2A (KW-6002, ZM-241385) ARs antagonists were also examined. The differences in the binding ability of both models were noted during the study. The beta2-adrenergic-based A2A AR model was much more capable to stabilize the ligand in the binding site cavity than the corresponding rhodopsin-based A2A AR model, however, such differences were not so clear in case of A1 AR models. It was suggested that for the A1 AR it is possible to use the crystal structure of rhodopsin as a template as well as beta2-adrenergic receptor, but for A2A AR, with the now available beta2-adrenergic receptor X-ray structure, docking studies should be avoided on the rhodopsin-based model. However, taking into account that the beta2AR shares about 31% of the residues with the AR in comparison to 21% in case of bRho, we suggest using beta2-adrenergic-based models for the A1 and A2A ARs for further in silico ligand screening also because of their generally better ability to stabilize ligands inside the binding pocket.     

Monitoring the Reversibility of GPCR Signaling by Combining Photochromic Ligands with Label-free Impedance Analysis

G protein-coupled cell surface receptors (GPCR) trigger complex intracellular signaling cascades upon agonist binding. Classic pharmacological assays provide information about binding affinities, activation or blockade at different stages of the signaling cascade, but real time dynamics and reversibility of these processes remain often disguised. We show that combining photochromic NPY receptor ligands, which can be toggled in their receptor activation ability by irradiation with light of different wavelengths, with whole cell label-free impedance assays allows observing the cell response to receptor activation and its reversibility over time. The concept demonstrated on NPY receptors may be well applicable to many other GPCRs providing a deeper insight into the time course of intracellular signaling processes.     

Discovery of novel regulatory peptides by reverse pharmacology: spotlight on chemerin and the RF-amide peptides metastin and QRFP

Reverse pharmacology is a screening technology that matches G protein-coupled receptors (GPCRs) with unknown cognate ligands in cell-based screening assays by detection of agonist-induced signaling pathways. One decade spent pursuing orphan GPCR screening by this technique assigned over 30 ligand/receptor pairs and revealed previously known or novel undescribed ligands, mostly of a peptidic nature. In this review, we describe the discovery, characterization of the structural composition, biological function, physiological role and therapeutic potential of three recently identified peptidic ligands. These are metastin, QRFP in a context of five RF-amide genes described in humans and the chemoattractant, chemerin. Metastin was initially characterized as a metastasis inhibitor. Investigations using ligand/receptor pairing revealed that metastin was involved in a variety of physiological processes, including endocrine function during pregnancy and gonad development. The novel RF-amide QRFP is implicated in food intake and aldosterone release from the adrenal cortex in the rat. Chemerin, first described as TIG2, is upregulated in tazarotene-treated psoriatic skin. By GPCR screening, bioactive chemerin was isolated from ovarial carcinoma fluid as well as hemofiltrate. Characterization as a chemoattractant for immature dendritic cells and analysis of the expression profile of metastin and its receptor suggested a physiological role of chemerin as a mediator of the immune response, inflammatory processes and bone development.     

INPHARMA-Based Determination of Ligand Binding Modes at α1-Adrenergic Receptors Explains the Molecular Basis of Subtype Selectivity

The structural poses of ligands that bind weakly to protein receptors are challenging to define. In this work we have studied ligand interactions with the adrenoreceptor (AR) subtypes, α_1A-AR and α_1B-AR, which belong to the G protein-coupled receptor (GPCR) superfamily, by employing the solution-based ligand-observed NMR method interligand NOEs for pharmacophore mapping (INPHARMA). A lack of receptor crystal structures and of subtype-selective drugs has hindered the definition of the physiological roles of each subtype and limited drug development. We determined the binding pose of the weakly binding α_1A-AR-selective agonist A-61603 relative to an endogenous agonist, epinephrine, at both α_1A-AR and α_1B-AR. The NMR experimental data were quantitatively compared, by using SpINPHARMA, to the back-calculated spectra based on ligand poses obtained from all-atom molecular dynamics simulations. The results helped mechanistically explain the selectivity of (R)-A-61603 towards α_1A-AR, thus demonstrating an approach for targeting subtype selectivity in ARs.     

The Structural Basis of Peptide Binding at Class A G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) represent the largest membrane protein family and a significant target class for therapeutics. Receptors from GPCRs’ largest class, class A, influence virtually every aspect of human physiology. About 45% of the members of this family endogenously bind flexible peptides or peptides segments within larger protein ligands. While many of these peptides have been structurally characterized in their solution state, the few studies of peptides in their receptor-bound state suggest that these peptides interact with a shared set of residues and undergo significant conformational changes. For the purpose of understanding binding dynamics and the development of peptidomimetic drug compounds, further studies should investigate the peptide ligands that are complexed to their cognate receptor.     

A single amino acid determines preference between phospholipids and reveals length restriction for activation ofthe S1P<Subscript>4</Subscript> receptor

Background

Sphingosine-1-phosphate and lysophosphatidic acid (LPA) are ligands for two related families of G protein-coupled receptors, the S1P and LPA receptors, respectively. The lysophospholipid ligands of these receptors are structurally similar, however recognition of these lipids by these receptors is highly selective. A single residue present within the third transmembrane domain (TM) of S1P receptors is thought to determine ligand selectivity; replacement of the naturally occurring glutamic acid with glutamine (present at this position in the LPA receptors) has previously been shown to be sufficient to change the specificity of S1P₁ from S1P to 18:1 LPA.

Results

We tested whether mutation of this "ligand selectivity" residue to glutamine could confer LPA-responsiveness to the related S1P receptor, S1P₄. This mutation severely affected the response of S1P₄ to S1P in a [³⁵S]GTPγS binding assay, and imparted sensitivity to LPA species in the order 14:0 LPA > 16:0 LPA > 18:1 LPA. These results indicate a length restriction for activation of this receptor and demonstrate the utility of using LPA-responsive S1P receptor mutants to probe binding pocket length using readily available LPA species. Computational modelling of the interactions between these ligands and both wild type and mutant S1P₄ receptors showed excellent agreement with experimental data, therefore confirming the fundamental role of this residue in ligand recognition by S1P receptors.

Conclusions

Glutamic acid in the third transmembrane domain of the S1P receptors is a general selectivity switch regulating response to S1P over the closely related phospholipids, LPA. Mutation of this residue to glutamine confers LPA responsiveness with preference for short-chain species. The preference for short-chain LPA species indicates a length restriction different from the closely related S1P₁ receptor.

     

Computational analysis of ligand recognition mechanisms by prostaglandin E2 (subtype 2) and D2 receptors

The eight members of the prostanoid receptor family belong to the class A G protein-coupled receptors. We investigated the evolutionary relationship of the eight members by a molecular phylogenetic analysis and found that prostaglandin E₂ receptor subtype 2 (EP₂) and prostaglandin D₂ receptor (DP) were closely related. The structures of the ligands for the two receptors are similar to each other but are distinguished by the exchanged locations of the carbonyl oxygen and the hydroxy group in the cyclopentane ring. The ligand recognition mechanisms of the receptors were examined by an integrated approach using several computational methods, such as amino acid sequence comparison, homology modeling, docking simulation, and molecular dynamics simulation. The results revealed the similar location of the ligand between the two receptors. The common carboxy group of the ligands interacts with the Arg residue on the seventh transmembrane (TM) helix, which is invariant among the prostanoid receptors. EP₂ uses a Ser on TM1 to recognize the carbonyl oxygen in the cyclopentane ring of the ligand. The Ser is specifically conserved within EP₂. On the other hand, DP uses a Lys on TM2 to recognize the hydroxy group of the ?? chain of the ligand. The Lys is also specifically conserved within DP. The interaction network between the D(E)RY motif and TM6 was found in EP₂. However, DP lacked this network, due to the mutation in the D(E)RY motif. Based on these observations and the previously published mutational studies on the motif, the possibility of another activation mechanism that does not involve the interaction between the D(E)RY motif and TM6 will be discussed.