An Ultra High-Throughput,Whole-Animal Screen for Small Molecule Modulators of a Specific Genetic Pathway in Caenorhabditis elegans

High-throughput screening (HTS) is a powerful approach to drug discovery, but many lead compounds are found to be unsuitable for use in vivo after initial screening. Screening in small animals like C. elegans can help avoid these problems, but this system has been limited to screens with low-throughput or no specific molecular target. We report the first in vivo 1536-well plate assay for a specific genetic pathway in C. elegans. Our assay measures induction of a gene regulated by SKN-1, a master regulator of detoxification genes. SKN-1 inhibitors will be used to study and potentially reverse multidrug resistance in parasitic nematodes. Screens of two small commercial libraries and the full Molecular Libraries Small Molecule Repository (MLSMR) of ∼364,000 compounds validate our platform for ultra HTS. Our platform overcomes current limitations of many whole-animal screens and can be widely adopted for other inducible genetic pathways in nematodes and humans.     

Deletion of Dystrophin In-Frame Exon 5 Leads to a Severe Phenotype: Guidance for Exon Skipping Strategies

Duchenne and Becker muscular dystrophy severity depends upon the nature and location of the DMD gene lesion and generally correlates with the dystrophin open reading frame. However, there are striking exceptions where an in-frame genomic deletion leads to severe pathology or protein-truncating mutations (nonsense or frame-shifting indels) manifest as mild disease. Exceptions to the dystrophin reading frame rule are usually resolved after molecular diagnosis on muscle RNA. We report a moderate/severe Becker muscular dystrophy patient with an in-frame genomic deletion of DMD exon 5. This mutation has been reported by others as resulting in Duchenne or Intermediate muscular dystrophy, and the loss of this in-frame exon in one patient led to multiple splicing events, including omission of exon 6, that disrupts the open reading frame and is consistent with a severe phenotype. The patient described has a deletion of dystrophin exon 5 that does not compromise recognition of exon 6, and although the deletion does not disrupt the reading frame, his clinical presentation is more severe than would be expected for classical Becker muscular dystrophy. We suggest that the dystrophin isoform lacking the actin-binding sequence encoded by exon 5 is compromised, reflected by the phenotype resulting from induction of this dystrophin isoform in mouse muscle in vivo. Hence, exon skipping to address DMD-causing mutations within DMD exon 5 may not yield an isoform that confers marked clinical benefit. Additional studies will be required to determine whether multi-exon skipping strategies could yield more functional dystrophin isoforms, since some BMD patients with larger in-frame deletions in this region have been reported with mild phenotypes.     

Identification of Small Molecules with Type I Interferon Inducing Properties by High-Throughput Screening

The continuous emergence of virus that are resistant to current anti-viral drugs, combined with the introduction of new viral pathogens for which no therapeutics are available, creates an urgent need for the development of novel broad spectrum antivirals. Type I interferon (IFN) can, by modulating the cellular expression profile, stimulate a non-specific antiviral state. The antiviral and adjuvant properties of IFN have been extensively demonstrated; however, its clinical application has been so far limited. We have developed a human cell-based assay that monitors IFN-β production for use in a high throughput screen. Using this assay we screened 94,398 small molecules and identified 18 compounds with IFN-inducing properties. Among these, 3 small molecules (C3, E51 and L56) showed activity not only in human but also in murine and canine derived cells. We further characterized C3 and showed that this molecule is capable of stimulating an anti-viral state in human-derived lung epithelial cells. Furthermore, the IFN-induction by C3 is not diminished by the presence of influenza A virus NS1 protein or hepatitis C virus NS3/4A protease, which make this molecule an interesting candidate for the development of a new type of broad-spectrum antiviral. In addition, the IFN-inducing properties of C3 also suggest its potential use as vaccine adjuvant.     

High-Throughput Screening of Myometrial Calcium-Mobilization to Identify Modulators of Uterine Contractility

The uterine myometrium (UT-myo) is a therapeutic target for preterm labor, labor induction, and postpartum hemorrhage. Stimulation of intracellular Ca²⁺-release in UT-myo cells by oxytocin is a final pathway controlling myometrial contractions. The goal of this study was to develop a dual-addition assay for high-throughput screening of small molecular compounds, which could regulate Ca²⁺-mobilization in UT-myo cells, and hence, myometrial contractions. Primary murine UT-myo cells in 384-well plates were loaded with a Ca²⁺-sensitive fluorescent probe, and then screened for inducers of Ca²⁺-mobilization and inhibitors of oxytocin-induced Ca²⁺-mobilization. The assay exhibited robust screening statistics (Z´ = 0.73), DMSO-tolerance, and was validated for high-throughput screening against 2,727 small molecules from the Spectrum, NIH Clinical I and II collections of well-annotated compounds. The screen revealed a hit-rate of 1.80% for agonist and 1.39% for antagonist compounds. Concentration-dependent responses of hit-compounds demonstrated an EC₅₀ less than 10μM for 21 hit-antagonist compounds, compared to only 7 hit-agonist compounds. Subsequent studies focused on hit-antagonist compounds. Based on the percent inhibition and functional annotation analyses, we selected 4 confirmed hit-antagonist compounds (benzbromarone, dipyridamole, fenoterol hydrobromide and nisoldipine) for further analysis. Using an ex vivo isometric contractility assay, each compound significantly inhibited uterine contractility, at different potencies (IC₅₀). Overall, these results demonstrate for the first time that high-throughput small-molecules screening of myometrial Ca²⁺-mobilization is an ideal primary approach for discovering modulators of uterine contractility.     

A High-Throughput Assay for Small Molecule Destabilizers of the KRAS Oncoprotein

Mutations in the Ras family of small GTPases, particularly KRAS, occur at high frequencies in cancer and represent a major unmet therapeutic need due to the lack of effective targeted therapies. Past efforts directed at inhibiting the activity of the Ras oncoprotein have proved difficult. We propose an alternative approach to target Ras by eliminating Ras protein from cells with pharmacological means. In this study, we developed a cell-based, high-content screening platform to identify small molecules that could promote the degradation of the KRAS oncoprotein. We generated an EGFP-KRAS^G12V fluorescence reporter system and implemented it for automated screening in 1536-well plates using high-throughput cellular imaging. We screened a library of clinically relevant compounds at wide dose range and identified Ponatinib and AMG-47a as two candidate compounds that selectively reduced the levels of EGFP-KRAS^G12V protein but did not affect EGFP protein in cells. This proof-of-principle study demonstrates that it is feasible to use a high-throughput screen to identify compounds that promote the degradation of the Ras oncoprotein as a new approach to target Ras.     

Effects of Small Molecule Modulators on ATP Binding to Skeletal Ryanodine Receptor

Calcium release for muscle contraction in skeletal muscle is mediated in part by the ryanodine receptor 1, RyR1, Ca²⁺-channel and is strongly affected by intrinsic modulators like Ca²⁺, Mg²⁺ and ATP. We showed differential effects on ATP binding in the presence of Ca²⁺ or Mg²⁺ ions using ESR spectroscopy and a spin-labeled ATP analog, SL-ATP (Dias et al. Biochemistry 45: 9408–9415, 2006). We here report the effects of RyR1 modulators like ryanodine, caffeine and dantrolene on the ATP binding of RyR1 using the same technique. We present evidence that the exogenous effectors induce changes within RyR1 that lead to different ATP binding characteristics: In the presence of the activating modulator, caffeine, or in the presence of ryanodine, which causes a half-open state of the channel, binding of eight ATP per RyR1 was observed, even in the presence of inhibitory Ca²⁺, suggestive of a stable “open” channel conformation. In the presence of the inhibitory modulator dantrolene, ATP binding affinity decreased in the presence of activating Ca²⁺, while in the presence of inhibitory Ca²⁺, ATP binding affinity increased, but at the same time the number of accessible sites decreased to four, suggestive of a closed conformation of the channel. The results imply that modulation of ATP binding to RyR1 as well as the overall number of accessible ATP binding sites on the channel are crucial for regulation and are in direct correlation with the modified activity of the channel induced by pharmacological agents.     

Identification of Potent Chemotypes Targeting Leishmania major Using a High-Throughput,Low-Stringency,Computationally Enhanced,Small Molecule Screen

Patients with clinical manifestations of leishmaniasis, including cutaneous leishmaniasis, have limited treatment options, and existing therapies frequently have significant untoward liabilities. Rapid expansion in the diversity of available cutaneous leishmanicidal chemotypes is the initial step in finding alternative efficacious treatments. To this end, we combined a low-stringency Leishmania major promastigote growth inhibition assay with a structural computational filtering algorithm. After a rigorous assay validation process, we interrogated ∼200,000 unique compounds for L. major promastigote growth inhibition. Using iterative computational filtering of the compounds exhibiting >50% inhibition, we identified 553 structural clusters and 640 compound singletons. Secondary confirmation assays yielded 93 compounds with EC₅₀s ≤ 1 µM, with none of the identified chemotypes being structurally similar to known leishmanicidals and most having favorable in silico predicted bioavailability characteristics. The leishmanicidal activity of a representative subset of 15 chemotypes was confirmed in two independent assay formats, and L. major parasite specificity was demonstrated by assaying against a panel of human cell lines. Thirteen chemotypes inhibited the growth of a L. major axenic amastigote-like population. Murine in vivo efficacy studies using one of the new chemotypes document inhibition of footpad lesion development. These results authenticate that low stringency, large-scale compound screening combined with computational structure filtering can rapidly expand the chemotypes targeting in vitro and in vivo Leishmania growth and viability.     

Identification of Pharmacological Modulators of HMGB1-Induced Inflammatory Response by Cell-Based Screening

High mobility group box 1 (HMGB1), a highly conserved, ubiquitous protein, is released into the circulation during sterile inflammation (e.g. arthritis, trauma) and circulatory shock. It participates in the pathogenesis of delayed inflammatory responses and organ dysfunction. While several molecules have been identified that modulate the release of HMGB1, less attention has been paid to identify pharmacological inhibitors of the downstream inflammatory processes elicited by HMGB1 (C23-C45 disulfide C106 thiol form). In the current study, a cell-based medium-throughput screening of a 5000+ compound focused library of clinical drugs and drug-like compounds was performed in murine RAW264.7 macrophages, in order to identify modulators of HMGB1-induced tumor-necrosis factor alpha (TNFα) production. Clinically used drugs that suppressed HMGB1-induced TNFα production included glucocorticoids, beta agonists, and the anti-HIV compound indinavir. A re-screen of the NIH clinical compound library identified beta-agonists and various intracellular cAMP enhancers as compounds that potentiate the inhibitory effect of glucocorticoids on HMGB1-induced TNFα production. The molecular pathways involved in this synergistic anti-inflammatory effect are related, at least in part, to inhibition of TNFα mRNA synthesis via a synergistic suppression of ERK/IκB activation. Inhibition of TNFα production by prednisolone+salbutamol pretreatment was also confirmed in vivo in mice subjected to HMGB1 injection; this effect was more pronounced than the effect of either of the agents administered separately. The current study unveils several drug-like modulators of HMGB1-mediated inflammatory responses and offers pharmacological directions for the therapeutic suppression of inflammatory responses in HMGB1-dependent diseases.     

Identification of Small Molecule Activators of BMP Signaling

Bone Morphogenetic Proteins (BMPs) are morphogens that play a major role in regulating development and homeostasis. Although BMPs are used for the treatment of bone and kidney disorders, their clinical use is limited due to the supra-physiological doses required for therapeutic efficacy causing severe side effects. Because recombinant BMPs are expensive to produce, small molecule activators of BMP signaling would be a cost-effective alternative with the added benefit of being potentially more easily deliverable. Here, we report our efforts to identify small molecule activators of BMP signaling. We have developed a cell-based assay to monitor BMP signaling by stably transfecting a BMP-responsive human cervical carcinoma cell line (C33A) with a reporter construct in which the expression of luciferase is driven by a multimerized BMP-responsive element from the Id1 promoter. A BMP-responsive clone C33A-2D2 was used to screen a bioactive library containing ∼5,600 small molecules. We identified four small molecules of the family of flavonoids all of which induced luciferase activity in a dose-dependent manner and ventralized zebrafish embryos. Two of the identified compounds induced Smad1, 5 phosphorylation (P-Smad), Id1 and Id2 expression in a dose-dependent manner demonstrating that our assays identified small molecule activators of BMP signaling.     

A Genome-Wide CRISPR Library for High-Throughput Genetic Screening in Drosophila Cells

The simplicity of the CRISPR/Cas9 system of genome engineering has opened up the possibility of performing genome-wide targeted mutagenesis in cell lines,enabling screening for cellular phenotypes resulting from genetic aberrations.Drosophila cells have proven to be highly effective in identifying genes involved in cellular processes through similar screens using partial knockdown by RNAi.This is in part due to the lower degree of redundancy between genes in this organism,whilst still maintaining highly conserved gene networks and orthologs of many human disease-causing genes.The ability of CRISPR to generate genetic loss of function mutations not only increases the magnitude of any effect over currently employed RNAi techniques,but allows analysis over longer periods of time which can be critical for certain phenotypes.In this study,we have designed and built a genome-wide CRISPR library covering 13,501 genes,among which 8989 genes are targeted by three or more independent single guide RNAs(sg RNAs).Moreover,we describe strategies to monitor the population of guide RNAs by high throughput sequencing(HTS).We hope that this library will provide an invaluable resource for the community to screen loss of function mutations for cellular phenotypes,and as a source of guide RNA designs for future studies.     

Identification by Automated Screening of a Small Molecule that Selectively Eliminates Neural Stem Cells Derived from hESCs but Not Dopamine Neurons

Background

We have previously described fundamental differences in the biology of stem cells as compared to other dividing cell populations. We reasoned therefore that a differential screen using US Food and Drug Administration (FDA)-approved compounds may identify either selective survival factors or specific toxins and may be useful for the therapeutically-driven manufacturing of cells in vitro and possibly in vivo.

Methodology/Principal Findings

In this study we report on optimized methods for feeder-free culture of hESCs and hESC-derived neural stem cells (NSCs) to facilitate automated screening. We show that we are able to measure ATP as an indicator of metabolic activity in an automated screening assay. With this optimized platform we screened a collection of FDA-approved drugs to identify compounds that have differential toxicity to hESCs and their neural derivatives. Nine compounds were identified to be specifically toxic for NSCs to a greater extent than for hESCs. Six of these initial hits were retested and verified by large-scale cell culture to determine dose-responsive NSC toxicity. One of the compounds retested, amiodarone HCL, was further tested for possible effects on postmitotic neurons, a likely target for transplant therapy. Amiodarone HCL was found to be selectively toxic to NSCs but not to differentiated neurons or glial cells. Treated and untreated NSCs and neurons were then interrogated with global gene expression analysis to explore the mechanisms of action of amiodarone HCl. The gene expression analysis suggests that activation of cell-type specific cationic channels may underlie the toxicity of the drug.

Conclusions/Significance

In conclusion, we have developed a screening strategy that allows us to rapidly identify clinically approved drugs for use in a Chemistry, Manufacture and Control protocol that can be safely used to deplete unwanted contaminating precursor cells from a differentiated cell product. Our results also suggest that such a strategy is rich in the potential of identifying lineage specific reagents and provides additional evidence for the utility of stem cells in screening and discovery paradigms.     

An In Silico Analysis of the Binding Modes and Binding Affinities of Small Molecule Modulators of PDZ-Peptide Interactions

Inhibitors of PDZ-peptide interactions have important implications in a variety of biological processes including treatment of cancer and Parkinson’s disease. Even though experimental studies have reported characterization of peptidomimetic inhibitors of PDZ-peptide interactions, the binding modes for most of them have not been characterized by structural studies. In this study we have attempted to understand the structural basis of the small molecule-PDZ interactions by in silico analysis of the binding modes and binding affinities of a set of 38 small molecules with known K_i or K_d values for PDZ2 and PDZ3 domains of PSD-95 protein. These two PDZ domains show differential selectivity for these compounds despite having a high degree of sequence similarity and almost identical peptide binding pockets. Optimum binding modes for these ligands for PDZ2 and PDZ3 domains were identified by using a novel combination of semi-flexible docking and explicit solvent molecular dynamics (MD) simulations. Analysis of the binding modes revealed most of the peptidomimectic ligands which had high K_i or K_d moved away from the peptide binding pocket, while ligands with high binding affinities remained in the peptide binding pocket. The differential specificities of the PDZ2 and PDZ3 domains primarily arise from differences in the conformation of the loop connecting βB and βC strands, because this loop interacts with the N-terminal chemical moieties of the ligands. We have also computed the MM/PBSA binding free energy values for these 38 compounds with both the PDZ domains from multiple 5 ns MD trajectories on each complex i.e. a total of 228 MD trajectories of 5 ns length each. Interestingly, computational binding free energies show good agreement with experimental binding free energies with a correlation coefficient of approximately 0.6. Thus our study demonstrates that combined use of docking and MD simulations can help in identification of potent inhibitors of PDZ-peptide complexes.