Aqueous solubility is one of the most important properties in drug discovery, as it has profound impact on various drug properties, including biological activity, pharmacokinetics (PK), toxicity, and in vivo efficacy. Both kinetic and thermodynamic solubilities are determined during different stages of drug discovery and development. Since kinetic solubility is more relevant in preclinical drug discovery research, especially during the structure optimization process, we have developed predictive models for kinetic solubility with in-house data generated from 11,780 compounds collected from over 200 NCATS intramural research projects. This represents one of the largest kinetic solubility datasets of high quality and integrity. Based on the customized atom type descriptors, the support vector classification (SVC) models were trained on 80% of the whole dataset, and exhibited high predictive performance for estimating the solubility of the remaining 20% compounds within the test set. The values of the area under the receiver operating characteristic curve (AUC-ROC) for the compounds in the test sets reached 0.93 and 0.91, when the threshold for insoluble compounds was set to 10 and 50 μg/mL respectively. The predictive models of aqueous solubility can be used to identify insoluble compounds in drug discovery pipeline, provide design ideas for improving solubility by analyzing the atom types associated with poor solubility and prioritize compound libraries to be purchased or synthesized.

Chemical genomics aims to comprehensively define, and ultimately predict, the effects of small molecule compounds on biological systems. Chemical activity profiling approaches must consider chemical effects on all pathways operative in mammalian cells. To enable a strategic and maximally efficient chemical profiling of pathway space, we have created the NCATS BioPlanet, a comprehensive integrated pathway resource that incorporates the universe of 1,658 human pathways sourced from publicly available, manually curated sources, which have been subjected to thorough redundancy and consistency cross-evaluation. BioPlanet supports interactive browsing, retrieval, and analysis of pathways, exploration of pathway connections, and pathway search by gene targets, category, and availability of corresponding bioactivity assay, as well as visualization of pathways on a 3-dimensional globe, in which the distance between any two pathways is proportional to their degree of gene component overlap. Using this resource, we propose a strategy to identify a minimal set of 362 biological assays that can interrogate the universe of human pathways. The NCATS BioPlanet is a public resource, which will be continually expanded and updated, for systems biology, toxicology, and chemical genomics, available at http://tripod.nih.gov/bioplanet/.

Aldehyde dehydrogenases (ALDHs) are responsible for the metabolism of aldehydes (exogenous and endogenous) and possess vital physiological and toxicological functions in areas such as CNS, inflammation, metabolic disorders, and cancers. Overexpression of certain ALDHs (e.g., ALDH1A1) is an important biomarker in cancers and cancer stem cells (CSCs) indicating the potential need for the identification and development of small molecule ALDH inhibitors. Herein, a newly designed series of quinoline-based analogs of ALDH1A1 inhibitors is described. Extensive medicinal chemistry optimization and biological characterization led to the identification of analogs with significantly improved enzymatic and cellular ALDH inhibition. Selected analogs, e.g., 86 (NCT-505) and 91 (NCT-506), demonstrated target engagement in a cellular thermal shift assay (CETSA), inhibited the formation of 3D spheroid cultures of OV-90 cancer cells, and potentiated the cytotoxicity of paclitaxel in SKOV-3-TR, a paclitaxel resistant ovarian cancer cell line. Lead compounds also exhibit high specificity over other ALDH isozymes and unrelated dehydrogenases. The in vitro ADME profiles and pharmacokinetic evaluation of selected analogs are also highlighted.

A chemical genomics "Toxmatrix" method was developed to elucidate mechanisms of cytotoxicity using neuronal models. Quantitative high-throughput screening (qHTS) was applied to systematically screen each toxicant against a panel of 70 modulators, drugs or chemicals that act on a known target, to identify interactions that either protect or sensitize cells to each toxicant. Thirty-two toxicants were tested at 10 concentrations for cytotoxicity to SH-SY5Y human neuroblastoma cells, with results fitted to the Hill equation to determine an IC50 for each toxicant. Thirty-three toxicant:modulator interactions were identified in SH-SY5Y cells for 14 toxicants, as modulators that shifted toxicant IC50 values lower or higher. The target of each modulator that sensitizes cells or protects cells from a toxicant suggests a mode of toxicant action or cellular adaptation. In secondary screening, we tested modulator-toxicant pairs identified from the SH-SY5Y primary screening for interactions in three differentiated neuronal human cell lines: dSH-SY5Y, conditionally immortalized dopaminergic neurons (LUHMES), and neural stem cells. Twenty toxicant-modulator pairs showed pronounced interactions in one or several differentiated cell models. Additional testing confirmed that several modulators acted through their primary targets. For example, several chelators protected differentiated LUHMES neurons from four toxicants by chelation of divalent cations and buthionine sulphoximine sensitized cells to 6-hydroxydopamine and 4-(methylamino)phenol hemisulfate by blocking glutathione synthesis. Such modulators that interact with multiple neurotoxicants suggest these may be vulnerable toxicity pathways in neurons. Thus, the Toxmatrix method is a systematic high-throughput approach that can identify mechanisms of toxicity and cellular adaptation.

Aldehyde dehydrogenase enzymes (ALDHs) have a broad spectrum of biological activities through the oxidation of both endogenous and exogenous aldehydes. Increased expression of ALDH1A1 has been identified in a wide-range of human cancer stem cells and is associated with cancer relapse and poor prognosis, raising the potential of ALDH1A1 as a therapeutic target. To facilitate quantitative high-throughput screening (qHTS) campaigns for the discovery, characterization and structure-activity-relationship (SAR) studies of small molecule ALDH1A1 inhibitors with cellular activity, we show herein the miniaturization to 1536-well format and automation of a high-content cell-based ALDEFLUOR assay. We demonstrate the utility of this assay by generating dose-response curves on a comprehensive set of prior art inhibitors as well as hundreds of ALDH1A1 inhibitors synthesized in house. Finally, we established a screening paradigm using a pair of cell lines with low and high ALDH1A1 expression, respectively, to uncover novel cell-active ALDH1A1-specific inhibitors from a collection of over 1,000 small molecules.

Drug-induced QT prolongation leads to life-threatening cardiotoxicity, mostly through blockage of the human ether-à-go-go-related gene (hERG) encoded potassium ion (K(+) ) channels. The hERG channel is one of the most important antitargets to be addressed in the early stage of drug discovery process, in order to avoid more costly failures in the development phase. Using a thallium flux assay, 4,323 molecules were screened for hERG channel inhibition in a quantitative high throughput screening (qHTS) format. Here, we present support vector classification (SVC) models of hERG channel inhibition with the averaged area under the receiver operator characteristics curve (AUC-ROC) of 0.93 for the tested compounds. Both Jackknifing and bootstrapping have been employed to rebalance the heavily biased training datasets, and the impact of these two under-sampling rebalance methods on the performance of the predictive models is discussed. Our results indicated that the rebalancing techniques did not enhance the predictive power of the resulting models; instead, adoption of optimal cutoffs could restore the desirable balance of sensitivity and specificity of the binary classifiers. In an external validation set of 66 drug molecules, the SVC model exhibited an AUC-ROC of 0.86, further demonstrating the utility of this modeling approach to predict hERG liabilities.

The endoplasmic reticulum (ER) is involved in Ca2+ signaling and protein folding. ER Ca2+ depletion and accumulation of unfolded proteins activate the molecular chaperone GRP78 (glucose-regulated protein 78) which in turn triggers the ER stress response (ERSR) pathway aimed to restore ER homeostasis. Failure to adapt to stress, however, results in apoptosis. We and others have shown that malignant cells are more susceptible to ERSR-induced apoptosis than their normal counterparts, implicating the ERSR as a potential target for cancer therapeutics. Predicated on these findings, we developed an assay that uses a GRP78 biosensor to identify small molecule activators of ERSR in glioma cells. We performed a quantitative high-throughput screen (qHTS) against a collection of ~425,000 compounds and a comprehensive panel of orthogonal secondary assays was formulated for stringent compound validation. We identified novel activators of ERSR, including a compound with a 2,9-diazaspiro[5.5]undecane core, which depletes intracellular Ca2+ stores and induces apoptosis-mediated cell death in several cancer cell lines, including patient-derived and 3D cultures of glioma cells. This study demonstrates that our screening platform enables the identification and profiling of ERSR inducers with cytotoxic activity and advocates for characterization of these compound in in vivo models.

The structure-activity relationship (SAR) study of two chemotypes identified as inhibitors of the human NAD(+)-dependent 15-hydroxyprostaglandin dehydrogenase (HPGD, 15-PGDH) was conducted. Top compounds from both series displayed potent inhibition (IC50 <50 nM), demonstrate excellent selectivity towards HPGD and potently induce PGE2 production in A549 lung cancer and LNCaP prostate cancer cells.

Small-molecule compounds approved for use as drugs may be "repurposed" for new indications and studied to determine the mechanisms of their beneficial and adverse effects. A comprehensive collection of all small-molecule drugs approved for human use would be invaluable for systematic repurposing across human diseases, particularly for rare and neglected diseases, for which the cost and time required for development of a new chemical entity are often prohibitive. Previous efforts to build such a comprehensive collection have been limited by the complexities, redundancies, and semantic inconsistencies of drug naming within and among regulatory agencies worldwide; a lack of clear conceptualization of what constitutes a drug; and a lack of access to physical samples. We report here the creation of a definitive, complete, and nonredundant list of all approved molecular entities as a freely available electronic resource and a physical collection of small molecules amenable to high-throughput screening.

Chromo/fluorophoric properties often accompany the heterocyclic scaffolds and impurities that comprise libraries used for high-throughput screening (HTS). These properties affect assay outputs obtained with optical detection, thus complicating analysis and leading to false positives and negatives. Here, we report the fluorescence profile of more than 70,000 samples across spectral regions commonly utilized in HTS. The quantitative HTS paradigm was utilized to test each sample at seven or more concentrations over a 4-log range in 1,536-well format. Raw fluorescence was compared with fluorophore standards to compute a normalized response as a function of concentration and spectral region. More than 5% of library members were brighter than the equivalent of 10 nM 4-methyl umbelliferone, a common UV-active probe. Red-shifting the spectral window by as little as 100 nm was accompanied by a dramatic decrease in autofluorescence. Native compound fluorescence, fluorescent impurities, novel fluorescent compounds, and the utilization of fluorescence profiling data are discussed.