Global Analysis of Small Molecule Binding to Related Protein Targets

We report on the integration of pharmacological data and homology information for a large scale analysis of small molecule binding to related targets. Differences in small molecule binding have been assessed for curated pairs of human to rat orthologs and also for recently diverged human paralogs. Our analysis shows that in general, small molecule binding is conserved for pairs of human to rat orthologs. Using statistical tests, we identified a small number of cases where small molecule binding is different between human and rat, some of which had previously been reported in the literature. Knowledge of species specific pharmacology can be advantageous for drug discovery, where rats are frequently used as a model system. For human paralogs, we demonstrate a global correlation between sequence identity and the binding of small molecules with equivalent affinity. Our findings provide an initial general model relating small molecule binding and sequence divergence, containing the foundations for a general model to anticipate and predict within-target-family selectivity

The role of kinetic context in apparent biased agonism at GPCRs

Biased agonism describes the ability of ligands to stabilize different conformations of a GPCR linked to distinct functional outcomes and offers the prospect of designing pathway-specific drugs that avoid on-target side effects. This mechanism is usually inferred from pharmacological data with the assumption that the confounding influences of observational (that is, assay dependent) and system (that is, cell background dependent) bias are excluded by experimental design and analysis. Here we reveal that ‘kinetic context’, as determined by ligand-binding kinetics and the temporal pattern of receptor-signalling processes, can have a profound influence on the apparent bias of a series of agonists for the dopamine D2 receptor and can even lead to reversals in the direction of bias. We propose that kinetic context must be acknowledged in the design and interpretation of studies of biased agonism

Using Multiple Microenvironments to Find Similar Ligand-Binding Sites: Application to Kinase Inhibitor Binding

The recognition of cryptic small-molecular binding sites in protein structures is important for understanding off-target side effects and for recognizing potential new indications for existing drugs. Current methods focus on the geometry and detailed chemical interactions within putative binding pockets, but may not recognize distant similarities where dynamics or modified interactions allow one ligand to bind apparently divergent binding pockets. In this paper, we introduce an algorithm that seeks similar microenvironments within two binding sites, and assesses overall binding site similarity by the presence of multiple shared microenvironments. The method has relatively weak geometric requirements (to allow for conformational change or dynamics in both the ligand and the pocket) and uses multiple biophysical and biochemical measures to characterize the microenvironments (to allow for diverse modes of ligand binding). We term the algorithm PocketFEATURE, since it focuses on pockets using the FEATURE system for characterizing microenvironments. We validate PocketFEATURE first by showing that it can better discriminate sites that bind similar ligands from those that do not, and by showing that we can recognize FAD-binding sites on a proteome scale with Area Under the Curve (AUC) of 92%. We then apply PocketFEATURE to evolutionarily distant kinases, for which the method recognizes several proven distant relationships, and predicts unexpected shared ligand binding. Using experimental data from ChEMBL and Ambit, we show that at high significance level, 40 kinase pairs are predicted to share ligands. Some of these pairs offer new opportunities for inhibiting two proteins in a single pathway

Whole-exome sequencing uncovers frequent GNAS mutations in intraductal papillary mucinous neoplasms of the pancreas

Intraductal papillary mucinous neoplasm (IPMN) is a common pancreatic cystic neoplasm that is often invasive and metastatic, resulting in a poor prognosis. Few molecular alterations unique to IPMN are known. We performed whole-exome sequencing for a primary IPMN tissue, which uncovered somatic mutations in KCNF1, DYNC1H1, PGCP, STAB1, PTPRM, PRPF8, RNASE3, SPHKAP, MLXIPL, VPS13C, PRCC, GNAS, KRAS, RBM10, RNF43, DOCK2, and CENPF. We further analyzed GNAS mutations in archival cases of 118 IPMNs and 32 pancreatic ductal adenocarcinomas (PDAs), which revealed that 48 (40.7%) of the 118 IPMNs but none of the 32 PDAs harbored GNAS mutations. G-protein alpha-subunit encoded by GNAS and its downstream targets, phosphorylated substrates of protein kinase A, were evidently expressed in IPMN; the latter was associated with neoplastic grade. These results indicate that GNAS mutations are common and specific for IPMN, and activation of G-protein signaling appears to play a pivotal role in IPMN

TargetMine, an Integrated Data Warehouse for Candidate Gene Prioritisation and Target Discovery

Prioritising candidate genes for further experimental characterisation is a non-trivial challenge in drug discovery and biomedical research in general. An integrated approach that combines results from multiple data types is best suited for optimal target selection. We developed TargetMine, a data warehouse for efficient target prioritisation. TargetMine utilises the InterMine framework, with new data models such as protein-DNA interactions integrated in a novel way. It enables complicated searches that are difficult to perform with existing tools and it also offers integration of custom annotations and in-house experimental data. We proposed an objective protocol for target prioritisation using TargetMine and set up a benchmarking procedure to evaluate its performance. The results show that the protocol can identify known disease-associated genes with high precision and coverage. A demonstration version of TargetMine is available at http://targetmine.nibio.go.jp/

Structural genomics target selection for the New York consortium on membrane protein structure

The New York Consortium on Membrane Protein Structure (NYCOMPS), a part of the Protein Structure Initiative (PSI) in the USA, has as its mission to establish a high-throughput pipeline for determination of novel integral membrane protein structures. Here we describe our current target selection protocol, which applies structural genomics approaches informed by the collective experience of our team of investigators. We first extract all annotated proteins from our reagent genomes, i.e. the 96 fully sequenced prokaryotic genomes from which we clone DNA. We filter this initial pool of sequences and obtain a list of valid targets. NYCOMPS defines valid targets as those that, among other features, have at least two predicted transmembrane helices, no predicted long disordered regions and, except for community nominated targets, no significant sequence similarity in the predicted transmembrane region to any known protein structure. Proteins that feed our experimental pipeline are selected by defining a protein seed and searching the set of all valid targets for proteins that are likely to have a transmembrane region structurally similar to that of the seed. We require sequence similarity aligning at least half of the predicted transmembrane region of seed and target. Seeds are selected according to their feasibility and/or biological interest, and they include both centrally selected targets and community nominated targets. As of December 2008, over 6,000 targets have been selected and are currently being processed by the experimental pipeline. We discuss how our target list may impact structural coverage of the membrane protein space

Simulation-based cheminformatic analysis of organelle-targeted molecules: lysosomotropic monobasic amines

Cell-based molecular transport simulations are being developed to facilitate exploratory cheminformatic analysis of virtual libraries of small drug-like molecules. For this purpose, mathematical models of single cells are built from equations capturing the transport of small molecules across membranes. In turn, physicochemical properties of small molecules can be used as input to simulate intracellular drug distribution, through time. Here, with mathematical equations and biological parameters adjusted so as to mimic a leukocyte in the blood, simulations were performed to analyze steady state, relative accumulation of small molecules in lysosomes, mitochondria, and cytosol of this target cell, in the presence of a homogenous extracellular drug concentration. Similarly, with equations and parameters set to mimic an intestinal epithelial cell, simulations were also performed to analyze steady state, relative distribution and transcellular permeability in this non-target cell, in the presence of an apical-to-basolateral concentration gradient. With a test set of ninety-nine monobasic amines gathered from the scientific literature, simulation results helped analyze relationships between the chemical diversity of these molecules and their intracellular distributions

A Structure-Based Approach for Mapping Adverse Drug Reactions to the Perturbation of Underlying Biological Pathways

Adverse drug reactions (ADR), also known as side-effects, are complex undesired physiologic phenomena observed secondary to the administration of pharmaceuticals. Several phenomena underlie the emergence of each ADR; however, a dominant factor is the drug's ability to modulate one or more biological pathways. Understanding the biological processes behind the occurrence of ADRs would lead to the development of safer and more effective drugs. At present, no method exists to discover these ADR-pathway associations. In this paper we introduce a computational framework for identifying a subset of these associations based on the assumption that drugs capable of modulating the same pathway may induce similar ADRs. Our model exploits multiple information resources. First, we utilize a publicly available dataset pairing drugs with their observed ADRs. Second, we identify putative protein targets for each drug using the protein structure database and in-silico virtual docking. Third, we label each protein target with its known involvement in one or more biological pathways. Finally, the relationships among these information sources are mined using multiple stages of logistic-regression while controlling for over-fitting and multiple-hypothesis testing. As proof-of-concept, we examined a dataset of 506 ADRs, 730 drugs, and 830 human protein targets. Our method yielded 185 ADR-pathway associations of which 45 were selected to undergo a manual literature review. We found 32 associations to be supported by the scientific literature

A two-step target binding and selectivity support vector machines approach for virtual screening of dopamine receptor subtype-selective ligands

10.1371/journal.pone.0039076PLoS ONE76

G protein-coupled receptor-mediated calcium signaling in astrocytes

Astrocytes express a large variety of G~protein-coupled receptors (GPCRs) which mediate the transduction of extracellular signals into intracellular calcium responses. This transduction is provided by a complex network of biochemical reactions which mobilizes a wealth of possible calcium-mobilizing second messenger molecules. Inositol 1,4,5-trisphosphate is probably the best known of these molecules whose enzymes for its production and degradation are nonetheless calcium-dependent. We present a biophysical modeling approach based on the assumption of Michaelis-Menten enzyme kinetics, to effectively describe GPCR-mediated astrocytic calcium signals. Our model is then used to study different mechanisms at play in stimulus encoding by shape and frequency of calcium oscillations in astrocytes.Comment: 35 pages, 6 figures, 1 table, 3 appendices (book chapter