The activities of drug inactive ingredients on biological targets.

Excipients, considered "inactive ingredients," are a major component of formulated drugs and play key roles in their pharmacokinetics. Despite their pervasiveness, whether they are active on any targets has not been systematically explored. We computed the likelihood that approved excipients would bind to molecular targets. Testing in vitro revealed 25 excipient activities, ranging from low-nanomolar to high-micromolar concentration. Another 109 activities were identified by testing against clinical safety targets. In cellular models, five excipients had fingerprints predictive of system-level toxicity. Exposures of seven excipients were investigated, and in certain populations, two of these may reach levels of in vitro target potency, including brain and gut exposure of thimerosal and its major metabolite, which had dopamine D3 receptor dissociation constant K d values of 320 and 210 nM, respectively. Although most excipients deserve their status as inert, many approved excipients may directly modulate physiologically relevant targets.

Virtual discovery of melatonin receptor ligands to modulate circadian rhythms.

The neuromodulator melatonin synchronizes circadian rhythms and related physiological functions through the actions of two G-protein-coupled receptors: MT1 and MT2. Circadian release of melatonin at night from the pineal gland activates melatonin receptors in the suprachiasmatic nucleus of the hypothalamus, synchronizing the physiology and behaviour of animals to the light-dark cycle1-4. The two receptors are established drug targets for aligning circadian phase to this cycle in disorders of sleep5,6 and depression1-4,7-9. Despite their importance, few in vivo active MT1-selective ligands have been reported2,8,10-12, hampering both the understanding of circadian biology and the development of targeted therapeutics. Here we docked more than 150 million virtual molecules to an MT1 crystal structure, prioritizing structural fit and chemical novelty. Of these compounds, 38 high-ranking molecules were synthesized and tested, revealing ligands with potencies ranging from 470 picomolar to 6 micromolar. Structure-based optimization led to two selective MT1 inverse agonists-which were topologically unrelated to previously explored chemotypes-that acted as inverse agonists in a mouse model of circadian re-entrainment. Notably, we found that these MT1-selective inverse agonists advanced the phase of the mouse circadian clock by 1.3-1.5 h when given at subjective dusk, an agonist-like effect that was eliminated in MT1- but not in MT2-knockout mice. This study illustrates the opportunities for modulating melatonin receptor biology through MT1-selective ligands and for the discovery of previously undescribed, in vivo active chemotypes from structure-based screens of diverse, ultralarge libraries.

Interactions of Oral Molecular Excipients with Breast Cancer Resistance Protein, BCRP.

Mechanistic-understanding-based selection of excipients may improve formulation development strategies for generic drug products and potentially accelerate their approval. Our study aimed at investigating the effects of molecular excipients present in orally administered FDA-approved drug products on the intestinal efflux transporter, BCRP (ABCG2), which plays a critical role in drug absorption with potential implications on drug safety and efficacy. We determined the interactions of 136 oral molecular excipients with BCRP in isolated membrane vesicles and identified 26 excipients as BCRP inhibitors with IC50 values less than 5 µM using 3H-cholecystokinin octapeptide (3H-CCK8). These BCRP inhibitors belonged to three functional categories of excipients: dyes, surfactants, and flavoring agents. Compared with noninhibitors, BCRP inhibitors had significantly higher molecular weights and SLogP values. The inhibitory effects of excipients identified in membrane vesicles were also evaluated in BCRP-overexpressing HEK293 cells at similar concentrations. Only 1 of the 26 inhibitors of BCRP identified in vesicles inhibited BCRP-mediated 3H-oxypurinol uptake by more than 50%, consistent with the notion that BCRP inhibition depends on transmembrane or intracellular availability of the inhibitors. Collectively, the results of this study provide new information on excipient selection during the development of drug products with active pharmaceutical ingredients that are BCRP substrates.

Selectivity Challenges in Docking Screens for GPCR Targets and Antitargets.

To investigate large library docking's ability to find molecules with joint activity against on-targets and selectivity versus antitargets, the dopamine D2 and serotonin 5-HT2A receptors were targeted, seeking selectivity against the histamine H1 receptor. In a second campaign, κ-opioid receptor ligands were sought with selectivity versus the µ-opioid receptor. While hit rates ranged from 40% to 63% against the on-targets, they were just as good against the antitargets, even though the molecules were selected for their putative lack of binding to the off-targets. Affinities, too, were often as good or better for the off-targets. Even though it was occasionally possible to find selective molecules, such as a mid-nanomolar D2/5-HT2A ligand with 21-fold selectivity versus the H1 receptor, this was the exception. Whereas false-negatives are tolerable in docking screens against on-targets, they are intolerable against antitargets; addressing this problem may demand new strategies in the field.

In silico design of novel probes for the atypical opioid receptor MRGPRX2.

The primate-exclusive MRGPRX2 G protein-coupled receptor (GPCR) has been suggested to modulate pain and itch. Despite putative peptide and small-molecule MRGPRX2 agonists, selective nanomolar-potency probes have not yet been reported. To identify a MRGPRX2 probe, we first screened 5,695 small molecules and found that many opioid compounds activated MRGPRX2, including (-)- and (+)-morphine, hydrocodone, sinomenine, dextromethorphan, and the prodynorphin-derived peptides dynorphin A, dynorphin B, and α- and ß-neoendorphin. We used these to select for mutagenesis-validated homology models and docked almost 4 million small molecules. From this docking, we predicted ZINC-3573-a potent MRGPRX2-selective agonist, showing little activity against 315 other GPCRs and 97 representative kinases-along with an essentially inactive enantiomer. ZINC-3573 activates endogenous MRGPRX2 in a human mast cell line, inducing degranulation and calcium release. MRGPRX2 is a unique atypical opioid-like receptor important for modulating mast cell degranulation, which can now be specifically modulated with ZINC-3573.

Virtual Screening for UDP-Galactopyranose Mutase Ligands Identifies a New Class of Antimycobacterial Agents.

Galactofuranose (Galf) is present in glycans critical for the virulence and viability of several pathogenic microbes, including Mycobacterium tuberculosis, yet the monosaccharide is absent from mammalian glycans. Uridine 5'-diphosphate-galactopyranose mutase (UGM) catalyzes the formation of UDP-Galf, which is required to produce Galf-containing glycoconjugates. Inhibitors of UGM have therefore been sought, both as antimicrobial leads and as tools to delineate the roles of Galf in cells. Obtaining cell permeable UGM probes by either design or high throughput screens has been difficult, as has elucidating how UGM binds small molecule, noncarbohydrate inhibitors. To address these issues, we employed structure-based virtual screening to uncover new inhibitor chemotypes, including a triazolothiadiazine series. These compounds are among the most potent antimycobacterial UGM inhibitors described. They also facilitated determination of a UGM-small molecule inhibitor structure, which can guide optimization. A comparison of results from the computational screen and a high-throughput fluorescence polarization (FP) screen indicated that the scaffold hits from the former had been evaluated in the FP screen but missed. By focusing on promising compounds, the virtual screen rescued false negatives, providing a blueprint for generating new UGM probes and therapeutic leads.

Colloidal aggregation and the in vitro activity of traditional Chinese medicines.

Traditional Chinese Medicines (TCMs) have been the sole source of therapeutics in China for two millennia. In recent drug discovery efforts, purified components of TCM formulations have shown activity in many in vitro assays, raising concerns of promiscuity. Here, we investigated 14 bioactive small molecules isolated from TCMs for colloidal aggregation. At concentrations commonly used in cell-based or biochemical assay conditions, eight of these compounds formed particles detectable by dynamic light scattering and showed detergent-reversible inhibition against ß-lactamase and malate dehydrogenase, two counter-screening enzymes. When three of these compounds were tested against their literature-reported molecular targets, they showed similar reversal of their inhibitory activity in the presence of detergent. For three of the most potent aggregators, contributions to promiscuity via oxidative cycling were investigated; addition of 1 mM DTT had no effect on their activity, which is inconsistent with an oxidative mechanism. TCMs are often active at micromolar concentrations; this study suggests that care must be taken to control for artifactual activity when seeking their primary targets. Implications for the formulation of these molecules are considered.

In silico molecular comparisons of C. elegans and mammalian pharmacology identify distinct targets that regulate feeding.

Phenotypic screens can identify molecules that are at once penetrant and active on the integrated circuitry of a whole cell or organism. These advantages are offset by the need to identify the targets underlying the phenotypes. Additionally, logistical considerations limit screening for certain physiological and behavioral phenotypes to organisms such as zebrafish and C. elegans. This further raises the challenge of elucidating whether compound-target relationships found in model organisms are preserved in humans. To address these challenges we searched for compounds that affect feeding behavior in C. elegans and sought to identify their molecular mechanisms of action. Here, we applied predictive chemoinformatics to small molecules previously identified in a C. elegans phenotypic screen likely to be enriched for feeding regulatory compounds. Based on the predictions, 16 of these compounds were tested in vitro against 20 mammalian targets. Of these, nine were active, with affinities ranging from 9 nM to 10 µM. Four of these nine compounds were found to alter feeding. We then verified the in vitro findings in vivo through genetic knockdowns, the use of previously characterized compounds with high affinity for the four targets, and chemical genetic epistasis, which is the effect of combined chemical and genetic perturbations on a phenotype relative to that of each perturbation in isolation. Our findings reveal four previously unrecognized pathways that regulate feeding in C. elegans with strong parallels in mammals. Together, our study addresses three inherent challenges in phenotypic screening: the identification of the molecular targets from a phenotypic screen, the confirmation of the in vivo relevance of these targets, and the evolutionary conservation and relevance of these targets to their human orthologs.

Conformation guides molecular efficacy in docking screens of activated ß-2 adrenergic G protein coupled receptor.

A prospective, large library virtual screen against an activated ß2-adrenergic receptor (ß2AR) structure returned potent agonists to the exclusion of inverse-agonists, providing the first complement to the previous virtual screening campaigns against inverse-agonist-bound G protein coupled receptor (GPCR) structures, which predicted only inverse-agonists. In addition, two hits recapitulated the signaling profile of the co-crystal ligand with respect to the G protein and arrestin mediated signaling. This functional fidelity has important implications in drug design, as the ability to predict ligands with predefined signaling properties is highly desirable. However, the agonist-bound state provides an uncertain template for modeling the activated conformation of other GPCRs, as a dopamine D2 receptor (DRD2) activated model templated on the activated ß2AR structure returned few hits of only marginal potency.

Large-scale prediction and testing of drug activity on side-effect targets.

Discovering the unintended 'off-targets' that predict adverse drug reactions is daunting by empirical methods alone. Drugs can act on several protein targets, some of which can be unrelated by conventional molecular metrics, and hundreds of proteins have been implicated in side effects. Here we use a computational strategy to predict the activity of 656 marketed drugs on 73 unintended 'side-effect' targets. Approximately half of the predictions were confirmed, either from proprietary databases unknown to the method or by new experimental assays. Affinities for these new off-targets ranged from 1 nM to 30 µM. To explore relevance, we developed an association metric to prioritize those new off-targets that explained side effects better than any known target of a given drug, creating a drug-target-adverse drug reaction network. Among these new associations was the prediction that the abdominal pain side effect of the synthetic oestrogen chlorotrianisene was mediated through its newly discovered inhibition of the enzyme cyclooxygenase-1. The clinical relevance of this inhibition was borne out in whole human blood platelet aggregation assays. This approach may have wide application to de-risking toxicological liabilities in drug discovery.

Chemical informatics and target identification in a zebrafish phenotypic screen.

Target identification is a core challenge in chemical genetics. Here we use chemical similarity to computationally predict the targets of 586 compounds that were active in a zebrafish behavioral assay. Among 20 predictions tested, 11 compounds had activities ranging from 1 nM to 10,000 nM on the predicted targets. The roles of two of these targets were tested in the original zebrafish phenotype. Prediction of targets from chemotype is rapid and may be generally applicable.

The presynaptic component of the serotonergic system is required for clozapine's efficacy.

Clozapine, by virtue of its absence of extrapyramidal side effects and greater efficacy, revolutionized the treatment of schizophrenia, although the mechanisms underlying this exceptional activity remain controversial. Combining an unbiased cheminformatics and physical screening approach, we evaluated clozapine's activity at >2350 distinct molecular targets. Clozapine, and the closely related atypical antipsychotic drug olanzapine, interacted potently with a unique spectrum of molecular targets. This distinct pattern, which was not shared with the typical antipsychotic drug haloperidol, suggested that the serotonergic neuronal system was a key determinant of clozapine's actions. To test this hypothesis, we used pet1(-/-) mice, which are deficient in serotonergic presynaptic markers. We discovered that the antipsychotic-like properties of the atypical antipsychotic drugs clozapine and olanzapine were abolished in a pharmacological model that mimics NMDA-receptor hypofunction in pet1(-/-) mice, whereas haloperidol's efficacy was unaffected. These results show that clozapine's ability to normalize NMDA-receptor hypofunction, which is characteristic of schizophrenia, depends on an intact presynaptic serotonergic neuronal system.

Quantitative analyses of aggregation, autofluorescence, and reactivity artifacts in a screen for inhibitors of a thiol protease.

The perceived and actual burden of false positives in high-throughput screening has received considerable attention; however, few studies exist on the contributions of distinct mechanisms of nonspecific effects like chemical reactivity, assay signal interference, and colloidal aggregation. Here, we analyze the outcome of a screen of 197861 diverse compounds in a concentration-response format against the cysteine protease cruzain, a target expected to be particularly sensitive to reactive compounds, and using an assay format with light detection in the short-wavelength region where significant compound autofluorescence is typically encountered. Approximately 1.9% of all compounds screened were detergent-sensitive inhibitors. The contribution from autofluorescence and compounds bearing reactive functionalities was dramatically lower: of all hits, only 1.8% were autofluorescent and 1.5% contained reactive or undesired functional groups. The distribution of false positives was relatively constant across library sources. The simple step of including detergent in the assay buffer suppressed the nonspecific effect of approximately 93% of the original hits.

Predicting new molecular targets for known drugs.

Although drugs are intended to be selective, at least some bind to several physiological targets, explaining side effects and efficacy. Because many drug-target combinations exist, it would be useful to explore possible interactions computationally. Here we compared 3,665 US Food and Drug Administration (FDA)-approved and investigational drugs against hundreds of targets, defining each target by its ligands. Chemical similarities between drugs and ligand sets predicted thousands of unanticipated associations. Thirty were tested experimentally, including the antagonism of the beta(1) receptor by the transporter inhibitor Prozac, the inhibition of the 5-hydroxytryptamine (5-HT) transporter by the ion channel drug Vadilex, and antagonism of the histamine H(4) receptor by the enzyme inhibitor Rescriptor. Overall, 23 new drug-target associations were confirmed, five of which were potent (<100 nM). The physiological relevance of one, the drug N,N-dimethyltryptamine (DMT) on serotonergic receptors, was confirmed in a knockout mouse. The chemical similarity approach is systematic and comprehensive, and may suggest side-effects and new indications for many drugs.

Automated docking screens: a feasibility study.

Molecular docking is the most practical approach to leverage protein structure for ligand discovery, but the technique retains important liabilities that make it challenging to deploy on a large scale. We have therefore created an expert system, DOCK Blaster, to investigate the feasibility of full automation. The method requires a PDB code, sometimes with a ligand structure, and from that alone can launch a full screen of large libraries. A critical feature is self-assessment, which estimates the anticipated reliability of the automated screening results using pose fidelity and enrichment. Against common benchmarks, DOCK Blaster recapitulates the crystal ligand pose within 2 A rmsd 50-60% of the time; inferior to an expert, but respectrable. Half the time the ligand also ranked among the top 5% of 100 physically matched decoys chosen on the fly. Further tests were undertaken culminating in a study of 7755 eligible PDB structures. In 1398 cases, the redocked ligand ranked in the top 5% of 100 property-matched decoys while also posing within 2 A rmsd, suggesting that unsupervised prospective docking is viable. DOCK Blaster is available at http://blaster.docking.org .

Comprehensive mechanistic analysis of hits from high-throughput and docking screens against beta-lactamase.

High-throughput screening (HTS) is widely used in drug discovery. Especially for screens of unbiased libraries, false positives can dominate "hit lists"; their origins are much debated. Here we determine the mechanism of every active hit from a screen of 70,563 unbiased molecules against beta-lactamase using quantitative HTS (qHTS). Of the 1,274 initial inhibitors, 95% were detergent-sensitive and were classified as aggregators. Among the 70 remaining were 25 potent, covalent-acting beta-lactams. Mass spectra, counter-screens, and crystallography identified 12 as promiscuous covalent inhibitors. The remaining 33 were either aggregators or irreproducible. No specific reversible inhibitors were found. We turned to molecular docking to prioritize molecules from the same library for testing at higher concentrations. Of 16 tested, 2 were modest inhibitors. Subsequent X-ray structures corresponded to the docking prediction. Analog synthesis improved affinity to 8 microM. These results suggest that it may be the physical behavior of organic molecules, not their reactivity, that accounts for most screening artifacts. Structure-based methods may prioritize weak-but-novel chemotypes in unbiased library screens.

High-throughput assays for promiscuous inhibitors.

High-throughput screening (HTS) searches large libraries of chemical compounds for those that can modulate the activity of a particular biological target; it is the dominant technique used in early-stage drug discovery. A key problem in HTS is the prevalence of nonspecific or 'promiscuous' inhibitors. These molecules have peculiar properties, act on unrelated targets and can dominate the results from screening campaigns. Several explanations have been proposed to account for promiscuous inhibitors, including chemical reactivity, interference in assay read-out, high molecular flexibility and hydrophobicity. The diversity of these models reflects the apparently unrelated molecules whose behaviors they seek to explain. However, a single mechanism may explain the effects of many promiscuous inhibitors: some organic molecules form large colloid-like aggregates that sequester and thereby inhibit enzymes. Hits from HTS, leads for drug discovery and even several drugs appear to act through this mechanism at micromolar concentrations. Here, we report two rapid assays for detecting promiscuous aggregates that we tested against 1,030 'drug-like' molecules. The results from these assays were used to test two preliminary computational models of this phenomenon and as benchmarks to develop new models.

Virtual screening of chemical libraries.

Virtual screening uses computer-based methods to discover new ligands on the basis of biological structures. Although widely heralded in the 1970s and 1980s, the technique has since struggled to meet its initial promise, and drug discovery remains dominated by empirical screening. Recent successes in predicting new ligands and their receptor-bound structures, and better rates of ligand discovery compared to empirical screening, have re-ignited interest in virtual screening, which is now widely used in drug discovery, albeit on a more limited scale than empirical screening.

Molecular docking and high-throughput screening for novel inhibitors of protein tyrosine phosphatase-1B.

High-throughput screening (HTS) of compound libraries is used to discover novel leads for drug development. When a structure is available for the target, computer-based screening using molecular docking may also be considered. The two techniques have rarely been used together on the same target. The opportunity to do so presented itself in a project to discover novel inhibitors for the enzyme protein tyrosine phosphatase-1B (PTP1B), a tyrosine phosphatase that has been implicated as a key target for type II diabetes. A corporate library of approximately 400 000 compounds was screened using high-throughput experimental techniques for compounds that inhibited PTP1B. Concurrently, molecular docking was used to screen approximately 235 000 commercially available compounds against the X-ray crystallographic structure of PTP1B, and 365 high-scoring molecules were tested as inhibitors of the enzyme. Of approximately 400 000 molecules tested in the high-throughput experimental assay, 85 (0.021%) inhibited the enzyme with IC50 values less than 100 microM; the most active had an IC50 value of 4.2 microM. Of the 365 molecules suggested by molecular docking, 127 (34.8%) inhibited PTP1B with IC50 values less than 100 microM; the most active of these had an IC50 of 1.7 microM. Structure-based docking therefore enriched the hit rate by 1700-fold over random screening. The hits from both the high-throughput and docking screens were dissimilar from phosphotyrosine, the canonical substrate group for PTP1B; the two hit lists were also very different from each other. Surprisingly, the docking hits were judged to be more druglike than the HTS hits. The diversity of both hit lists and their dissimilarity from each other suggest that docking and HTS may be complementary techniques for lead discovery.