Recent Publications:

Nature Biotechnology	Comprehensive characterization of cytochrome P450 isozyme selectivity across chemical libraries [ PDF ] Veith H, Southall N, Huang R, James T, Fayne D, Artemenko N, Shen M, Inglese J, Austin CP, Lloyd DG, Auld DS The cytochrome P450 (CYP) gene family catalyzes drug metabolism and bioactivation and is therefore relevant to drug development. We determined potency values for 17,143 compounds against five recombinant CYP isozymes (1A2, 2C9, 2C19, 2D6 and 3A4) using an in vitro bioluminescent assay. The compounds included libraries of US Food and Drug Administration (FDA)-approved drugs and screening libraries. We observed cross-library isozyme inhibition (30–78%) with important differences between libraries. Whereas only 7% of the typical screening library was inactive against all five isozymes, 33% of FDA-approved drugs were inactive, reflecting the optimized pharmacological properties of the latter. Our results suggest that low CYP 2C isozyme activity is a common property of drugs, whereas other isozymes, such as CYP 2D6, show little discrimination between drugs and unoptimized compounds found in screening libraries. We also identified chemical substructures that differentiated between the five isozymes. The pharmacological compendium described here should further the understanding of CYP isozymes.
Current Topics in Medicinal Chemistry	The Pilot Phase of the NIH Chemical Genomics Center. [ PDF ] Thomas CJ, Auld DS, Huang R, Huang W, Jadhav A, Johnson RL, Leister W, Maloney DJ, Marugan JJ, Michael S, Simeonov A, Southall N, Xia M, Zheng W, Inglese J, Austin CP. The NIH Chemical Genomics Center (NCGC) was the inaugural center of the Molecular Libraries and Screening Center Network (MLSCN). Along with the nine other research centers of the MLSCN, the NCGC was established with a primary goal of bringing industrial technology and experience to empower the scientific community with small molecule compounds for use in their research. We intend this review to serve as 1) an introduction to the NCGC standard operating procedures, 2) an overview of several of the lessons learned during the pilot phase and 3) a review of several of the innovative discoveries reported during the pilot phase of the MLSCN.
Toxicological Sciences	Weighted Feature Significance (WFS): A Simple, Interpretable Model of Compound Toxicity Based on the Statistical Enrichment of Structural Features. [ PDF ] Huang R, Southall N, Xia M, Cho MH, Jadhav A, Nguyen DT, Inglese J, Tice RR, Austin CP In support of the U.S. Tox21 program, we have developed a simple and chemically intuitive model we call Weighted Feature Significance (WFS) to predict the toxicological activity of compounds, based on the statistical enrichment of structural features in toxic compounds. We trained and tested the model on: (1) data from quantitative high-throughput screening (qHTS) cytotoxicity and caspase activation assays conducted at the NIH Chemical Genomics Center (NCGC), (2) data from Salmonella typhimurium reverse mutagenicity assays conducted by the U.S. National Toxicology Program (NTP), and (3) hepatotoxicity data published in the Registry of Toxic Effects of Chemical Substances (RTECS). Enrichments of structural features in toxic compounds are evaluated for their statistical significance and compiled into a simple additive model of toxicity, and then used to score new compounds for potential toxicity. The predictive power of the model for cytotoxicity was validated using an independent set of compounds from the U.S. Environmental Protection Agency (EPA) tested also at the NCGC. We compared the performance of our WFS approach with classical classification methods such as Naïve Bayesian clustering and support vector machines. In most test cases, WFS showed similar or slightly better predictive power, especially in the prediction of hepatotoxic compounds, where WFS appeared to have the best performance among the three methods. The new algorithm has the important advantages of simplicity, power, interpretability, and ease of implementation.
Molecular Biosystems	A quantitative high-throughput screen for modulators of IL-6 signaling: a model for interrogating biological networks using chemical libraries. [ PDF ] Johnson RL, Huang R, Jadhav A, Southall N, Wichterman J, MacArthur R, Xia M, Bi K, Printen J, Austin CP, Inglese J. Small molecule modulators are critical for dissecting and understanding signaling pathways at the molecular level. Interleukin 6 (IL-6) is a cytokine that signals via the JAK-STAT pathway and is implicated in cancer and inflammation. To identify modulators of this pathway, we screened a chemical collection against an IL-6 responsive cell line stably expressing a beta-lactamase reporter gene fused to a sis-inducible element (SIE-bla cells). This assay was optimized for a 1536-well microplate format and screened against 11 693 small molecules using quantitative high-throughput screening (qHTS), a method that assays a chemical library at multiple concentrations to generate titration-response profiles for each compound. The qHTS recovered 564 actives with well-fit curves that clustered into 32 distinct chemical series of 13 activators and 19 inhibitors. A retrospective analysis of the qHTS data indicated that single concentration data at 1.5 and 7.7 microM scored 35 and 71% of qHTS actives, respectively, as inactive and were therefore false negatives. Following counter screens to identify fluorescent and non-selective series, we found four activator and one inhibitor series that modulated SIE-bla cells but did not show similar activity in reporter gene assays induced by EGF and hypoxia. Small molecules within these series will make useful tool compounds to investigate IL-6 signaling mediated by JAK-STAT activation.
Analytical Biochemistry	A new homogeneous high-throughput screening assay for profiling compound activity on the human ether-a-go-go-related gene channel. [ PDF ] Titus SA, Beacham D, Shahane SA, Southall N, Xia M, Huang R, Hooten E, Zhao Y, Shou L, Austin CP, Zheng W. Long QT syndrome, either inherited or acquired from drug treatments, can result in ventricular arrhythmia (torsade de pointes) and sudden death. Human ether-a-go-go-related gene (hERG) channel inhibition by drugs is now recognized as a common reason for the acquired form of long QT syndrome. It has been reported that more than 100 known drugs inhibit the activity of the hERG channel. Since 1997, several drugs have been withdrawn from the market due to the long QT syndrome caused by hERG inhibition. Food and Drug Administration regulations now require safety data on hERG channels for investigative new drug (IND) applications. The assessment of compound activity on the hERG channel has now become an important part of the safety evaluation in the process of drug discovery. During the past decade, several in vitro assay methods have been developed and significant resources have been used to characterize hERG channel activities. However, evaluation of compound activities on hERG have not been performed for large compound collections due to technical difficulty, lack of throughput, and/or lack of biological relevance to function. Here we report a modified form of the FluxOR thallium flux assay, capable of measuring hERG activity in a homogeneous 1536-well plate format. To validate the assay, we screened a 7-point dilution series of the LOPAC 1280 library collection and reported rank order potencies of ten common hERG inhibitors. A correlation was also observed for the hERG channel activities of 10 known hERG inhibitors determined in this thallium flux assay and in the patch clamp experiment. Our findings indicate that this thallium flux assay can be used as an alternative method to profile large-volume compound libraries for compound activity on the hERG channel.
Biochemistry	The Identification of Aminothienopyridazine Inhibitors of Tau Assembly by Quantitative High-Throughput Screening. Crowe A, Huang W, Ballatore C, Johnson RL, Hogan AM, Huang R, Wichtermann J, McCoy J, Huryn DM, Auld DS, Smith Iii AB, Inglese J, Trojanowski JQ, Austin CP, Brunden KR, Lee VM. Inclusions comprised of fibrils of the microtubule (MT)-associated protein tau are found in the brains of those with Alzheimer's disease (AD) and other neurodegenerative tauopathies. The pathology that is observed in these diseases is believed to result from the formation of toxic tau oligomers or fibrils, and/or from the loss of normal tau function due to its sequestration into insoluble deposits. Hence, small molecules that prevent tau oligomerization and/or fibrillization might have therapeutic value. Indeed, examples of such compounds have been published but nearly all have properties that render them unsuitable as drug candidates. For these reasons, we conducted quantitative high-throughput screening (qHTS) of ~292,000 compounds to identify drug-like inhibitors of tau assembly. The fibrillization of a truncated tau fragment that contains four MT-binding domains was monitored in an assay that employed complementary thioflavine T fluorescence and fluorescence polarization methods. Previously described classes of inhibitors as well as new scaffolds were identified, including novel aminothienopyridazines (ATPZ's). A number of ATPZ analogs were synthesized and structure-activity relationships were defined. Further characterization of representative ATPZ compounds showed they do not interfere with tau-mediated MT assembly, and they are significantly more effective at preventing the fibrillization of tau than the Abeta(1-42) peptide which forms AD senile plaques. Thus, the ATPZ molecules described here represent a novel class of tau assembly inhibitors that merit further development for testing in animal models of AD-like tau pathology.
Journal of Biomolecular Screening	A New Homogeneous High Throughput Screening Assay for Profiling Compound Activity on the hERG Channel. Titus S, Beacham D, Shahane S, Southall N, Xia M, Huang R, Hooten E, Zhao Y, Shou L, Austin CP, Zheng W. Long QT syndrome, either inherited or acquired from drug treatments, can result in ventricular arrhythmia (Torsade de Pointes) and sudden death. hERG channel inhibition by drugs is now recognized as a common reason for the acquired form of Long QT syndrome. It has been reported that over 100 known drugs inhibit the activity of the hERG channel. Since 1997, several drugs have been withdrawn from the market due to the long QT syndrome caused by hERG inhibition. FDA regulations now require safety data on hERG channels for investigative new drug (IND) applications. The assessment of compound activity on the hERG channel has now become an important part of the safety evaluation in the process of drug discovery. In the last decade, several in vitro assay methods have been developed, and significant resources have been used to characterize hERG channel activities. However, evaluation of compound activities on hERG have not been performed for large compound collections due to technical difficulty, lack of throughput, or lack of biological relevance to function. We report here a modified form of the FluxOR(TM) thallium flux assay, capable of measuring hERG activity in a homogeneous 1536-well plate format. To validate the assay, we screened a 7 point dilution series of the LOPAC1280 collection and reported rank order potencies of eight common hERG inhibitors. A correlation was also observed for the hERG channel activities of 10 known hERG inhibitors determined in this thallium flux assay and in the patch clamp experiment. Our findings indicate that this thallium flux assay can be used as an alternative method to profile large volume compound libraries for compound activity on the hERG channel.
Toxicological Sciences	Identification of Chemical Compounds that Induce HIF-1{alpha} Activity. Xia M, Huang R, Sun Y, Semenza GL, Aldred SF, Witt KL, Inglese J, Tice RR, Austin CP. Cellular metabolism depends on the availability of oxygen and the major regulator of oxygen homeostasis is hypoxia-inducible factor 1 (HIF-1), a highly conserved transcription factor that plays an essential role in cellular and systemic homeostatic responses to hypoxia. HIF-1 is a heterodimeric transcription factor composed of hypoxia-inducible HIF-1alpha and constitutively-expressed HIF-1beta. Under hypoxic conditions, the two subunits dimerize, allowing translocation of the HIF-1 complex to the nucleus where it binds to hypoxia response elements (HRE) and activates expression of target genes implicated in angiogenesis, cell growth, and survival. The HIF-1 pathway is essential to normal growth and development, and is involved in the pathophysiology of cancer, inflammation, and ischemia. Thus, there is considerable interest in identifying compounds that modulate the HIF-1 signaling pathway. To assess the ability of environmental chemicals to stimulate the HIF-1 signaling pathway, we screened a National Toxicology Program collection of 1408 compounds using a cell-based beta-lactamase HRE reporter gene assay in a quantitative high throughput screening (qHTS) format. Twelve active compounds were identified. These compounds were tested in a confirmatory assay for induction of vascular endothelial growth factor, a known hypoxia target gene, and confirmed compounds were further tested for their ability to mimic the effect of a reduced-oxygen environment on hypoxia-regulated promoter activity. Based on this testing strategy, three compounds (o-phenanthroline, iodochlorohydroxyquinoline, cobalt sulfate heptahydrate) were confirmed as hypoxia mimetics, while 2 compounds (7-diethylamino-4-methylcoumarin and 7,12-dimethylbenz(a)anthracence) were found to interact with HIF-1 in a manner different from hypoxia. These results demonstrate the effectiveness of qHTS in combination with secondary assays for identification of HIF-1alpha inducers and for distinguishing among inducers based on their pattern of activated hypoxic target genes. Identification of environmental compounds having HIF-1alpha activation activity in cell-based assays may be useful for prioritizing chemicals for further testing as hypoxia-response inducers in vivo.
Bioorganic & Medicinal Chemistry Letters	Exploration and Optimization of Substituted Triazolothiadiazines and Triazolopyridazines as PDE4 Inhibitors. Skoumbourdis AP, et al. An expansion of structure-activity studies on a series of substituted 7H-[1,2,4]triazolo[3,4-b][1,3,4]thiadiazine PDE4 inhibitors and the introduction of a related [1,2,4]triazolo[4,3-b]pyridazine based inhibitor of PDE4 is presented. The development of SAR included strategic incorporation of known substituents on the critical catachol diether moiety of the 6-phenyl appendage on each heterocyclic core. From these studies, (R)-3-(2,5-dimethoxyphenyl)-6-(4-methoxy-3-(tetrahydrofuran-3-yloxy)phenyl)-7H-[1,2,4]triazolo[3,4-b][1,3,4]thiadiazine (10) and (R)-3-(2,5-dimethoxyphenyl)-6-(4-methoxy-3-(tetrahydrofuran-3-yloxy)phenyl)-[1,2,4]triazolo[4,3-b]pyridazine (18) were identified as highly potent PDE4A inhibitors. Each of these analogues was submitted across a panel of 21 PDE family members and was shown to be highly selective for PDE4 isoforms (PDE4A, PDE4B, PDE4C, PDE4D). Both 10 and 18 were then evaluated in divergent cell-based assays to assess their relevant use as probes of PDE4 activity. Finally, docking studies with selective ligands (including 10 and 18) were undertaken to better understand this chemotypes ability to bind and inhibit PDE4 selectively.
Proceedings of the National Academy of Sciences	Identification of compounds that potentiate CREB signaling as possible enhancers of long-term memory. Xia M, Huang R, Guo V, Southall N, Cho MH, Inglese J, Austin CP, Nirenberg M. Many studies have implicated the cAMP Response Element Binding (CREB) protein signaling pathway in long-term memory. To identify small molecule enhancers of CREB activation of gene expression, we screened approximately 73,000 compounds, each at 7-15 concentrations in a quantitative high-throughput screening (qHTS) format, for activity in cells by assaying CREB mediated beta-lactamase reporter gene expression. We identified 1,800 compounds that potentiated CREB mediated gene expression, with potencies as low as 16 nM, comprising 96 structural series. Mechanisms of action were systematically determined, and compounds that affect phosphodiesterase 4, protein kinase A, and cAMP production were identified, as well as compounds that affect CREB signaling via apparently unidentified mechanisms. qHTS folowed by interrogation of pathway targets is an efficient paradigm for lead generation for chemical genomics and drug development.
Bioorganic & Medicinal Chemistry Letters	Identification of a potent new chemotype for the selective inhibition of PDE4. Skoumbourdis AP, Huang R, Southall N, Leister W, Guo V, Cho MH, Inglese J, Nirenberg M, Austin CP, Xia M, Thomas CJ. A series of substituted 3,6-diphenyl-7H-[1,2,4]triazolo[3,4-b][1,3,4]thiadiazines were prepared and analyzed as inhibitors of phosphodiesterase 4 (PDE4). Synthesis, structure–activity relationships, and the selectivity of a highly potent analogue against related phosphodiesterase isoforms are presented.
Chemical Research Toxicology	Characterization of Diversity in Toxicity Mechanism Using In Vitro Cytotoxicity Assays in Quantitative High Throughput Screening. Huang R, Southall N, Cho MH, Xia M, Inglese J, Austin CP. Assessing the potential health risks of environmental chemical compounds is an expensive undertaking which has motivated the development of new alternatives to traditional in vivo toxicological testing. One approach is to stage the evaluation, beginning with less expensive and higher throughput in vitro testing before progressing to more definitive trials. In vitro testing can be used to generate a hypothesis about a compound's mechanism of action, which can then be used to design an appropriate in vivo experiment. Here we begin to address the question of how to design such a battery of in vitro cell-based assays by combining data from two different types of assays, cell viability and caspase activation, with the aim of elucidating mechanism of action. Because caspase activation is a transient event during apoptosis, it is not possible to design a single end-point assay protocol that would identify all instances of compound-induced caspase activation. Nevertheless, useful information about compound mechanism of action can be obtained from these assays in combination with cell viability data. Unsupervised clustering in combination with Dunn's cluster validity index is a robust method for identifying mechanisms of action without requiring any a priori knowledge about mechanisms of toxicity. The performance of this clustering method is evaluated by comparing the clustering results against literature annotations of compound mechanisms.
Toxicology in Vitro	A bioluminescent cytotoxicity assay for assessment of membrane integrity using a proteolytic biomarker. Cho MH, Niles A, Huang R, Inglese J, Austin CP, Riss T, Xia M. Measurement of cell membrane integrity has been widely used to assess chemical cytotoxity. Several assays are available for determining cell membrane integrity including differential labeling techniques using neutral red and trypan blue dyes or fluorescent compounds such as propidium iodide. Other common methods for assessing cytotoxicity are enzymatic "release" assays which measure the extra-cellular activities of lactate dehydrogenase (LDH), adenylate kinase (AK), or glyceraldehyde-3-phosphate dehydrogenase (GAPDH) in culture medium. However, all these assays suffer from several practical limitations, including multiple reagent additions, scalability, low sensitivity, poor linearity, or requisite washes and medium exchanges. We have developed a new cytotoxicity assay which measures the activity of released intracellular proteases as a result of cell membrane impairment. It allows for a homogenous, one-step addition assay with a luminescent readout. We have optimized and miniaturized this assay into a 1536-well format, and validated it by screening a library of known compounds from the National Toxicology Program (NTP) using HEK 293 and human renal mesangial cells by quantitative high-throughput screening (qHTS). Several known and novel membrane disrupters were identified from the library, which indicates that the assay is robust and suitable for large scale library screening. This cytotoxicity assay, combined with the qHTS platform, allowed us to quickly and efficiently evaluate compound toxicities related to cell membrane integrity.
Journal of Medicinal Chemistry	Fluorescence Spectroscopic Profiling of Compound Libraries. Simeonov S, Jadhav A, Thomas CJ, Wang Y, Huang R, Southall N, Shinn P, Smith J, Austin CP, Inglese J. Chromo/fluorophoric properties often accompany the conjugated, aromatic and heterocyclic features of many of the scaffolds and impurities that make up library samples used for high throughput screening (HTS). These properties impart highly variable effects on assay outputs employing optical detection, thus complicating the interpretation of data and leading to false positives and negatives. Here, we report the comprehensive fluorescence profile of >70,000 samples across multiple spectral regions commonly utilized in HTS assays. The quantitative HTS (qHTS) paradigm was utilized to test each sample at seven or more concentration points over a 4-log concentration range in 1536-well format, with raw fluorescence response collected using a CCD-based imager. The resulting output was compared with fluorophore standards to compute a normalized fluorescence response (termed fluorophore-equivalent concentration, FEC) for each sample, concentration, and relevant spectral region. The greatest fraction of fluorescent compounds appeared in the UV-end of the light spectrum, where over 5% of library members matched or exceeded 10 nM FEC of 4-methylumbelliferone and AlexaFluor 350, while approximately 1.8% of the library matched or exceeded 100 nM FEC of these standards. Red-shifting the spectral window by as little as 100 nm was accompanied by a dramatic decrease in autofluorescence. Native compound fluorescence, scaffold overlap with known fluorophores, fluorescent impurities, novel fluorescent compounds, and the ability to discriminate generalities of fluorescent interferences and devise strategies to identify them are discussed.
Environmental Health Perspectives	Compound Cytotoxicity Profiling Using Quantitative High-Throughput Screening. Xia M, Huang R, Witt KL, Southall N, Fostel J, Cho MH, Jadhav A, Smith CS, Inglese J, Portier CJ, Tice RR, Austin CP. Background: The propensity of compounds to produce adverse health effects in humans is generally evaluated using animal-based test methods. Such methods can be relatively expensive, low-throughput, and associated with pain suffered by the treated animals. In addition, differences in species biology may confound extrapolation to human health effects. Objective: The U.S. National Toxicology Program and the NIH Chemical Genomics Center are collaborating to identify a battery of cell-based screens to prioritize compounds for further toxicological evaluation. Methods: 1,408 compounds previously tested in one or more traditional toxicological assays were profiled for cytotoxicity using quantitative high-throughput screening (qHTS) in 13 human and rodent cell types derived from six common targets of xenobiotic toxicity (liver, blood, kidney, nerve, lung, skin). Selected cytotoxicants were further tested to define response kinetics. Results: qHTS of these compounds produced robust and reproducible results which allowed cross-compound, cross-cell type, and cross-species comparisons. Some compounds were cytotoxic to all cell types at similar concentrations, while others exhibited species- or cell typespecific cytotoxicity. Closely related cell types and analogous cell types in human and rodent frequently showed different patterns of cytotoxicity. Some compounds inducing similar levels of cytotoxicity showed distinct time-dependence in kinetic studies, consistent with known mechanisms of toxicity. Conclusions: The generation of high-quality cytotoxicity data on this large library of known compounds using qHTS demonstrates the potential of this methodology to profile a much broader array of assays and compounds, which, in aggregate, may be valuable for prioritizing compounds for further toxicological evaluation, identifying compounds with particular mechanisms of action, and potentially predicting in vivo biological response
Molecular Cancer Therapeutics	Anticancer medicines in development: assessment of bioactivity profiles within the National Cancer Institute anticancer screening data. Covell DG, Huang R, Wallqvist A. We present an analysis of current anticancer compounds that are in phase I, II, or III clinical trials and their structural analogues that have been screened in the National Cancer Institute (NCI) anticancer screening program. Bioactivity profiles, measured across the NCI 60 cell lines, were examined for a correspondence between the type of cancer proposed for clinical testing and selective sensitivity to appropriately matched tumor subpanels in the NCI screen. These results find strongest support for using the NCI anticancer screen to select analogue compounds with selective sensitivity to the leukemia, colon, central nervous system, melanoma, and ovarian panels, but not for renal, prostate, and breast panels. These results are extended to applications of two-dimensional structural features to further refine compound selections based on tumor panel sensitivity obtained from tumor screening results.
Journal of Chemical Information and Modeling	Chemoinformatic analysis of NCI preclinical tumor data: evaluating compound efficacy from mouse xenograft data, NCI-60 screening data, and compound descriptors Wallqvist A, Huang R, Covell DG. We provide a chemoinformatic examination of the NCI public human tumor xenograft data to explore relationships between small molecules, treatment modality, efficacy, and toxicity. Efficacy endpoints of tumor weight reduction (TW) and survival time increase (ST) compared to tumor bearing control mice were augmented by a toxicity measure, defined as the survival advantage of treated versus control animals (TX). These endpoints were used to define two independent therapeutic indices (TIs) as the ratio of efficacy (TW or ST) to toxicity (TX). Linear models predictive of xenograft endpoints were successfully constructed (0.67 < r(2) < or = 0.74)(observed_versus_predicted) using a model comprised of variables in treatment modality, chemoinformatic descriptors, and in vitro cell growth inhibition in the NCI 60-cell assay. Cross-validation analysis based on randomly chosen training subsets found these predictive correlations to be robust. Model-based sensitivity analysis found chemistry and growth inhibition to provide the best, and treatment modality the worst, indicators of xenograft endpoint. The poor predictive power derived from treatment alone appears to be of less importance to xenograft outcome for compounds having strongly similar chemical and biological features. ROC-based model validation found a 70% positive predictive value for distinguishing FDA approved oncology agents from available xenograft tested compounds. Additional chemoinformatic applications are provided that relate xenograft outcome to biological pathways and putative mechanism of compound action. These results find a strong relationship between xenograft efficacy and pathways comprised of genes having highly correlated mRNA expressions. Our analysis demonstrates that chemoinformatic studies utilizing a combination of xenograft data and in vitro preclinical testing offer an effective means to identify compound classes with superior efficacy and reduced toxicity.
Molecular Cancer Therapeutics	Targeting changes in cancer: assessing pathway stability by comparing pathway gene expression coherence levels in tumor and normal tissues. Huang R, Wallqvist A, Covell DG. The purpose of this study is to examine gene expression changes occurring in cancer from a pathway perspective by analyzing the level of pathway coherence in tumor tissues in comparison with their normal counterparts. Instability in pathway regulation patterns can be considered either as a result of or as a contributing factor to genetic instability and possibly cancer. Our analysis has identified pathways that show a significant change in their coherence level in tumor tissues, some of which are tumor type specific, indicating novel targets for cancer type-specific therapies. Pathways are found to have a general tendency to lose their gene expression coherence in tumor tissues when compared with normal tissues, especially for signaling pathways. The selective growth advantage of cancer cells over normal cells seems to originate from their preserved control over vital pathways to ensure survival and altered signaling, allowing excessive proliferation. We have additionally investigated the tissue-related instability of pathways, providing valuable clues to the cellular processes underlying the tumorigenesis and/or growth of specific cancer types. Pathways that contain known cancer genes (i.e., "cancer pathways") show significantly greater instability and are more likely to become incoherent in tumor tissues. Finally, we have proposed strategies to target instability (i.e., pathways that are prone to changes) by identifying compound groups that show selective activity against pathways with a detectable coherence change in cancer. These results can serve as guidelines for selecting novel agents that have the potential to specifically target a particular pathway that has relevance in cancer.
Journal of Medicinal Chemistry	Targeting changes in cancer: assessing pathway stability by comparing pathway gene expression coherence levels in tumor and normal tissues. Huang R, Wallqvist A, Covell DG. This paper examines two biological models of anticancer activity, cytotoxicity and hollow fiber (HF) activity, for chemotherapeutic agents evaluated as part of the National Cancer Institute's (NCI's) drug screening effort. Our analysis proposes strategies to globally assess compounds tested in the NCI's 60-cell (NCI60) in vitro anticancer screen in terms of structural features, biological activity, target specificity, and mechanism of action by data integration via our self-organizing maps of structural and biological response patterns. We have built statistical models to predict compound potency and HF activity based on physicochemical properties. Our results find that it is the combination of different structural properties that determines a compound's biological activity. A direct correlation is also found between compound potency and specificity, indicating that specific targeting, rather than promiscuous poisoning, gives rise to potency. Finally, we offer a strategy to exploit this relationship for future mining of novel anticancer candidates.
Journal of Chemical Information and Modeling	Evaluating Chemical Structure Similarity as an Indicator of Cellular Growth Inhibition Wallqvist A, Huang R, Thanki N, Covell DG. Chemical variations of small compounds are commonly used to probe biological systems and potentially discover lead-like compounds with selective target activity. Molecular probes are either generated by synthesis or acquired through directed searches of commercially available compound libraries. The data generated when testing the probes in various biological systems constitutes a structure/activity analysis. The ability to detect variations and classify biological responses requires the analysis of a compound in multiple assays. While the concept of a structure/activity relationship is straightforward, its implementation can vary considerably depending on the biological system under study and the probe library selected for testing. The analysis presented here will focus on the accumulated compound library used to screen for growth inhibition across the National Cancer Institute's panel of 60 tumor cells. The considerable chemical and biological diversity inherent in these data offers an opportunity to establish a quantifiable connection between chemical structure and biological activity. We find that the connection between structure and biological response is not symmetric, with biological response better at predicting chemical structure than vice versa. Structurally and functionally similar compounds can have distinguishable biological responses reflecting different mechanisms of action.
Genomics	Comprehensive analysis of pathway or functionally related gene expression in the National Cancer Institute's anticancer screen. Huang R, Wallqvist A, Covell DG. We have analyzed the level of gene coregulation, using gene expression patterns measured across the National Cancer Institute's 60 tumor cell panels (NCI(60)), in the context of predefined pathways or functional categories annotated by KEGG (Kyoto Encyclopedia of Genes and Genomes), BioCarta, and GO (Gene Ontology). Statistical methods were used to evaluate the level of gene expression coherence (coordinated expression) by comparing intra- and interpathway gene-gene correlations. Our results show that gene expression in pathways, or groups of functionally related genes, has a significantly higher level of coherence than that of a randomly selected set of genes. Transcriptional-level gene regulation appears to be on a "need to be" basis, such that pathways comprising genes encoding closely interacting proteins and pathways responsible for vital cellular processes or processes that are related to growth or proliferation, specifically in cancer cells, such as those engaged in genetic information processing, cell cycle, energy metabolism, and nucleotide metabolism, tend to be more modular (lower degree of gene sharing) and to have genes significantly more coherently expressed than most signaling and regular metabolic pathways. Hierarchical clustering of pathways based on their differential gene expression in the NCI(60) further revealed interesting interpathway communications or interactions indicative of a higher level of pathway regulation. The knowledge of the nature of gene expression regulation and biological pathways can be applied to understanding the mechanism by which small drug molecules interfere with biological systems.
Molecular Cancer Therapeutics	Drugs aimed at targeting characteristic karyotypic phenotypes of cancer cells. Wallqvist A, Huang R, Covell DG, Roschke AV, Gelhaus KS, Kirsch IR. The karyotypic features of cancer cells have not been a particular focus of anticancer drug targeting either as guidance for treatment or as specific drug targets themselves. Cancer cell lines typically have considerable, characteristic, and variable chromosomal aberrations. Here, we consider small-molecule screening data across the National Cancer Institute's 60 tumor cell line drug screening panel (NCI-60) analyzed for specific association with karyotypic variables (numerical and structural complexity and heterogeneity) determined for these same cell lines. This analysis is carried out with the aid of a self-organizing map allowing for a simultaneous assessment of all screened compounds, revealing an association between karyotypic variables and a unique part of the cytotoxic response space. Thirteen groups of compounds based on related specific chemical structural motifs are identified as possible leads for anticancer drug discovery. These compounds form distinct groups of molecules associated with relatively unexplored regions of the NCI-60 self-organizing map where anticancer agents currently standard in the clinic are not present. We suggest that compounds identified in this study may represent new classes of potential anticancer agents.
Proteins	Linking tumor cell cytotoxicity to mechanism of drug action: an integrated analysis of gene expression, small-molecule screening and structural databases. Covell DG, Wallqvist A, Huang R, Thanki N, Rabow AA, Lu XJ. An integrated, bioinformatic analysis of three databases comprising tumor-cell-based small molecule screening data, gene expression measurements, and PDB (Protein Data Bank) ligand-target structures has been developed for probing mechanism of drug action (MOA). Clustering analysis of GI50 profiles for the NCI's database of compounds screened across a panel of tumor cells (NCI60) was used to select a subset of unique cytotoxic responses for about 4000 small molecules. Drug-gene-PDB relationships for this test set were examined by correlative analysis of cytotoxic response and differential gene expression profiles within the NCI60 and structural comparisons with known ligand-target crystallographic complexes. A survey of molecular features within these compounds finds thirteen conserved Compound Classes, each class exhibiting chemical features important for interactions with a variety of biological targets. Protein targets for an additional twelve Compound Classes could be directly assigned using drug-protein interactions observed in the crystallographic database. Results from the analysis of constitutive gene expressions established a clear connection between chemo-resistance and overexpression of gene families associated with the extracellular matrix, cytoskeletal organization, and xenobiotic metabolism. Conversely, chemo-sensitivity implicated overexpression of gene families involved in homeostatic functions of nucleic acid repair, aryl hydrocarbon metabolism, heat shock response, proteasome degradation and apoptosis. Correlations between chemo-responsiveness and differential gene expressions identified chemotypes with nonselective (i.e., many) molecular targets from those likely to have selective (i.e., few) molecular targets. Applications of data mining strategies that jointly utilize tumor cell screening, genomic, and structural data are presented for hypotheses generation and identifying novel anticancer candidates.
The Pharmacogenomics Journal	Linking pathway gene expressions to the growth inhibition response from the National Cancer Institute's anticancer screen and drug mechanism of action. Huang R, Wallqvist A, Thanki N, Covell DG. Novel strategies are proposed to quantitatively analyze and relate biological pathways to drug responses using gene expression and small-molecule growth inhibition data (GI(50)) derived from the National Cancer Institute's 60 cancer cells (NCI(60)). We have annotated groups of drug GI(50) responses with pathways defined by the Kyoto Encyclopedia of Genes and Genomes (KEGG) and BioCarta, and functional categories defined by Gene Ontology (GO), through correlations between pathway gene expression patterns and drug GI(50) profiles. Drug-gene-pathway relationships may then be utilized to find drug targets or target-specific drugs. Significantly correlated pathways and the gene products involved represent interesting targets for further exploration, whereas drugs that are significantly correlated with only certain pathways are more likely to be target specific. Separate pathway clustering finds that pathways engaged in the same biological process tend to have similar drug correlation patterns. The biological and statistical significances of our method are established by comparison to known small-molecule inhibitor-gene target relationships reported in the literature and by standard randomization procedures. The results of our pathway, gene expression and drug-induced growth inhibition associations, can serve as a basis for proposing testable hypotheses about potential anticancer drugs, their targets, and mechanisms of action.
Biochemical Pharmacology	Anticancer metal compounds in NCI's tumor-screening database: putative mode of action. Huang R, Wallqvist A, Covell DG Clustering analysis of tumor cell cytotoxicity profiles for the National Cancer Institute (NCI)'s open compound repository has been used to catalog over 1100 metal or metalloid containing compounds with potential anticancer activity. The molecular features and corresponding reactivity of these compounds have been analyzed in terms of properties of their metals, their associated organic components (ligands) and their capacity to inhibit tumor cell growth. Cytotoxic responses are influenced by both the identity of the metal and the properties of its coordination ligand, with clear associations between structural similarities and cytotoxicity. Assignments of mechanisms of action (MOAs) for these compounds could be segregated into four broad response classes according to preference for binding to biological sulfhydryl groups, chelation, generation of reactive oxygen species (ROS), and production of lipophilic ions. Correlations between specific cytotoxic responses and differential gene expression profiles within the NCI's tumor cell panel serve as a validation for candidate biological targets and putative MOA classes. In addition, specific sensitivity toward subsets of metal containing agents has been found for certain tumor cell panels. Taken together, our results expand the knowledge base available for evaluating, designing and developing new metal-based anticancer drugs that may provide the basis for target-specific therapeutics.
Anticancer Drugs	Cytotoxicity of RH1: NAD(P)H:quinone acceptor oxidoreductase (NQO1)-independent oxidative stress and apoptosis induction. Tudor G, Alley M, Nelson CM, Huang R, Covell DG, Gutierrez P, Sausville EA. The elevated expression of the flavoprotein NAD(P)H:quinone acceptor oxidoreductase (NQO1) (EC 1.6.99.2) in many human solid tumors, along with its ability to activate quinone-based anticancer agents, makes it an excellent target for enzyme-directed drug development. Previous studies have shown a significant statistical correlation between NQO1 enzymatic activity and the cytotoxicity of certain antitumor quinones. RH1 [2,5-diaziridinyl-3-(hydroxymethyl)-6-methyl-1,4-benzoquinone], presently in late preclinical and entering early clinical development, has been previously considered to be an excellent substrate for activation by NQO1. In this study we investigate the cytotoxicity of RH1 in cell lines selected from the NCI's 60 tumor cell line panel, expressing varying levels of NQO1 activity. Exposure time- and concentration-dependent cytotoxicity was seen, apparently independent from levels of NQO1 activity in these cells. Furthermore, the NQO1 inhibitor dicoumarol had no impact on the sensitivity profiles of RH1 response. The HL-60 myeloid leukemia cells, which do not have detectable NQO1 activity, were further investigated. RH1 treatment of HL-60 cells generated high levels of free radicals, which was accompanied by robust redox cycling, oxygen consumption and induction of apoptosis. These results are in agreement with previous data suggesting that, in addition to its activation by NQO1, RH1-induced cytotoxicity might involve alternative pathways for activation of this compound. Furthermore, the high cytotoxicity of RH1 in the leukemia/lymphoma subpanel of the NCI in vitro cell line screen would suggest an empirical rationale for the utilization of this compound in the treatment of these malignancies.