Identification of small molecule inhibitors of anthrax lethal factor.

The virulent spore-forming bacterium Bacillus anthracis secretes anthrax toxin composed of protective antigen (PA), lethal factor (LF) and edema factor (EF). LF is a Zn-dependent metalloprotease that inactivates key signaling molecules, such as mitogen-activated protein kinase kinases (MAPKK), to ultimately cause cell death. We report here the identification of small molecule (nonpeptidic) inhibitors of LF. Using a two-stage screening assay, we determined the LF inhibitory properties of 19 compounds. Here, we describe six inhibitors on the basis of a pharmacophoric relationship determined using X-ray crystallographic data, molecular docking studies and three-dimensional (3D) database mining from the US National Cancer Institute (NCI) chemical repository. Three of these compounds have K(i) values in the 0.5-5 microM range and show competitive inhibition. These molecular scaffolds may be used to develop therapeutically viable inhibitors of LF.

Anticancer activity of BIM-46174, a new inhibitor of the heterotrimeric Galpha/Gbetagamma protein complex.

A large number of hormones and local agonists activating guanine-binding protein-coupled receptors (GPCR) play a major role in cancer progression. Here, we characterize the new imidazo-pyrazine derivative BIM-46174, which acts as a selective inhibitor of heterotrimeric G-protein complex. BIM-46174 prevents the heterotrimeric G-protein signaling linked to several GPCRs mediating (a) cyclic AMP generation (Galphas), (b) calcium release (Galphaq), and (c) cancer cell invasion by Wnt-2 frizzled receptors and high-affinity neurotensin receptors (Galphao/i and Galphaq). BIM-46174 inhibits the growth of a large panel of human cancer cell lines, including anticancer drug-resistant cells. Exposure of cancer cells to BIM-46174 leads to caspase-3-dependent apoptosis and poly(ADP-ribose) polymerase cleavage. National Cancer Institute COMPARE analysis for BIM-46174 supports its novel pharmacologic profile compared with 12,000 anticancer agents. The growth rate of human tumor xenografts in athymic mice is significantly reduced after administration of BIM-46174 combined with either cisplatin, farnesyltransferase inhibitor, or topoisomerase inhibitors. Our data validate the feasibility of targeting heterotrimeric G-protein functions downstream the GPCRs to improve anticancer chemotherapy.

A common pharmacophore for a diverse set of colchicine site inhibitors using a structure-based approach.

Modulating the structure and function of tubulin and microtubules is an important route to anticancer therapeutics, and therefore, small molecules that bind to tubulin and cause mitotic arrest are of immense interest. A large number of synthetic and natural compounds with diverse structures have been shown to bind at the colchicine site, one of the major binding sites on tubulin, and inhibit tubulin assembly. Using the recently determined X-ray structure of the tubulin:colchicinoid complex as the template, we employed docking studies to determine the binding modes of a set of structurally diverse colchicine site inhibitors. These binding models were subsequently used to construct a comprehensive, structure-based pharmacophore that in combination with molecular dynamics simulations confirms and extends our understanding of binding interactions at the colchicine site.

Homology model of the CDK1/cyclin B complex.

We describe a refined homology model of a CDK1/cyclin B complex that was previously used for the structure-based optimization of the Paullone class of inhibitors. The preliminary model was formed from the homologous regions of the deposited CDK2/cyclin A crystal structure. Further refinement of the CDK1/cyclin B complex was accomplished using molecular mechanics and hydropathic analysis with a protocol of constraints and local geometry searches. For the most part, our CKD1/cyclin B homology model is very similar to the high resolution CDK2/cyclin A crystal structure regarding secondary and tertiary features. However, minor discrepancies between the two kinase structures suggest the possibility that ligand design may be specifically tuned for either CDK1 or CDK2. Our examination of the CDK1/cyclin B model includes a comparison with the CDK2/cyclin A crystal structure in the PSTAIRE interface region, connecting portions to the ATP binding domain, as well as the ATP binding site itself.

Epoxide-containing side chains enhance antiproliferative activity of paullones.

The introduction of side chains bearing epoxide motifs into the molecular scaffold of kenpaullone and 9-trifluoromethylpaullone led to improved antiproliferative activity of the novel derivatives for human tumor cell lines. The syntheses were accomplished applying Stille coupling for the introduction of unsaturated side chains into the 2-position of the paullones and subsequently employing a hydrogen peroxide/nitrile mixture for the epoxidation of C,C-double bonds.

Evaluation and comparison of 3D-QSAR CoMSIA models for CDK1, CDK5, and GSK-3 inhibition by paullones.

With a view to the rational design of selective GSK-3beta inhibitors, 3D-QSAR CoMSIA models were developed for the inhibition of the three serine/threonine kinases CDK1/cyclin B, CDK5/p25, and GSK-3beta by compounds from the paullone inhibitor family. The models are based on the kinase inhibition data of 52 paullone entities, which were aligned by a docking routine into the ATP-binding cleft of a CDK1/cyclin B homology model. Variation of grid spacing and column filtering were used during the optimization of the models. The predictive ability of the models was shown by a leave-one-out cross-validation and the prediction of an independent set of test compounds, which were synthesized especially for this purpose. Besides paullones with the basic indolo[3,2-d][1]benzazepine core, the test set comprised novel thieno[3',2':2,3]azepino[4,5-b]indoles, pyrido[2',3':2,3]azepino[4,5-b]indoles, and a pyrido[3',2':4,5]pyrrolo[3,2-d][1]benzazepine. The best statistical values for the CoMSIA were obtained for the CDK1-models (r(2)() = 0.929 and q(2)() = 0.699), which were clearly superior to the models for CDK5 (r(2)() = 0.874 and q(2)() = 0.652) and GSK-3 (r(2)() = 0.871 and q(2)() = 0.554).

Models for evaluating agents intended for the prophylaxis, mitigation and treatment of radiation injuries. Report of an NCI Workshop, December 3-4, 2003.

To develop approaches to prophylaxis/protection, mitigation and treatment of radiation injuries, appropriate models are needed that integrate the complex events that occur in the radiation-exposed organism. While the spectrum of agents in clinical use or preclinical development is limited, new research findings promise improvements in survival after whole-body irradiation and reductions in the risk of adverse effects of radiotherapy. Approaches include agents that act on the initial radiochemical events, agents that prevent or reduce progression of radiation damage, and agents that facilitate recovery from radiation injuries. While the mechanisms of action for most of the agents with known efficacy are yet to be fully determined, many seem to be operating at the tissue, organ or whole animal level as well as the cellular level. Thus research on prophylaxis/protection, mitigation and treatment of radiation injuries will require studies in whole animal models. Discovery, development and delivery of effective radiation modulators will also require collaboration among researchers in diverse fields such as radiation biology, inflammation, physiology, toxicology, immunology, tissue injury, drug development and radiation oncology. Additional investment in training more scientists in radiation biology and in the research portfolio addressing radiological and nuclear terrorism would benefit the general population in case of a radiological terrorism event or a large-scale accidental event as well as benefit patients treated with radiation.

Correlation of nucleoside and nucleobase transporter gene expression with antimetabolite drug cytotoxicity.

Antimetabolite drugs that inhibit nucleic acid metabolism are widely used in cancer chemotherapy. Nucleoside and nucleobase transporters are important for the cellular uptake of nucleic acids and their corresponding anticancer analogue drugs. Thus, these transporters may play a role both in antimetabolite drug sensitivity, by mediating the uptake of nucleoside analogues, and in antimetabolite drug resistance, by mediating the uptake of endogenous nucleosides that may rescue cells from toxicity. Therefore, we examined the relation of the expression of nucleoside and nucleobase transporters to antimetabolite cytotoxicity. We measured the RNA levels of all eight known nucleoside and nucleobase transporters in 50 cell lines included in the National Cancer Institute's Anticancer Drug Screen panel. RNA levels of concentrative nucleoside transporters (CNTs), equilibrative nucleoside transporters (ENTs) and nucleobase transporters (NCBTs) were determined by quantitative RT-PCR using real-time fluorescence acquisition. This method was validated by measuring the expression of the MDR1 gene, and correlating our results with independently determined measurements of MDR1 RNA levels and protein function in these cell lines. We then correlated the pattern of RNA levels to the pattern of cytotoxicity of anticancer drugs in the NCI drug screen database using the COMPARE analysis. Several hypothesized relations between transporter gene expression and cytotoxicity, based upon known interactions between certain nucleoside analogues and transporter proteins, were not observed, suggesting that expression of individual transporters may not be a significant determinant of the cytotoxicity of these drugs. The most closely correlated drug cytotoxicity patterns to transporter gene expression patterns (where increased expression corresponds to increase sensitivity) included those between CNT1 and O6-methylguanine and between ENT2 and hydroxyurea. We also observed that p53 status influenced correlations between ENT1 transporter gene RNA levels and sensitivity to the drugs tiazafurin, AZQ and 3-deazauridine. One of three drugs identified by correlation of cytotoxicity patterns with ENT1 RNA levels, 3-deazauridine, inhibited uptake of the classic ENT1 substrate uridine, demonstrating a physical interaction between an identified drug and the transporter. These studies demonstrate that it is possible to correlate genetic information to functional databases to determine the influence of transport gene expression on drug sensitivity and to identify transporter-drug interactions.

COMPARE: a web accessible tool for investigating mechanisms of cell growth inhibition.

For more than 10 years the National Cancer Institute (NCI) has tested compounds for their ability to inhibit the growth of human tumor cell lines in culture (NCI screen). Work of Ken Paull [J. Natl. Cancer Inst. 81 (1989) 1088] demonstrated that compounds with similar mechanism of cell growth inhibition show similar patterns of activity in the NCI screen. This observation was developed into an algorithm called COMPARE and has been successfully used to predict mechanisms for a wide variety of compounds. More recently, this method has been extended to associate patterns of cell growth inhibition by compounds with measurements of molecular entities (such as gene expression) in the cell lines in the NCI screen. The COMPARE method and associated data are freely available on the Developmental Therapeutics Program (DTP) web site (http://dtp.nci.nih.gov/). Examples of the use of COMPARE on these web pages will be explained and demonstrated. Published by Elsevier Science Inc.

Identification of small-molecule inhibitors of the colorectal cancer oncogene Krüppel-like factor 5 expression by ultrahigh-throughput screening.

The transcription factor Krüppel-like factor 5 (KLF5) is primarily expressed in the proliferative zone of the mammalian intestinal epithelium, where it regulates cell proliferation. Studies showed that inhibition of KLF5 expression reduces proliferation rates in human colorectal cancer cells and intestinal tumor formation in mice. To identify chemical probes that decrease levels of KLF5, we used cell-based ultrahigh-throughput screening (uHTS) to test compounds in the public domain of NIH, the Molecular Libraries Probe Production Centers Network library. The primary screen involved luciferase assays in the DLD-1/pGL4.18hKLF5p cell line, which stably expressed a luciferase reporter driven by the human KLF5 promoter. A cytotoxicity counterscreen was done in the rat intestinal epithelial cell line, IEC-6. We identified 97 KLF5-selective compounds with EC(50) < 10 µmol/L for KLF5 inhibition and EC(50) > 10 µmol/L for IEC-6 cytotoxicity. The two most potent compounds, CIDs (PubChem Compound IDs) 439501 and 5951923, were further characterized on the basis of computational, Western blot, and cell viability analyses. Both of these compounds, and two newly synthesized structural analogs of CID 5951923, significantly reduced endogenous KLF5 protein levels and decreased viability of several colorectal cancer cell lines without any apparent impact on IEC-6 cells. Finally, when tested in the NCI-60 panel of human cancer cell lines, compound CID 5951923 was selectively active against colon cancer cells. Our results show the feasibility of uHTS in identifying novel compounds that inhibit colorectal cancer cell proliferation by targeting KLF5.

Homology model of RSK2 N-terminal kinase domain, structure-based identification of novel RSK2 inhibitors, and preliminary common pharmacophore.

Ribosomal S6 kinase 2 (RSK2) is a serine/threonine kinase that plays a role in human cancer and Coffin-Lowry syndrome and is comprised of two nonidentical kinase domains, each domain with its own ATP-binding site. RSK2 can be inactivated by different types of small organic molecules. Potent RSK2 inhibitors include the two classic bisindole maleimide PKC inhibitors, Ro31-8220 and GF109203X, and the natural product SL0101 that was shown to bind specifically to the ATP pocket of the N-terminal domain (NTD). In this paper, we present an atomic model of the RSK2 NTD (residues 68-323), which was built to simultaneously bind the distinctive molecular scaffolds of SL0101, Ro31-8220, and GF109203X. The RSK2 NTD model was used to identify two novel RSK2 inhibitors from the National Cancer Institute open chemical repository and to develop a preliminary structure-based pharmacophore model.

Structure-aided optimization of kinase inhibitors derived from alsterpaullone.

In order to perform computer-aided design of novel alsterpaullone derivatives, the vicinity of the entrance to the ATP-binding site was scanned for areas that could be useful as anchoring points for additional protein-ligand interactions. Based on the alignment of alsterpaullone in a CDK1/cyclin B homology model, substituents were attached to the 2-position of the parent scaffold to enable contacts within the identified areas. Synthesis of the designed structures revealed three derivatives (3-5) with kinase-inhibitory activity similar to alsterpaullone. The novel 2-cyanoethylalsterpaullone (7) proved to be the most potent paullone described so far, exhibiting inhibitory concentrations for CDK1/ cyclin B and GSK-3beta in the picomolar range.

Conformational sampling of the botulinum neurotoxin serotype A light chain: implications for inhibitor binding.

Botulinum neurotoxins (BoNTs) are the most potent of the known biological toxins, and consequently are listed as category A biowarfare agents. Currently, the only treatments against BoNTs include preventative antitoxins and long-term supportive care. Consequently, there is an urgent need for therapeutics to counter these enzymes--post exposure. In a previous study, we identified a number of small, nonpeptidic lead inhibitors of BoNT serotype A light chain (BoNT/A LC) metalloprotease activity, and we identified a common pharmacophore for these molecules. In this study, we have focused on how the dynamic movement of amino acid residues in and surrounding the substrate binding cleft of the BoNT/A LC might affect inhibitor binding modes. The X-ray crystal structures of two BoNT/A LCs (PDB refcodes=3BTA and 1E1H) were examined. Results from these analyses indicate that the core structural features of the examined BoNT/A LCs, including alpha-helices and beta-sheets, remained relatively unchanged during 1 ns dynamics trajectories. However, conformational flexibility was observed in surface loops bordering the substrate binding clefts in both examined structures. Our analyses indicate that these loops may possess the ability to decrease the solvent accessibility of the substrate binding cleft, while at the same time creating new residue contacts for the inhibitors. Loop movements and conformational/positional analyses of residues within the substrate binding cleft are discussed with respect to BoNT/A LC inhibitor binding and our common pharmacophore for inhibition. The results from these studies may aid in the future identification/development of more potent small molecule inhibitors that take advantage of new binding contacts in the BoNT/A LC.

Synthesis of paullones with aminoalkyl side chains.

Paullones 3 and 4 with aminoalkyl side chains in 2- or 3-position were synthesized as derivatives of kenpaullone 1. Both 3 and 4 showed the characteristic CDK1-inhibitory activity of the paullones and a modest antiproliferative activity on cultured human tumor cell lines. Hence, 3 and 4 appear to be suitable tools for affinity studies directed to find additional intracellular paullone targets.

CDK1-inhibitory activity of paullones depends on electronic properties of 9-substituents.

Multiple linear regression analysis was employed in an effort to establish a quantitative structure-activity relationship model for the CDK1-inhibitory activity of a series of 9-substituted paullones. While the electronic properties of the 9-substituents proved to be of high relevance for CDK1 inhibition, both lipophilic and a steric parameters could not be included in a meaningful equation for the calculation of biological properties. The equation solely based on the electronic parameter was successfully used for the prediction of the CDK1-inhibitory activity of a small test set comprising novel paullones with sulfur-containing 9-substituents. Among these new derivatives, 2-methoxy-9-methylsulfonylpaullone proved to be superior to the standard alsterpaullone with respect to CDK1 inhibition.

Poly(ethylene glycol) prodrugs of the CDK inhibitor, alsterpaullone (NSC 705701): synthesis and pharmacokinetic studies.

Two methods were devised to conjugate PEG to alsterpaullone (NSC 705701) via the N of the indole ring portion of the molecule. In the first approach, activation of the indole was accomplished by reaction with p-nitrophenyl chloroformate to produce a reactive carbamate that was then condensed with a mono blocked diamine to form a urea bond followed by deblocking and conjugation to PEG. The second route used the anion of the indole and produced a carbamate bond. Both compounds were highly water soluble, were stable in buffer, and released alsterpaullone in vitro and in vivo. Studies were conducted in mice to investigate the influence of PEGylation on the plasma pharmacokinetics of alsterpaullone. The total plasma clearance rate was decreased up to 32-fold, and the biological halflife lengthened up to 8-fold when alsterpaullone was injected i.v. as a PEG-conjugate and compared to injection of the unconjugated compound. The most pronounced effect on the pharmacokinetics of alsterpaullone was produced by a 40-kDa PEG urea-linked conjugate. When the 40- and 20-kDa urea-linked conjugates were administered by i.p. injection, high relative bioavailability (46% and 99%, respectively) of alsterpaullone was observed.

Identification of a potent and selective pharmacophore for Cdc25 dual specificity phosphatase inhibitors.

Small molecules provide powerful tools to interrogate biological pathways but many important pathway participants remain refractory to inhibitors. For example, Cdc25 dual-specificity phosphatases regulate mammalian cell cycle progression and are implicated in oncogenesis, but potent and selective inhibitors are lacking for this enzyme class. Thus, we evaluated 10,070 compounds in a publicly available chemical repository of the National Cancer Institute for in vitro inhibitory activity against oncogenic, full-length, recombinant human Cdc25B. Twenty-one compounds had mean inhibitory concentrations of <1 microM; >75% were quinones and >40% were of the para-naphthoquinone structural type. Most notable was NSC 95397 (2,3-bis-[2-hydroxyethylsulfanyl]-[1,4]naphthoquinone), which displayed mixed inhibition kinetics with in vitro K(i) values for Cdc25A, -B, and -C of 32, 96, and 40 nM, respectively. NSC 95397 was more potent than any inhibitor of dual specificity phosphatases described previously and 125- to 180-fold more selective for Cdc25A than VH1-related dual-specificity phosphatase or protein tyrosine phosphatase 1b, respectively. Modification of the bis-thioethanol moiety markedly decreased enzyme inhibitory activity, indicating its importance for bioactivity. NSC 95397 showed significant growth inhibition against human and murine carcinoma cells and blocked G(2)/M phase transition. A potential Cdc25 site of interaction was postulated based on molecular modeling with these quinones. We propose that inhibitors based on this chemical structure could serve as useful tools to probe the biological function of Cdc25.

Novel small molecule inhibitors of botulinum neurotoxin A metalloprotease activity.

Botulinum neurotoxins (BoNTs) are among the most lethal biological substances to have been weaponized and are listed as biodefense category A agents. Currently, no small molecule (non-peptidic) therapeutics exist to counter this threat; hence, identifying and developing compounds that inhibit BoNTs is a high priority. In the present study, a high-throughput assay was used to identify small molecules that inhibit the metalloprotease activity of BoNT serotype A light chain (BoNT/A LC). All inhibitors were further verified using a HPLC-based assay. Conformational analyses of these compounds, in conjunction with molecular docking studies, were used to predict structural features that contribute to inhibitor binding and potency. Based on these results, a common pharmacophore for BoNT/A LC inhibitors is proposed. This is the first study to report small molecules (non-peptidics) that inhibit BoNT/A LC metalloprotease activity in the low microM range.