SUMMARY
The actin-bundling protein fascin is a key mediator of tumor invasion and metastasis and its activity drives filopodia formation, cell-shape changes and cell migration. Small-molecule inhibitors of fascin block tumor metastasis in animal models. Conversely, fascin deficiency might underlie the pathogenesis of some developmental brain disorders. To identify fascin-pathway modulators we devised a cell-based assay for fascin function and used it in a bidirectional drug screen. The screen utilized cultured fascin-deficient mutant Drosophila neurons, whose neurite arbors manifest the ‘filagree’ phenotype. Taking a repurposing approach, we screened a library of 1040 known compounds, many of them FDA-approved drugs, for filagree modifiers. Based on scaffold distribution, molecular-fingerprint similarities, and chemical-space distribution, this library has high structural diversity, supporting its utility as a screening tool. We identified 34 fascin-pathway blockers (with potential anti-metastasis activity) and 48 fascin-pathway enhancers (with potential cognitive-enhancer activity). The structural diversity of the active compounds suggests multiple molecular targets. Comparisons of active and inactive compounds provided preliminary structure-activity relationship information. The screen also revealed diverse neurotoxic effects of other drugs, notably the ‘beads-on-a-string’ defect, which is induced solely by statins. Statin-induced neurotoxicity is enhanced by fascin deficiency. In summary, we provide evidence that primary neuron culture using a genetic model organism can be valuable for early-stage drug discovery and developmental neurotoxicity testing. Furthermore, we propose that, given an appropriate assay for target-pathway function, bidirectional screening for brain-development disorders and invasive cancers represents an efficient, multipurpose strategy for drug discovery.
INTRODUCTION
A highly conserved actin-bundling protein, fascin has diverse roles in the developmental and physiological regulation of cellular morphology and function (Kureishy et al., 2002; Jayo and Parsons, 2010; Sedeh et al., 2010; Hashimoto et al., 2011). It is also implicated in human disease pathogenesis, under both loss-of-function and gain-of-function conditions, which motivated us to develop a fascin bioassay for drug discovery. Note that fascin is unrelated to either the fasciclins or neurofascin, which are members of the immunoglobulin cell-adhesion molecule superfamily. Fascin drives the formation of cell-membrane protrusions, including lamellipodia (Yamashiro et al., 1998), microspikes (Svitkina et al., 2003), filopodia (Vignjevic et al., 2006) and invadopodia (Li et al., 2010a), in part because F-actin bundles increase mechanical stiffness (Tseng et al., 2005; Vignjevic et al., 2006). In addition, fascin-mediated actin bundling and crosslinking, which are regulated by phosphorylation (Ono et al., 1997; Aratyn et al., 2007) and the extracellular matrix (ECM), enhance cell migration (Ono et al., 1997; Yamashiro et al., 1998; Anilkumar et al., 2003; Jawhari et al., 2003) and ECM degradation (Li et al., 2010a).
Mammals have three fascin-coding genes, of which Fascin-2 and Fascin-3 are expressed in narrow domains (Tubb et al., 2000; Tubb et al., 2002; Shin et al., 2010), whereas Fascin-1 is broadly and dynamically expressed. Fascin-1 is abundant early in development, especially in the central nervous system (CNS) and migrating cells, and is then downregulated as cells mature (De Arcangelis et al., 2004; Zhang et al., 2008; Zanet et al., 2009; Tang et al., 2010). In this paper, ‘fascin’ refers to the product of the Fascin-1 genes (FSCN1 in humans, MIM#602689; Fscn1 in mouse; and singed in Drosophila, FBgn0003447).
Fascin has a pivotal role in tumor invasion and metastasis (Machesky and Li, 2010), leading to the proposal that fascin-blocking drugs might prevent the spread of malignant cancers (Yoder et al., 2005; Hashimoto et al., 2011). Because most cancer-related deaths are due to metastases, there is an urgent need for development of anti-metastasis agents (Sporn, 1996; Sleeman and Steeg, 2010). For carcinomas from numerous organs, high fascin expression levels are associated with increased invasiveness and earlier patient death (Machesky and Li, 2010). Fascin is also involved in tissue infiltration by circulating tumor cells (Kim et al., 2009a). Similar associations have been reported for malignant glioblastomas (Peraud et al., 2003; Roma and Prayson, 2005; Gunal et al., 2008), which are prone to extensive dispersion within the CNS (Giese et al., 2003).
TRANSLATIONAL IMPACT
Clinical issue
The findings reported here address two major clinical problems that do not initially seem related. First, invasive cancers, including most brain tumors, cause death due to invasion or migration, which are not inhibited by currently available cancer treatments. Second, developmental brain disorders are not treatable with drugs that enhance cognitive function. For both of these unmet medical needs, a major obstacle has been the lack of cellular bioassays for compound screening. The actin-bundling protein fascin links these two challenging clinical conditions: excess fascin promotes tumor invasion and metastasis, whereas insufficient fascin disrupts brain development. Thus, the fascin pathway represents a highly desirable drug target.
Results
The filagree phenotype of fascin-deficient mutant Drosophila neurons enabled the authors to develop a bidirectional in vitro cellular bioassay to screen for drugs that modify the fascin pathway. A library of 1040 known compounds (NINDS-II) was chosen on the basis of high molecular diversity, and was screened with the aim of identifying drugs that could be repurposed for new indications. Of these compounds, 81 were active as fascin-pathway modifiers. There was wide pharmacological and chemical-structure diversity in each set of active compounds (34 blockers and 48 enhancers), strongly suggesting that each set has multiple targets along the fascin pathway. Comparison of closely related compounds that differ in activity provided structure-activity relationship (SAR) hypotheses that can be tested in follow-up studies. Notably, all four of the statin compounds in the library caused a unique, reversible neurotoxic morphological effect characterized by intraneurite nodules containing aggregations of organelles [‘beads-on-a-string’ (BOS)]. Fascin deficiency enhances the sensitivity of neurons to BOS.
Implications and future directions
These findings introduce a conceptually simple cell-based fascin bioassay and apply it to identify many compounds and preliminary SAR information that can be pursued for drug development, either by repurposing or lead optimization. Fascin-pathway blockers could serve as anti-invasion and anti-metastasis agents for patients with malignant carcinomas or gliomas. Fascin-pathway enhancers could improve neurocognitive function and behavior in a subset of children with developmental brain disorders. The potential of fascin-pathway modifiers warrants testing in various mammalian models of fascin-sensitive human disease. In addition, bidirectional cell-based screening could be applied to other important biological pathways with dual functions in brain development and tumor invasion.
The authors propose that the statin-induced BOS they observed represents a cellular correlate of statin-associated cognitive side effects experienced by some patients on statins. Because cholesterol is an essential nutrient in Drosophila, this system provides a unique opportunity to study statin-mediated neurotoxicity separately from cholesterol biosynthesis. The authors’ demonstration that genetic background impacts the sensitivity of neurons to statin-induced BOS suggests experiments to identify potential human genetic risk factors for statin-associated cognitive deficits. More generally, this primary neuron culture system in Drosophila holds great promise as a neurotoxicity screening platform with the ability to identify biologically relevant gene-by-environment interactions.
The causal role of fascin in tumor phenotypes is supported by laboratory studies in which blocking of fascin expression reduced the invasive and/or metastatic properties of colon carcinoma (Hashimoto et al., 2007), glioblastoma multiforme (Hwang et al., 2008), gastric carcinoma (Fu et al., 2009; Kim et al., 2010) and breast carcinoma (Chen et al., 2010; Al-Alwan et al., 2011). Moreover, migrastatin-family compounds, which are potent inhibitors of tumor invasion and metastasis in the laboratory (Shan et al., 2005), bind to fascin and inhibit its actin-bundling activity (Chen et al., 2010). Thus, the experimental and clinical data strongly suggest that a fascin bioassay would be of great value for discovery of drugs with anti-invasion and/or anti-metastasis activity.
Fascin is also required for normal brain development (Yamakita et al., 2009) (R.K. and L.L.R., unpublished data), plausibly by regulating neuronal differentiation (Deinhardt et al., 2011; Marín-Vicente et al., 2011). Fascin insufficiency or dysregulation might underlie disorders of brain development and plasticity, resulting in intellectual disability (Kraft et al., 2006). Fascin regulation is likely to be faulty in the brain-development disorders Rubinstein-Taybi syndrome (Roelfsema and Peters, 2007) and tuberous sclerosis (Ess, 2006), which are caused by mutations in CREBBP (MIM#600140) and TSC1 (MIM#605284) or TSC2 (MIM#191092), respectively. FSCN1 is an upregulated transcriptional target of CREB binding protein (CREBBP) (Megiorni et al., 2005), whereas fascin protein is a target of the TSC1-TSC2 complex (Gan et al., 2008).
The connection between fascin and brain plasticity has also been revealed by unbiased proteomics screens. Reduced fascin levels were found in two mouse models of absence epilepsy (Ryu et al., 2007; Ryu et al., 2008) and after long-term memory induction (Li et al., 2010b). These data suggest that fascin levels are downregulated by neural activity, perhaps to permit synapse-structure changes. By contrast, fascin levels were elevated in a polytransgenic mouse model of Down syndrome (Shin et al., 2007). This could be a molecular feature of the brain-development disorder or early-onset neurodegeneration, or both. In a neuron-culture model of neuroprotection, rapid induction of ischemic tolerance was associated with ubiquitylation of fascin, which was subsequently degraded, as well as with transient retraction of dendritic spines. This is consistent with the idea that fascin removal allows dissolution of actin bundles, thereby accelerating synapse remodeling (Meller, 2009). Finally, in a canine model of aging, slowing of cognitive decline by environment and diet was associated with reduced levels of fascin carbonyl, a marker of oxidation (Opii et al., 2008). This suggests that preventing fascin oxidation contributes to better cognitive performance. These reports point to a role for fascin in regulating neuronal differentiation and synaptic plasticity, which are disrupted in brain-development disorders (Johnston, 2004) and vulnerable during aging (Burke and Barnes, 2006). Hence, pharmacological enhancers of fascin expression or function could be beneficial for diverse neurodevelopmental and cognitive or behavioral conditions.
Loss-of-function Fascin-1 mutations are available in Drosophila melanogaster and Mus musculus. The targeted Fscn1 disruption in the mouse causes structural brain abnormalities and high rates of neonatal death (Yamakita et al., 2009). Drosophila has a single fascin-coding gene (Bryan et al., 1993; Kureishy et al., 2002), named singed for the gnarled bristles of mutant flies (Bender, 1960; Cant et al., 1994; Tilney et al., 1995; Wulfkühle et al., 1998). Wild-type singed function is also essential for oogenesis (Cant et al., 1994), blood cell migration (Zanet et al., 2009) and some aspects of brain development (R.K. and L.L.R., unpublished data), but not for viability (Bender, 1960). We discovered that singed mutations cause a neuronal morphogenesis defect, ‘filagree’, when developing mutant CNS neurons are cultured in vitro (Kraft et al., 2006). Filagree neurite arbors have exaggerated clockwise curvature and erratic variation in neurite caliber, but no reduction in length, branching or axon-dendrite ratio (Kraft et al., 2006). Filagree neurons have striking disruptions of their actin cytoskeleton which, like other singed phenotypes, can be explained by failure of actin bundling (Cant et al., 1994; Guild et al., 2003; Kraft et al., 2006; Zanet et al., 2009). Because decreasing fascin function causes increasing severity of the filagree phenotype, we reasoned that cultured neurons could provide a cellular bioassay for fascin function.
The use of a Drosophila cell-based fascin bioassay for drug discovery is justified by the phylogenetic conservation of fascin (Hashimoto et al., 2011) and of the pathways underlying cognition (Greenspan and Dierick, 2004; Inlow and Restifo, 2004; Bolduc and Tully, 2009), tumor invasion and metastasis (Miles et al., 2011). Here we present the results of a bidirectional screen for pharmacological modifiers of the filagree neuronal morphogenesis defect.
RESULTS
Drug screen design using the fascin bioassay
Filagree is a highly penetrant, quantifiable, pan-neuronal phenotype of fascin-deficient singed-mutant neurons cultured from the developing CNS of Drosophila mature larvae or young pupae (Fig. 1A–C) (Kraft et al., 2006). This developmental interval corresponds to high levels of singed transcript accumulation in wild-type CNS (Kaitlin L. Bergfield, R.K. and L.L.R., unpublished). The filagree phenotype is so distinctive and so consistent that trained observers can easily distinguish singed-mutant versus wild-type population cultures, typically containing 1500–2000 neurons each. In fact, both humans and a computational method for neurite-curvature quantification can classify individual photomicrographs of mutant or wild-type neurons with =90% accuracy (Kraft et al., 2006). To anticipate how much restoration of function would be required to rescue the fascin-deficient filagree phenotype, we asked whether filagree is recessive, like other singed phenotypes. We quantified the neurite-curvature distributions of wild-type and snX2/+ cultured larval CNS neurons, which have 100% and 50% of normal fascin function, respectively (Fig. 1D). The two distributions were statistically indistinguishable, whereas neurons from sn3/snX2 (a near-null genotype) and snX2/+ were very different (Welch’s t-test, P<1×10−18; Fig. 1D). Therefore, the results show that the singed neurite-curvature phenotype is recessive and that rescue of singed-mutant neurite-arbor morphology requires no more than 50% of wild-type fascin function.
Genetic and pharmacological modification of neurite curvature. (A–C,F–I) Phase-contrast images (60×) of neurons cultured for 3 d.i.v. from the CNS of wandering third instar larvae. Magnification is the same throughout. (A–C) Loss-of-function singed mutations increase neurite curvature and disrupt proximal-to-distal tapering. (A) Wild-type (sn+/sn+), OreR-C laboratory strain. (B) A partial loss-of-function (hypomorphic) mutation (sn3/sn3) causes a moderate filagree phenotype. (C) A null mutation (snX2/Y) causes a severe filagree phenotype. (D,E) Mean curvature distributions of neurite arbors of cultured neurons with differing genotypes, plotted with soft binning. Normal-distribution curves were fit to each population and scaled (y-axis) to the corresponding histograms. (D) Genetically marked γ-MB neurons. The increased curvature of sn3/snX2 (blue) neurons is easily seen. The similarity of snX2/+(red) and sn+/sn+ (green) curvature distributions demonstrates that filagree is recessive. (E) Random CNS neurons. The mean neurite curvature distribution of sn3/sn3 mutant neurons (magenta) is significantly increased over that of wild-type (sn+/sn+, OreR-C) neurons (green). (F–I) Exposure to drugs in vitro can modify the filagree phenotype of fascin-deficient sn3/sn3 neurons. (F,G) Two examples of filagree decreasers (fascin-pathway enhancers): estradiol propionate, 50 μM (E) and Anisindione, 50 μM (G). Neurite curvature, especially of terminal neurites, is reduced and the smooth proximal-to-distal tapering is restored. (H,I) Two examples of filagree increasers (fascin-pathway blockers): griseofulvin, 50 μM (H) and acetyltryptophan, 50 μM (I). Note the exaggerated neurite curvature and frequent expansions of neurite width.
A bidirectional drug-screen design requires that the baseline neuronal morphology be intermediate and distinct from both wild type (Fig. 1A) and singed-null (snX2/Y; Fig. 1C), allowing detection of both drug-induced worsening and drug-induced rescue of the filagree defect (Fig. 1B). The neurite-curvature distributions of the severe hypomorph,sn3/sn3 (Fig. 1B), and the wild-type (Fig. 1A) were significantly different (Welch’s t-test, P<1×10−13; Fig. 1E), with greater overlap of distributions than between wild type and sn3/snX2 (Kraft et al., 2006). This is consistent with the relative function among the genotypes. Moreover, the sn3 mutation leaves the fascin open reading frame intact, allowing small amounts of normal protein to be produced (Paterson and O’Hare, 1991; Cant et al., 1994). Thus, sn3/sn3 neurons have a desirable neuronal phenotype, and contain some wild-type fascin protein that can serve as a drug target.
The striking nature of the filagree defect makes it suitable for use in a modifier screen based on direct phenotype observation. Using an approach that parallels genetic-modifier screens commonly performed with model organisms, we conducted a chemical screen to identify small molecules that induce obvious changes in neuronal phenotype. Compounds were tested at 10 μM and 50 μM by adding them to the cultures at the time of plating the dissociated neurons, and evaluating their effects on the filagree phenotype after 3 days in vitro (d.i.v.). Thus, the neurons were exposed to the drug during the entire period of neurite-arbor morphogenesis. To facilitate the drug-screening process, we replaced quantitative analysis of immunostained images of randomly sampled neurons with holistic scoring by phase-contrast microscopy of populations of living cells, typically 1500–2000 neurons per ∼50-mm2 culture well.
Selection and diversity analysis of the NINDS-II compound collection
Inspired by striking examples of new uses for existing drugs, we took a ‘repurposing’ approach, also called ‘repositioning’ or ‘indication switch’ (Ashburn and Thor, 2004; Dueñas-González et al., 2008). For example, minoxidil was developed as an oral anti-hypertensive drug, but is now most popular as a topical treatment (e.g. Rogaine®) for hair loss (Zins, 1988). More dramatic is the repurposing of thalidomide, much-maligned for its teratogenic effects when taken during pregnancy (Ito and Handa, 2012), but now used to treat multiple myeloma, a bone marrow malignancy (Palumbo et al., 2008). Repurposing screens have the advantage of testing compounds whose pharmacology and side-effect profiles are partially understood.
Because fascin acts in diverse pathways (Jayo and Parsons, 2010; Hashimoto et al., 2011), modification of the filagree defect could occur through several potential drug targets, including but not limited to fascin itself. This called for a chemically diverse screening library. The new indications, including brain tumors and neurodevelopmental disorders, made good representation of neuroactive compounds highly desirable. A collaboration between the National Institute of Neurological Disorders and Stroke (NINDS) and the private sector (Heemskerk, 2005) produced a screening library of 1040 compounds (NINDS-II; supplementary material Table S1), comprised primarily of FDA-approved drugs for diverse indications, as well as natural products and laboratory reagents; many of these compounds are neuroactive. Screens of the NINDS-II library with various assays and protocols identified drugs with potential utility for Huntington’s disease, spinomuscular atrophy, amyotrophic lateral sclerosis (ALS), stroke and familial dysautonomia (Aiken et al., 2004; Piccioni et al., 2004; Slaugenhaupt et al., 2004; Stavrovskaya et al., 2004; Rothstein et al., 2005; Vincent et al., 2005; Wang et al., 2005a; Wang et al., 2005b; Desai et al., 2006). The antibiotic ceftriaxone was protective in diverse neurodegeneration assays, leading to Phase III clinical trials for ALS (Traynor et al., 2006) (ClinicalTrials.gov identifier NCT00349622). This library does not contain any migrastatin-family compounds.
We evaluated the chemical diversity of the NINDS-II compounds (Fig. 2) based on: (i) the types and distribution of molecular scaffolds (Bemis and Murcko, 1996; Singh et al., 2009), (ii) measures of molecular similarity based on shared substructural features (Willett et al., 1998) and (iii) the chemical-space distribution of the compounds (Maggiora and Shanmugasundaram, 2011). The 1040 compounds of the collection are distributed over 617 scaffolds (Fig. 2B), of which 551 are cyclic systems. The highly populated scaffolds (Fig. 2C) are indoles, pyridines, quinolines and sterols, common among drug-like small molecules. Indicative of high diversity, 77% (424) of the cyclic scaffolds are singletons, i.e. populated by a single compound. From the perspective of the compounds, ∼60% of them cover ∼23% of the cyclic scaffolds. In other words, this relatively small collection allows sampling of a relatively large number of scaffolds. The molecular fingerprint, based on 166 structural features (Fig. 2D), of each compound was compared pair-wise with all others, and the Tanimoto coefficients were computed as measures of similarity (Maggiora and Shanmugasundaram, 2011). Overall molecular diversity is inversely proportional to the average Tanimoto similarity for all pair-wise comparisons. For the NINDS-II library, low similarity values (mean 0.295; median 0.284) reflect high diversity and compare favorably to those of the ∼1500-compound DrugBank database (Wishart et al., 2008; Singh et al., 2009). The distribution of compounds in ‘chemical space,’ where the similarity of any pair is inversely related to the distance between them (Maggiora and Shanmugasundaram, 2011), was represented in a three-dimensional (3D) plot (Fig. 2E,F). The compounds are spread throughout the space, reflecting the high diversity of the collection. In summary, the high diversity of the modestly sized NINDS-II collection make it an efficient tool for conducting the first compound screen on the basis of the first fascin bioassay.
Molecular diversity of the NINDS-II compound collection screened in this study. (A) Relationship between a molecule and its scaffold, using tranilast as an example. The scaffold or cyclic system is obtained after iteratively removing the side chains and is identified by a five-character chemotype identifier code. The WSLAM scaffold is a singleton, i.e. tranilast is the only compound in the collection with this scaffold. (B) Scaffold-recovery curve. The x-axis represents the scaffolds, organized left-to-right from the most common (benzene) to the singletons, plotted as a cumulative fraction of the total (552, with all 66 acyclic systems considered a single scaffold). The y-axis shows the cumulative fraction of the 1040 compounds occupied by a given fraction of the scaffolds: ∼14% of the scaffolds contain 50% of the NINDS-II library (red lines), and 60% of the library is distributed over ∼24% of the scaffolds (green lines). Beyond this point, each scaffold is a singleton and the relationship becomes linear. (C) Common scaffolds in the NINDS-II library. 2D representations of cyclic systems found with a frequency ≤5; the number is under each structure. Highlighted are the pyridines (red), the quinolines (blue), and the scaffolds containing a four-ring sterol (green; 36 compounds). (D) Representation of chemical fingerprint analysis based on MACCS keys, using the antipsychotic drug droperidol as an example. Each of 166 possible substructural features is coded using a bit score: present at least once (1) or absent (0). Four of droperidol’s features (CN2O, fluorine, five-membered heterocycle and C3N) are shown with their positions in the MACCS keys bitstring. (E,F) 3D representations of the NINDS-II library in chemical space, obtained by principal component analysis (PCA) of the similarity matrix computed using MACCS keys and Tanimoto similarity. Each compound is shown as a sphere. The chemical space distribution of the compounds is extensive, with regional variation in population density. The 17 antipsychotic drugs (blue), are localized to a medium-sized sector of the space and interspersed among compounds with other activities. This is consistent with these drugs occupying 14 scaffolds and representing several different pharmacological mechanisms of action. (E) View centered on the vertical (PCA 2) axis. (F) The view obtained by rotation around the vertical axis. One antipsychotic drug is hidden behind another compound.
A bidirectional drug screen reveals diverse fascin-pathway enhancers and blockers
Of the 1040 compounds tested, 81 (7.8%) were active in the fascin bioassay at 10 and/or 50 μM concentrations, based on holistic scoring. We identified both drug-induced decreasers and increasers of the filagree phenotype (Figs 1, 3, 4; Table 1). Filagree decreasers (fascin-pathway enhancers; supplementary material Table S2) rescued the neurite-arbor shape defect of fascin-deficient mutant neurons, allowing them to extend neurites with normal trajectory and tapering (Fig. 1, compare F and G with A). Filagree increasers (fascin-pathway blockers; supplementary material Table S3) worsened the clockwise curvature and further disrupted neurite tapering, as seen in the fascin-null-mutant neurons (Fig. 1, compare H and I with C).
Compounds active in the fascin bioassay
Summary of the fascin bioassay-based drug screen. (A) Compounds were screened by holistic scoring for their ability to modify the filagree phenotype, revealing two classes of activity, filagree decreasers (fascin-pathway enhancers) and filagree increasers (fascin-pathway blockers). One compound, cloxyquin, had dose-dependent activity (asterisk) and is counted in both groups. Confirmed active compounds were classified as ‘high potency’ (active at 10 and 50 μM) or ‘low potency’ (active at only 50 μM). Further classification distinguished between simple actives and those with cytotoxicity at the high dose or morphological effects (Plus). (B,C) Phase-contrast photomicrographs after 3 d.i.v. Cloxyquin increased the filagree phenotype at 10 μM (B), but decreased the filagree phenotype at 50 μM (C). (D) Therapeutic classes of filagree-modifying compounds, based on clinical and laboratory data. A compound may be in more than one therapeutic class. ‘Other’ includes compounds used for drug overdose, local anesthesia, gout, eye surgery, allergy, as homeopathic remedy or as insecticide. See supplementary material Tables S2 and S3 for information on chemical class, pharmacology and therapeutic use.
Chemical space distribution of compounds active in the fascin bioassay. 3D views of the chemical space occupied by NINDS-II, centered on the vertical (PCA 2) axis and highlighting different subsets of fascin-pathway modifiers. (A,B) Each class of active compounds is distributed throughout a large fraction of the chemical space occupied by the whole collection. Each compound is represented as a sphere with (A) the 48 fascin-pathway enhancers in green and (B) the 34 fascin-pathway blockers in magenta against the gray background of the inactive compounds. (C,D) Fascin-pathway modifiers of the same potency class, plotted in contrasting colors without the whole collection. Within each class, there is some intermingling of compounds with opposite activities, and some clustering of compounds with the same activity. (C) High-potency fascin-pathway enhancers (green; 24) and blockers (magenta; 12). (D) Low-potency fascin-pathway enhancers (green; 24) and blockers (magenta; 22). See supplementary material Fig. S2 for chemical structures of all active compounds.
Compounds with activity at 10 μM were considered ‘high potency’ whereas those with activity only at 50 μM were considered ‘low potency’. The 48 fascin-pathway enhancers were split equally between low- and high-potency compounds; the 34 fascin-pathway blockers included 22 low- and 12 high-potency compounds (Fig. 3A). Cloxyquin was the only compound with a bimodal effect, acting as a filagree increaser at 10 μM, but as a filagree decreaser at 50 μM (Fig. 3B,C). Within each potency class, we made a distinction between ‘simple’ actives (21 fascin-pathway enhancers, 22 fascin-pathway blockers) and those with other effects on the neurons beyond filagree modification (Fig. 3A,D; see below). The compounds of each activity group span a broad range of therapeutic indications (Fig. 3D; supplementary material Tables S2, S3).
When active compounds were re-tested for quantification of neurite curvature, we found that holistic scoring is relatively insensitive and very likely to underestimate true drug potencies. For example, the antifungal agent griseofulvin, initially identified as a low-potency fascin-pathway blocker (Fig. 1H), was active at both 10 μM (P=0.011) and 50 μM (P=0.0002), causing a dose-dependent increase in neurite curvature. Surprisingly, the mushroom body γ-neurons (γ-MB neurons) in the griseofulvin-treated cultures failed to show a significant change in neurite curvature at either 10 μM (P=0.13) or 50 μM (P=0.08). This also raises the question of whether any of the ten filagree-modifying compounds whose activity failed to replicate in confirmation tests (Fig. 3A) might have been too close to the detection threshold by holistic scoring, but nonetheless have some activity. This highlights the differential drug sensitivity of classes of neurons in the fascin-deficient brain. It is therefore possible that we did not detect drugs whose effects were restricted to neuron subtypes.
The 3D chemical-space distribution of the 48 fascin-pathway enhancers (filagree decreasers) and 34 fascin-pathway blockers (filagree increasers) are shown plotted against the background of the compound collection (Fig. 4, A and B, respectively). Each set of active compounds is quite diverse, as indicated by its distribution throughout the chemical space of the collection without obvious localization to any particular neighborhood. A notable exception is the pair of high-potency fascin-pathway blockers, colchicine and colchiceine (Fig. 4C), which are closely related alkaloids that differ only by a methyl group (Tanimoto similarity 0.956). When the fascin-pathway enhancers and blockers were plotted together in chemical space, many instances of intermingling between the two were seen among both high- and low-potency active compounds (Fig. 4, C and D, respectively). The chemical structures of the high-potency simple actives are shown in Fig. 5, along with 2D chemical-space plots. The ten simple high-potency fascin-pathway enhancers (Fig. 5C) and the five simple high-potency fascin-pathway blockers (Fig. 5D) each belong to a different scaffold.
Structures of high-potency simple fascin-pathway enhancers and blockers. (A,B) 2D chemical-space distributions of the high-potency simple filagree modifiers plotted against the background of the NINDS-II compound collection in gray. Each compound is represented as a circle; due to the projection, some compounds are hidden from view. (A) Ten fascin-pathway enhancers (green). (B) Five fascin-pathway blockers (magenta). (C,D) 2D chemical structures of the high-potency simple filagree modifiers. Within each class, there is considerable structural diversity. (C) Fascin-pathway enhancers. (D) Fascin-pathway blockers.
Consistent with this level of structural diversity, the active compounds represent diverse chemical classes and wide-ranging pharmacological mechanisms of action (supplementary material Tables S2, S3). For example, drugs that directly impact neurotransmitter systems were found among both fascin-pathway enhancers (acetylcholine, adiphenine HCl, mephenesin, metaraminol bitartrate, spaglumic acid) and fascin-pathway blockers (amantadine, baclofen, imipramine, spiperone). Several of the steroid hormones in the library were fascin-pathway enhancers (cyproterone, estradiol acetate, estradiol propionate, methylprednisolone). Among the fascin-pathway modifiers were two dozen antimicrobial agents, including antibacterial, antiviral, antifungal and anti-amoebic compounds.
Neither pharmacological nor chemical class reliably predict the direction of activity in the fascin bioassay. Of two sulfanilamide-type antibacterial agents, one (sulfamethoxazole) was a fascin-pathway blocker and the other (sulfamethazine) was a fascin-pathway enhancer. Four non-steroidal anti-inflammatory drugs that inhibit cyclooxygenase (COX) were fascin-pathway blockers (diflunisal, naproxen, salicin, suprofen), whereas other COX inhibitors, including diclofenac and ibuprofen, were inactive. The diversity of the active compounds suggests that multiple drug targets and mechanisms are operative, which is consistent with a screen based on a pathway that is regulated at multiple levels (Jayo and Parsons, 2010). These considerations led to further analysis of the chemical structures of active compounds and their inactive relatives.
Structure-activity relationships of fascin-pathway enhancers and blockers
Due to the complexity of any cell-based assay, the activity of a compound can be affected by its transport, metabolism, binding to the molecular target(s) or a combination thereof. Nonetheless, small subsets of structurally similar compounds in the NINDS-II collection allowed us to conduct a preliminary analysis of structure-activity relationships (SARs) (Fig. 6) based on the concept of activity landscapes. Analogous to geographical landscapes, the latitude and longitude correspond to the position of a molecule in a 2D chemical space, and the altitude corresponds to its activity in a given assay. Because structurally similar compounds tend to have similar biological activities (Johnson and Maggiora, 1990), regions of active molecules in an activity landscape tend to resemble rolling hills. In some cases, however, similar molecules have significantly different activities, leading to cliff-like features. These activity cliffs contain significant information on SARs because they identify relatively small yet biologically relevant structural features associated with the dramatic changes in activity (Maggiora and Shanmugasundaram, 2011). We evaluated cases in which high-potency fascin-pathway enhancers (Fig. 6A–C) or blockers (Fig. 6D–F) were structurally very similar to low-potency or inactive compounds. By comparing the structures at the top of each cliff with those at its base, we could infer which structural features are necessary for activity. In some cases, the SAR hypotheses were also informed by the physical and/or chemical properties of the active and inactive compounds.
SAR analysis based on activity cliffs. (A–F) High-potency fascin-pathway modifiers (gray background) compared with structurally similar low-potency (light gray background) or inactive (white background) compounds. The highlighted substructures can be inferred to influence activity in the fascin bioassay. (A–C) Fascin-pathway enhancers (green highlights). (A) Estradiol propionate: the size of the side chain (oval) on the carboxy group on the sterol D ring is associated with activity, suggesting that in estradiol acetate the side chain is too small, whereas in estradiol valerate and cypionate, it is too large. (B) Sulfamethazine: activity is associated with dimethylation (circles) of the pyrimidine ring. Removal (asterisk) of one (sulfamerazine) or both (sulfadiazine and sulfapyridine) methyl groups is associated with loss of activity. (C) Mephenesin: replacement of the methyl group with a methoxy group in guaifenesin (circles) is associated with lack of activity. (D–F) Fascin-pathway blockers (magenta highlights). (D) Imipramine: in three other tricyclic compounds, single changes are associated with lack of activity. Desipramine is lacking a methyl group of the alkyl nitrogen (circles); clomipramine has a chloride on one of the aromatic rings (asterisks); trimipramine has a methyl substitution in the middle of the alkyl chain (arrows). (E) Suprofen: replacement of a hydrogen with a methyl group (circles) creates suprofen methyl ester, which is inactive. (F) Nifedipine: two changes, the position of the nitro group on the aromatic ring (asterisks) and addition of a methyl group (arrows), are associated with lack of activity of nitrendipine.
Comparing estradiol propionate with estradiol acetate (diminished activity) and two inactive estradiols (similarity to estradiol propionate ≥0.94), the key structural feature appears to be the substituent on the carboxy moiety (Fig. 6A, green ovals). Removing a methyl group lowers the activity, whereas lengthening the alkyl chain (estradiol valerate) or adding a terminal cyclopentane ring (estradiol cypionate) renders the compounds inactive. Because the steroid moiety remains unchanged, it can be inferred that the size of the carbon chain on the carboxy moiety is a crucial structural determinant of fascin-pathway enhancement.
For the high-potency fascin-pathway enhancer sulfamethazine, the focus is the pyrimidine ring, which is dimethylated in the active compound (Fig. 6B, green circles). Three inactive sulfonamides, sulfamerazine, sulfadiazine and sulfapyridine (with similarities to sulfamethazine of 0.98, 0.96 and 0.88, respectively), have reduced numbers of methyl substituents (Fig. 6B, asterisks). In sulfapyridine there is also a carbon atom replacing one of the nitrogens in the pyrimidine ring (Fig. 6B, green arrow). Because removing one or two methyl substituents renders the compound inactive, the importance of removing the ring nitrogen cannot be interpreted. Nevertheless, we can infer that dimethylation of the pyrimidine ring is important for activity. The high-potency fascin-pathway enhancer mephenesin was compared with the structurally related but inactive guaifenesin (similarity 0.87), which differs only by the replacement of the methyl substituent in the phenyl ring by a methoxy group (Fig. 6C, green circles). Thus, it can be inferred that activity is sensitive to substitution at that position of the phenyl ring.
The high-potency fascin-pathway blocker imipramine differs by a single change from each of three inactive tricyclic compounds, with similarities to imipramine of 0.93, 0.91 and 0.88 (Fig. 6D). Desipramine, an imipramine metabolite produced in the mammalian liver, lacks a methyl group on the alkyl nitrogen (Fig. 6D, circles). On the basis of experimentally determined dissociation constants, desipramine is a considerably stronger base than imipramine (Shalaeva et al., 2008). Hence, it has a harder time losing a proton to cross the plasma membrane and enter the neuron. This suggests that access to an intracellular target could be a factor controlling activity of tricyclic compounds in the fascin bioassay. For clomipramine, the single substitution of a chlorine atom on one aromatic ring (Fig. 6D, asterisks) eliminates activity. Lastly, trimipramine differs only by the presence of a methyl substitution on the alkyl chain (Fig. 6D, arrows). From these observations, it can be inferred that three distinct features are essential for imipramine activity.
In Fig. 6E, the high-potency fascin-pathway blocker suprofen is compared with its inactive methyl ester (similarity 0.82). Replacement of a hydrogen atom by a methyl group might disrupt activity because the acid moiety in suprofen is undoubtedly ionized in solution. Thus, the active form of suprofen is most probably an anion, whereas the inactive ester is uncharged. The high-potency fascin-pathway blocker nifedipine has two structural differences from its inactive analog nitrendipine (similarity 0.90). First, nifedipine has a methyl ester whereas nitrendipine has an ethyl ester (Fig. 6F, arrows). Second, the nitro group (Fig. 6F, asterisks) is in the ortho position in nifedipine, but in the meta position in nitrendipine. This change might affect the 3D conformation of nifedipine or indicate that a specific orientation of the nitro group is required for activity.
Other drug effects on cultured neurons
Roughly half of the actives had effects beyond modifying the filagree phenotype (Fig. 3A,D; Table 1; supplementary material Tables S2, S3), which we classified as cytotoxicity or morphological effects. Cytotoxicity was defined as apparent cell death without a neurite arbor (Fig. 7A,B) or, rarely, degeneration of neurites after formation of an arbor (data not shown). A subset of high-potency fascin-pathway modifiers caused overt cytotoxicity at the 50-μM concentration. Some inactive compounds also caused cytotoxicity. Distinctive morphological effects, the most common of which was reduced neurite outgrowth (Fig. 7E,F,H–K), were caused by both active and inactive compounds. Less common, but more dramatic, were changes in neuronal cell-body size or shape (Fig. 7C,D) and bizarre alterations in neurite morphology (Fig. 7C–E,G–K).
Neurotoxic morphological defects revealed in primary neuron culture. (A–M) Phase-contrast images (60×) of CNS neurons from sn3/sn3 (A–K) or wild-type (OreR-C; L–M) wandering third instar larvae, cultured for 3 d.i.v. with the following drugs (scale bar in M for images A–M): (A) tannic acid, 50 μM; (B) 3-3′-diindolylmethane; (C) ginkgolic acid, 50 μM; (D) tegaserod maleate, 50 μM; (E) azadirachtin, 50 μM; (F) 4’-demethylepipodopyllotoxin, 10 μM; (G) usnic acid, 10 μM; (H) atorvastatin, 50 μM; (I) rosuvastatin, 50 μM; (J) lovastatin, 50 μM; (K) pravastatin, 50 μM; (L) atorvastatin, 50 μM; (M) no-drug control. (N–O) The statin-induced BOS defect is modulated by genetic background. Box-plot distributions of the bead density (beads per 100 μm) along neurites of wild-type (+/+; gray) or fascin-deficient mutant (sn3/sn3; black) larval neurons cultured for 3 d.i.v. with pravastatin (N) or rosuvastatin (O) at the indicated concentrations. The median is indicated by the triangle. The top and bottom of the box represent the 75th and 25th percentiles, respectively; the top and bottom ‘whiskers’ represent the 90th and 10th percentiles, respectively. The number of neurons analyzed is indicated in parentheses below each group. At both concentrations, the mutant neurons are more sensitive than the wild-type neurons to the BOS effect of the statins (Mann-Whitney Rank Sum test). n.s., not significant; *P<0.05; **P<0.01; ****P<0.0001. (P–S) BOS beads contain mitochondria (scale bar in P for images P–S). Wild-type larval neurons expressing a GFP-tagged mitochondrial protein cultured for 3 d.i.v. and imaged live by phase-contrast (P,R) or fluorescence microscopy (Q,S). (Grid lines on the dish floor are visible behind the neuron in R.) In the no-drug control cultures (P,Q) GFP-tagged mitochondria are distributed as clusters or linear aggregates scattered throughout the neurite arbor, with a range of sizes and densities. In the atorvastatin-treated cultures (R,S), the GFP-tagged mitochondria are localized almost exclusively to the beads.
Within the NINDS-II library, 22% of the compounds are experimental (i.e. research reagents) rather than marketed pharmaceuticals or drugs that have been tested in clinical trials (supplementary material Table S1, Status column). Experimental compounds were over-represented in the cytotoxic-at-high-dose (54%) and morphological-effects (31%) groups, and were under-represented among the simple fascin-pathway modifiers (16%). It is not surprising that compounds with well-established toxicity, such as the sodium pump inhibitor ouabain and the protein synthesis inhibitors cycloheximide and anisomycin, had overt cytotoxic effects on cultured neurons. Similarly, the ionophore lasalocid, a high-potency fascin-pathway enhancer with cytotoxicity at 50 μM, is a known neurotoxin (Safran et al., 1996). Hence, much of the neurotoxicity among the active compounds was caused by molecules that are not drugs in the strict sense of the repurposing strategy. In many cases, synthetic chemistry can be used to modify such parent compounds to produce analogues that retain the desired activity but have reduced toxicity (e.g. Harrap, 1985).
Statin-induced BOS defect is modified by genetic background
Beads-on-a-string (BOS) is a striking morphological defect that was induced solely by the four statins in the library: atorvastatin, lovastatin, rosuvastatin and pravastatin (Fig. 7H–K). Wild-type neurons are also susceptible to the BOS defect when cultured with any of these four statins (Fig. 7L and data not shown). BOS has two components, ovoid intracellular nodules distributed along the length of the neurites and reduced neurite outgrowth (Fig. 7, compare L with M). These four statins differ in their potency of HMG-CoA reductase inhibition (White, 2002), and this paralleled their potency to induce BOS. Pravastatin and lovastatin induced BOS only at 50 μM, whereas rosuvastatin and atorvastatin did so at both 10 μM and 50 μM (Fig. 7N,O and data not shown). BOS was reversible when the culture media was replaced with drug-free media; the beads resolved and neurite outgrowth was markedly accelerated (data not shown). Using a transgenic GFP-tagged mitochondrial protein (Pilling et al., 2006), we found that the beads contained aggregations of mitochondria (Fig. 7P–S). The appearance of statin-induced beads is consistent with disrupted microtubule transport (Pilling et al., 2006), which would also contribute to reduced neurite outgrowth.
Because statin-induced myopathy is affected by genetic variants (Niemi, 2010), we tested the hypothesis that sensitivity of Drosophila neurons to statin-induced BOS is influenced by genotype. We compared the bead density of wild-type and fascin-deficient singed-mutant neurons in response to incubation with a low-potency (pravastatin) and a high-potency (rosuvastatin) statin. For both drugs at both 10 μM and 50 μM, fascin deficiency significantly enhanced the sensitivity of cultured neurons to statin-induced BOS (Fig. 7N,O).
DISCUSSION
The fascin pathway as a drug-discovery target
Motivated by the importance of fascin function in tumor invasion, as well as in brain development and plasticity, we designed and conducted a screen for pharmacological modifiers of the fascin pathway. By using primary cultured neurons in a cell-based assay, we avoided the disadvantages of immortalized cell lines and made use of the fact that dissociated mutant neurons can reveal cellular phenotypes with a magnitude far greater than that in the intact brain (Kraft et al., 2006; Liu et al., 2007; Chen and Herrup, 2008; Sawallisch et al., 2009). The filagree phenotype of fascin-deficient neurons, which is associated with a marked disruption of actin cytoskeleton organization (Kraft et al., 2006), causes a striking exaggeration of the intrinsic tendency of cultured neurites to turn clockwise (Tamada et al., 2010). This cell-based assay, with its easily observable morphological read-out, allowed us to cast a wide net, targeting the pathway without trying to predict a priori which molecule (fascin or one of its regulators or mediators) would be the best drug target for modulating fascin function. By choosing a well-characterized hypomorphic singed mutation, we were able to screen simultaneously for fascin-pathway blockers and enhancers in cells containing small amounts of wild-type fascin protein that could potentially serve as a drug target. Despite the limited sensitivity of holistic scoring, by screening a diverse collection we identified diverse high-potency compounds in both activity classes.
Fascin-pathway blockers identified in the Drosophila neuronal fascin bioassay are predicted to inhibit tumor invasion. Malignant glioblastoma represents a particularly challenging clinical problem because these cells migrate avidly along many routes beyond the boundaries of the main tumor mass. Chemotherapy with the new anti-proliferation drug temozolomide has increased patient survival, but only by months (Mangiola et al., 2010). It is likely that long-term survival will require a second drug to inhibit glioma dispersion (Giese et al., 2003; Salhia et al., 2006). Hence, we tested seven fascin-pathway blockers in a radial migration assay with two human glioblastoma multiforme cell lines, and found that three of the compounds caused dose-dependent inhibition of migration in both lines (L.L.R., W. S. McDonough, A.K. and M. E. Berens, unpublished results). This is encouraging preliminary cross-species evidence that a subset of compounds identified through the Drosophila fascin bioassay will have predictable effects on the phylogenetically conserved fascin pathway of mammalian cells. We note that, for both RET-dependent thyroid cancer (Das and Cagan, 2010) and fragile X syndrome (Bhogal and Jongens, 2010), Drosophila models have provided evidence of drug efficacy on the road to clinical trials (ClinicalTrials.gov database, NCT00514046 and NCT00965432, respectively).
The tricyclic antidepressant imipramine was a high-potency fascin-pathway blocker in our assay. An intriguing study of imipramine blue reported that this triphenylmethane derivative has anti-invasion properties against malignant glioma cells in vitro and in vivo (Munson et al., 2012). Although that study design was based on inhibition of NADPH oxidase by imipramine blue, glioma cells treated in vitro showed dramatic reorganization of their actin cytoskeleton, with a marked loss of actin bundle-based protrusions and extensions (Munson et al., 2012), consistent with the loss of fascin function (Vignjevic et al., 2006).
What can we infer about molecular mechanisms of action of the fascin-pathway modulators? Some possibilities can be proposed or ruled out based on the nature of the fascin bioassay. First, the molecular diversity of the fascin-pathway modulators argues strongly against a common target for either the enhancers or the blockers. Candidate targets, in addition to fascin per se, include CREB, CREBBP, Stat3, PKCα, Rab35, Rac and myosin X (Hashimoto et al., 2011).
Second, the ‘known’ pharmacology of the active compounds (supplementary material Tables S2, S3) might not be relevant, and this is implicit in the rationale of drug repurposing. For example, because the dissociated neuron cultures are low density by design, with limited cell-cell contact, the effects on neurite trajectory of neuroactive compounds like acetylcholine, adiphenine, baclofen and spiperone cannot be mediated through classical neurotransmission at functional synapses. Rather, their effects are better classified with the morphogenetic phenomena by which neurotransmitters influence brain development (Herlenius and Lagercrantz, 2001; Verney, 2003). Similarly, the antibiotic ceftriaxone, best known for blocking bacterial cell wall synthesis, has neuroprotective activity due to upregulation of glutamate transporter GT1 expression (Rothstein et al., 2005). However, this mechanism is unlikely to explain the activity of ceftriaxone in the fascin bioassay because GT1 is a glial protein, and our cell culture system is neuron-only.
Third, blockers that target the residual wild-type fascin in the singed-mutant neurons most probably do not have the same mechanism as the migrastatin family. This is because Drosophila and mammalian fascin differ at amino acid position 474, which is essential for migrastatin binding (Chen et al., 2010). That site is on the surface of fascin (Jansen et al., 2011), a location with potential disadvantages for drug targeting. By contrast, high-resolution X-ray crystallographic studies identified internal pockets within the two actin-binding sites of fascin that, on the basis of glycerol and polyethylene glycol binding, could accommodate small-molecule blockers (Jansen et al., 2011).
Fourth, the fascin-pathway modulators with additional effects (Table 1) might represent distinct targets. For instance, 16 of the 18 fascin-pathway enhancers with other morphological effects caused reduced neurite outgrowth. Note that very small neurite arbors can manifest the filagree phenotype (Kraft et al., 2006). Therefore, our observations suggest that slowing neurite outgrowth compensates for fascin deficiency, allowing straighter neurite trajectory. Consistent with this, in vitro treatment of fascin-deficient neurons with a physiological hormone that boosts neurite outgrowth increases neurite curvature (Kraft et al., 2006). Although slowing the rate of neurite outgrowth might not be appropriate for therapeutic intervention, the observed relationship between neurite trajectory and outgrowth rate provides an important clue to fascin biology.
Repurposing and structure-activity relationships
Repurposing drug screens are done in the hope that the time to clinical efficacy trials will be considerably shorter than for a new molecular entity and, at least in theory, that retail prices will therefore be lower (Dueñas-González et al., 2008). Nonetheless, SAR investigation continues to be useful for probing the molecular requirements for activity. For example, the inferences from the activity cliffs (Fig. 6) represent SAR hypotheses that can be tested in the fascin bioassay using commercially available compounds that are similar to the original actives. Moreover, the actives can also be the starting point for structural modification for lead optimization, to enhance potency or to enhance blood-brain barrier penetration and CNS distribution (Denora et al., 2009).
Primary neuron culture for neurotoxicity screening in a genetic model organism
Insects lack the post-squalene sterol biosynthetic pathway and thus cannot synthesize cholesterol (Chapman, 1998). Hence, the BOS defect cannot be due to statin effects on cholesterol levels. HMG-CoA reductase inhibition also blocks the synthesis of isoprenoids that serve as membrane anchors for small GTPases such as Ras, Rho and Rac (Cordle et al., 2005; Wang et al., 2008). Isoprenyl membrane anchors are necessary for the activation of these signaling molecules, which serve essential roles in nervous system development and synaptic plasticity (Hall and Lalli, 2010; Samuel and Hynds, 2010; Ye and Carew, 2010). In mammals, statin neurotoxicity is manifest as reduced neurite outgrowth in vitro due to blockade of the isoprenoid pathway (Schulz et al., 2004) and, in the case of sympathetic neurons, specifically by inhibition of RhoA activation (Kim et al., 2009b). Thus, the BOS neurotoxic effect could involve the functional inhibition of small GTPases. We propose that BOS reflects disruption of microtubule-based transport, resulting in large aggregations of organelles, including mitochondria (Fig. 7P–S) and synaptic vesicles (data not shown). The synergistic effect of statin exposure and fascin deficiency highlights the importance of fascin for actin-myosin interactions (Ishikawa et al., 2003; Nagy et al., 2008; Tamada et al., 2010) that link actin and microtubule cytoskeletons (Cao et al., 2004).
BOS stands out from other neurotoxic effects (Fig. 7A–G) because of growing evidence that statins can cause significant, reversible neurocognitive or psychiatric side effects (King et al., 2003; Wagstaff et al., 2003; Golomb et al., 2004; Galatti et al., 2006; Tatley and Savage, 2007; Evans and Golomb, 2009). Recognition and reporting of these statin-caused adverse events, previously under-recognized and under-reported (Golomb et al., 2007), are likely to increase now that the FDA has added confusion and memory loss to statin warning labels (for more details, see http://www.fda.gov/Drugs/DrugSafety/ucm293101.htm). There is some urgency to address statin-induced cognitive side effects. In a randomized clinical trial, statin treatment disrupted learning in patients, who did not perceive that they had developed a cognitive deficit (Muldoon et al., 2004). In addition, there is a push from industry to approve statins for over-the-counter sales (Tinetti, 2008) and for use in children (de Ferranti and Ludwig, 2008). The latter is particularly worrisome because the developing brain might well be at higher risk for sides effect from statins.
There is presently no way to predict which patients will have statin-induced neurocognitive side effects. By contrast, statin-induced myopathy is associated with a single-nucleotide polymorphism in SLCO1B1 (MIM#604843), which encodes a hepatic statin transporter (Niemi, 2010). Our demonstration that fascin modulates BOS (Fig. 7N–O) suggests that Drosophila genetics and neuron culture could be applied systematically to identify candidate genes whose human orthologs control sensitivity to statin-induced cognitive side effects. More generally, tens of thousands of environmental compounds await testing for developmental neurotoxicity (Landrigan, 2010). Primary culture of developing Drosophila neurons could provide a starting point for neurotoxicity screening. There are many potential benefits for screening with an invertebrate in vitro model system (Coecke et al., 2007), especially one so well-suited to examining gene-X-environment interactions.
Brain development and malignant tumors: two sides of the same genetic coin
Fascin is not unique in being essential for brain development, yet also having pathological consequences when overexpressed or overactive in tumor cells. For example, the proto-oncogene MET (MIM#64860) encodes a receptor tyrosine kinase whose activation drives tumor progression and drug resistance (Stella et al., 2010). On the other hand, MET controls neuronal differentiation and is a risk gene for both schizophrenia (Cannon, 2010) and autism, with the autism-associated variants disrupting MET or reducing MET expression (Judson et al., 2011). Similarly, the group I PAK (p21-activated kinase) genes, which regulate cytoskeletal dynamics, are drivers of tumor invasion (Molli et al., 2009). However, they play crucial roles in brain development and synaptic plasticity, with cognitive and behavioral deficits resulting from PAK3 (MIM#300142) loss-of-function mutations (Kreis and Barnier, 2009).
Extrapolating further, one would predict that patients with developmental brain disorders due to loss-of-function mutations should be at lower risk for invasive cancers. A relevant case report describes a child with fragile X syndrome whose untreatable glioma had an indolent course for more than 8 years despite high-grade histopathology (Kalkunte et al., 2007). In fact, this biological duality might explain why patients with schizophrenia, a highly heritable, late-onset brain-development disorder (Karlsgodt et al., 2008), have lower rates of cancer despite engaging in high-risk behaviors such as smoking and drinking (Cannon, 2010). Thus, two seemingly unconnected medical conditions are linked by genes with dual roles in brain development and tumor biology. We propose that the bidirectional screen presented here represents a powerful new approach for simultaneous drug discovery for brain-development disorders and cancer.
METHODS
Drosophila culture and genetics
Fly stocks were maintained at room temperature on corn meal-yeast-agar medium obtained from the University of Arizona fly food facility. Experimental cultures were reared at 25°C and 60–80% relative humidity on corn flour-yeast-agar medium as previously described (Kraft et al., 2006). OreR-C was the wild-type laboratory strain. The sources of singed mutant stocks, the w(z) strain of w1118 (in a CantonS background), and the γ-MB driver/reporter stock (P[Gal4]201Y UAS-lacZ) were as described previously (Kraft et al., 1998; Kraft et al., 2006). The neuron-specific driver line elav-Gal4C155 (Lin and Goodman, 1994) was provided by Mani Ramaswami (then at University of Arizona, Tucson, AZ). The GAL4-responsive mitochondria-targeted GFP transgenic line w1118; P(UAS-mitoGFP) (Pilling et al., 2006) was provided by Konrad Zinsmaier (University of Arizona, Tucson, AZ). To identify mitochondria in living neurons, we generated elav-GAL4C155/w1118; P(UAS-mitoGFP)/+ progeny larvae.
γ-MB neurons with wild-type fascin/singed function were identified in vitro by crossing the γ-MB driver/reporter line to w(z) to generate w sn+/ w sn+; P[Gal4]201Y UAS-lacZ / +. To study the effect of fascin/singed mutation severity on the filagree phenotype of γ-MB neurons, we performed three crosses in parallel using the 201Y driver to produce female larval progeny of the desired genotypes. Crossing y w sn3/ Y; P[Gal4]201Y UAS-lacZ males to sn3/sn3 virgin females yielded y w sn3/sn3; P[Gal4]201Y UAS-lacZ / + daughters. In(1) dl-49, snX2/Y males were crossed either to y w sn3; P[Gal4]201Y UAS-lacZ females to produce y w sn3/snX2; P[Gal4]201Y UAS-lacZ/+ daughters, or to w; P[Gal4]201Y UAS-lacZ females to produce snX2/+; P[Gal4]201Y UAS-lacZ/+ daughters.
Chemical informatics
Scaffold analysis
The analysis of the scaffold distribution described in this work is based upon specific molecular chemotypes (Medina-Franco et al., 2006) called cyclic systems, originally defined by Xu and Johnson (Xu and Johnson, 2002). This approach is consistent with general guidelines for defining molecular scaffolds (Langdon et al., 2010). A cyclic system consists of a set of rings and the chains of atoms that link them to one another. To obtain the scaffold, all substituent groups are removed from the rings and the linkers except endocyclic carbonyls and imines. Heteroatoms are retained and all hydrogen atoms attached to them are considered side chains and therefore removed. The Meqi program (version 2.41), originally developed by Xu and Johnson (Xu and Johnson, 2002), was obtained from Mark Johnson (Pannanugget Consulting, Kalamazoo, MI) and used to determine the specific cyclic system (described by a five-character alphanumeric chemotype identifier) present in each compound in the NINDS-II library. The Scaffold Recovery Curve was generated as previously described (Lipkus et al., 2008; Medina-Franco et al., 2009).
Molecular similarity and library diversity
Molecular similarity refers to the relation between pairs of molecules, whereas molecular diversity is based on pair-wise relations of a population of molecules (e.g. a compound collection). Molecular fingerprints of each compound were used to measure all pair-wise similarities (Willett et al., 1998; Maggiora and Shanmugasundaram, 2011). In general, this approach involves designating the presence or absence of specific structural features, which can be considered molecular fragments. This information is encoded in ‘bit vectors’, each of whose components are assigned a value of ‘1’ if the fragment is present at least once in the molecule and ‘0’ if it is not. Thus, multiple occurrences of specific fragments are not explicitly counted. The work presented here employed MACCS key fingerprints (Anderson, 1984), the current version of which represents 166 molecular fragments and related structural data. MACCS Structural Keys software, now available through Accelrys, Inc. (San Diego, CA), was used to determine the fingerprints of all compounds in the library. The molecular fingerprints were then used to evaluate the Tanimoto similarities using the program Molecular Operating Environment (MOE, v.2010), available from the Chemical Computing Group, Inc. (Montreal, Quebec, Canada; www.chemcomp.com). The Tanimoto coefficient between two molecules, A and B, is given by T(A,B)=c/(a+b–c), where a is the number of fragments present in molecule A, b is the number of fragments present in molecule B and c is the number of fragments common to both A and B. Details of this procedure, which is widely used in chemical informatics applications, are available (Willett et al., 1998; Maggiora and Shanmugasundaram, 2011).
Principal component analysis of chemical-space distribution
Molecular similarity measures do not provide explicit information about the chemical space that molecules inhabit. To obtain a graphical representation of this space we created 3D plots using principal component analysis (PCA) (Jolliffe, 2002) based on the Tanimoto similarity matrix associated with the molecules in the NINDS-II collection. In the rectangular data matrix usually employed in PCA, the rows correspond to the objects under study and the columns to the descriptors, features or parameters used to describe the objects. By contrast, the similarity matrix for a set of n molecules is an n×n square symmetric array of Tanimoto similarity values where the rows and columns both correspond to the molecules under study. To obtain the principal components, the similarity matrix was treated as if it were a typical rectangular data matrix and the usual principal component procedure was applied. The columns were mean-centered, matrix multiplication was applied, and the product matrix was diagonalized yielding the eigenvectors and eigenvalues. The molecules (the objects in this study) were then plotted in the principal-component coordinate system using the first two or three eigenvectors or the principal components. Typically, the first three principal components represent most (65–75%) of the variance of the data. In this study, they represent ∼67% of the variance, whereas the first two represent ∼55% of the variance. The principal component data and plots were obtained using the Spotfire® dynamic data visualization program (v.9.1.2) from TIBCO Software, Inc. (Somerville, MA; spotfire.tibco.com).
Primary neuronal cell culture preparation and evaluation
Primary neuronal cell cultures from the central nervous system (CNS) of wandering third instar larvae were prepared as previously described (Kraft et al., 1998; Kraft et al., 2006) with minor modifications. This method generates pure-neuron cultures within several days of incubation. Enzymatic treatment of the dissected tissue prior to dissociation was done with a commercial blend of purified collagenase I, collagenase II and the neutral protease dispase (Liberase Blendzyme 1, #1988409; Roche, Indianapolis, IN) at a concentration of 72 μg/ml in Rinaldini’s saline. Schneider’s Insect Medium (#11720; Invitrogen, San Diego, CA) was supplemented with 10% fetal bovine serum of US origin (#SH30071; lot #APC20860; Hyclone, Logan, UT) and 50 μg/ml bovine insulin (#I-6634; Sigma, St Louis, MO). For experiments analyzed by random-neuron sampling and quantification of morphometric parameters, the floor of the culture dish was made with a photoetched (‘gridded’) glass cover slip (#1916-91818; Bellco, Vineland, NJ). For cultures to be evaluated by holistic scoring (see below), a round German #1 glass cover slip (#1943-10012; Bellco) was used. In either case, the resulting 8-mm-diameter wells were coated with Conconavalin A (#C-2010; Sigma) and laminin (#354232; BD Biosciences, Franklin Lakes, NJ).
Unless otherwise specified, CNS of female larvae were explanted by microdissection and the entire CNS was dissociated, with the cells distributed into six culture dishes, followed by incubation at 25°C for 3 days (70–80 hours). After this time, each culture dish typically contained 1500–2000 neurons bearing scorable neurite arbors. The phenotypic severity of the filagree defect of fascin-deficient singed-mutant neurons is similar whether they are from late larval or early pupal CNS (Kraft et al., 2006). We chose the larval stage for the screen because the CNS dissection is far easier. Neuronal cultures prepared from individual CNS regions (optic lobes, brain, ventral ganglion) of sn3/sn3 larvae had indistinguishable filagree severity (data not shown).
Although neurite curvature can be quantified using image-analysis software (Kraft et al., 2006; Narro et al., 2007), the method requires immunostaining, random sampling, image acquisition and image processing. Alternatively, classification of neurite-curvature phenotypes can be performed much more quickly by holistic scoring of all neurons in a culture well (generally 1500–2000 neurons distributed over a ∼50 mm2 area). Classification of wild-type (+/+), intermediate-filagree (sn3/sn3), and severe-filagree (snX2/Y) cultures using phase-contrast microscopy of live neurons is 100% accurate when performed by trained observers. Training was done using non-drug-treated cultures of known genotype after 3 d.i.v. Training consisted of (i) a learning phase, with repeated observations by phase-contrast microscopy of live cultured neurons with known genotypes (+/+, sn3/sn3, snX2/Y), exhibiting wild-type, moderate filagree, or severe filagree phenotypes, respectively; and (ii) a testing phase, using coded dishes or images of individual microscopic fields. Training continued until dishes were assigned the correct genotype with 100% accuracy, and individual fields with ∼90% accuracy (see Kraft et al., 2006).
To guide the selection of compound concentrations for screening, we tested the effect of the solvent dimethyl sulfoxide (DMSO, USP grade; #D-2438, Sigma) on neurite-arbor morphogenesis, selecting genetically marked γ-MB neurons (Kraft et al., 1998; Kraft et al., 2006) to minimize biological variation. There was no systematic concentration-dependent effect of DMSO on neurite curvature of fascin-deficient sn3/sn3 larval γ-neurons (supplementary material Fig. S1). However, at a concentration equivalent to a 100-μM drug treatment, DMSO significantly reduced neurite arbor size (data not shown). We therefore tested the compounds at 10 μM and 50 μM.
Screening the compound collection
The NINDS Custom Collection II (NINDS-II) library of 1040 biologically active compounds was purchased from Microsource Discovery Systems (Gaylordsville, CT). The collection was provided in 13 microplate-formatted storage racks, each with ten columns of eight polypropylene storage tubes sealed with a single non-separable CapStrip (#4415; Matrix Technology, Hudson, NH). Each compound was provided as a 10-mM solution in DMSO and uniquely identified by its rack number, row and column location. Each column of eight compounds was screened at the same time. The racks were stored at −80°C and each rack was removed only long enough to retrieve or replace a single column of compounds. Because the melting point for DMSO is 19°C, each column of tubes held together by a CapStrip was floated upright in a room-temperature water bath for 15 minutes to facilitate thawing the contents before use.
For each experiment, primary cultures were prepared from CNS of wild-type (OreR-C) and homozygous sn3 mutant larvae as described above, with each dissociated CNS preparation distributed into six 8-mm-diameter culture wells. Two of the sn3 mutant dishes served as no-drug controls for the baseline filagree phenotype. OreR-C cultures served as no-drug controls for wild-type neurite outgrowth and arbor morphology, and to verify that the culture medium and dish coating supported typical neuron survival and morphology. The cell suspension from each CNS was distributed equally into six culture dishes and the total volume in each dish was brought to 100 μl with culture medium. While the cells were allowed to settle and adhere for 2 hours at 25°C, eight compounds were each diluted in culture medium as needed for testing at final concentrations of 10 and 50 μM (final volume 1 ml). Liquid handling was facilitated by using an eight-channel electronic pipettor (Matrix) and barrier tips. Each dish of cells was supplemented with 900 μl of culture medium for final drug concentrations of 0, 10, 50 or 100 μM. The dishes were sealed with Parafilm-M® to minimize evaporation and then incubated with humidification at 25°C.
To minimize microbial contamination risk during the compound screening, we adopted a number of exceptional measures. A dissecting microscope and fiber-optic light source were moved into a dedicated tissue culture room that was decontaminated nightly with germicidal ultraviolet light. Dissections and dissociation were performed wearing sterile surgical gloves and a laboratory coat. The drug-containing racks were handled and the tubes thawed by a second person. The CapStrip was replaced after the column of drug vials was accessed.
After 3 d.i.v., the live cultured neurons were evaluated independently by two or three trained observers using phase-contrast microscopy on a Nikon Diaphot 300 inverted microscope using a 40× objective (numerical aperture, 0.7). Holistic scoring was performed by systematic scanning of each entire dish to assess the overall effect of each drug treatment on the severity of the moderate filagree phenotype of the sn3/sn3 mutant neurons. A compound at a given dose was scored as decreasing the filagree phenotype (fascin-pathway enhancer) or increasing the filagree phenotype (fascin-pathway blocker) when a clear preponderance of neurons in the dish were observed to exhibit an unambiguous shift of the phenotype in the same direction. If no change was evident or there was more than the usual variability in the filagree phenotype, the score was recorded as ‘no effect’. In addition to scoring the filagree phenotype, other effects on neuronal morphology were noted, again only when a significant majority of neurons in the culture displayed a consistent feature. An ‘other effect’ scoring scheme evolved over the course of the drug screening. A minimum of 3–5 minutes was required to examine each dish and additional time was necessary when an effect was noticed and its distribution throughout the dish was evaluated.
Scores for each compound at the low and high dose were recorded by each observer, after which the outcomes between observers were compared and discussed. In most cases, the scores were concordant. In cases of disagreement, the individuals would re-assess the cultures together and come to a consensus score. Drug-induced changes in filagree (increased or decreased severity) or in other aspects of neuronal morphology were documented by image acquisition of representative examples using a 60× oil-immersion objective (numerical aperture, 1.4) and a Hamamatsu ORCA285 digital camera with HCImage software (Hamamatsu Photonic Systems, Bridgewater, NJ).
Experiments were repeated for compounds causing filagree modification or other effects on neuron morphology so that the scores could be confirmed. In cases where the scores in the second test were not consistent with the original experiment, the compounds were repeated a third time and the majority score from the three experiments was taken as the definitive outcome. Compounds were also retested in instances when low neuron density or microbial contamination had compromised scoring. Once a coded compound was successfully evaluated, its name was revealed. A survey of published biomedical literature and publicly available databases was then conducted to collect pharmacology, toxicology and usage information.
Immunostaining and morphometric analysis of 2D neuron images
Statin compounds were re-tested on wild-type (OreR-C) neurons, on neurons expressing a mitochondria-targeted GFP and on singed-mutant neurons. Reciprocal drug wash-out experiments were performed on wild-type neurons in which the culture medium was replaced after 1 d.i.v. after a random sample of 25 neurons was imaged by phase-contrast microscopy. Those neurons were re-imaged live at 2 and 3 d.i.v. To quantify BOS severity, bead counting was based on phase-contrast images of live neurons. A ‘bead’ was defined as a phase-dark round or ovoid swelling within a neurite with a diameter greater than that of the flanking neurite on both sides. Because neurite expansions are often seen at branch points and neurite tips, potential beads at these locations were difficult to assess and therefore not counted.
Neurite length and curvature measurements were based on images of immunofluorescently labeled neurons. Fixation and immunostaining were performed after 3 d.i.v. as previously described (Kraft et al., 1998; Kraft et al., 2006). To visualize neuronal membranes with high resolution, the Drosophila Nervana 2 protein (Sun and Salvaterra, 1995) was labeled with a polyclonal goat anti-horseradish peroxidase antiserum (Sigma) at 1:500 and detected with an Alexa-Fluor®-488-conjugated donkey anti-goat antiserum (Invitrogen) at 1:500. To visualize β-galactosidase expression in genetically marked γ-MB neurons, a preabsorbed polyclonal rabbit anti-β-galactosidase antiserum (Cappel, West Chester, PA) was used at 1:5000 and detected with a Cy3-conjugated donkey anti-rabbit antiserum (Jackson ImmunoResearch, West Grove, PA) at 1:500. Labeled neurons were identified on the Nikon Diaphot 300 inverted microscope with the 60× oil-immersion objective using epifluorescence illumination. The Alexa Fluor 488 fluorescent signal was detected with filter cube Chroma #41001 (exciter 460-500 nm, dichroic 505 nm, band-pass emitter 510–560 nm) and the Cy3 signal with Nikon G-2A filter cube (exciter 510–560 nm, dichroic 580 nm, long-pass emitter 590 nm). Images were collected with the Hamamatsu ORCA285 digital camera and HCImage software.
Acquisition of images of γ-MB neurons and their generic neuron neighbors was performed as described (Kraft et al., 2006). The nearest-neighbor neurons constitute a random sample of the heterogeneous non-γ CNS neuron types in the culture. For quantitative study of neurons in non-genetically marked cultures, a rigorous method of random sampling was done by following the alphanumeric grid of the dish in a systematic manner, up and down a ‘staircase’, acquiring images of neurons along the ‘steps’ until at least 50 neurons were imaged.
Computer-assisted analysis of the fluorescent neuron images was performed using NeuronMetrics™ software, which converts each neurite arbor into a one-pixel-wide skeletal representation from which morphometric parameters, including length, are calculated (Narro et al., 2007). For each neuron, the number of beads per 100 μm (bead density) was calculated from the total bead count and total neurite length. The nonparametric Mann-Whitney rank sum test was used to compare the distributions of bead-density values between populations of the same genotype treated with different drug concentrations, and between populations with different genotypes treated with the same drug concentration. The statistical tests and graphical representation of the data were performed using SigmaStat and SigmaPlot (Systat Software, San Jose, CA). The neurite-arbor skeletons were also used for fully automated determination of mean neurite curvature for each neuron (Kraft et al., 2006). Distributions of neurite curvature were plotted as histograms with soft binning, to reduce the effect of discretization error, and a normal-distribution curve was plotted over the data. In soft binning, each value is considered the mean of a normal distribution, with the standard deviation being set to the bin size. We then distribute the value of that point to each of the bins in proportion to the value of the normal distribution at the bin centers. Comparisons of curvature distributions between groups were made using the Welch’s t-test (Myers and Well, 1991), which accounts for different variances and sample sizes.
Acknowledgements
The authors thank Mindy Escobar, Benjamin Reddoch, Michelle O’Shea, David Andrew and Katie Olson for assistance in the early stages of the drug screen. We are also grateful to Anna Burns and Lynne Oland for expert oversight of the Department of Neuroscience tissue culture facility, to Konrad Zinsmaier and Mani Ramaswami for Drosophila stocks, to Andrew Predoehl for help running the neurite-curvature software, to Nina Ossanna and Laurence Hurley for advice about drug discovery, and to Mark Johnson (Pannanugget Consulting) for providing the Meqi program.
Footnotes
-
FUNDING
The project was funded by awards to L.L.R. from the National Institutes of Health [grant number NS055774], to L.L.R. and G.M.M. from the Center for Insect Science (University of Arizona), and to J.L.M.-F. from the State of Florida and the Menopause & Women’s Health Research Center (Port St Lucie, FL). C.B. was funded by NIGMS (Minority Access to Research Careers).
-
COMPETING INTERESTS
G.M.M. was previously employed by Upjohn and Pharmacia.
-
AUTHOR CONTRIBUTIONS
L.L.R. conceived and developed the experimental concepts and directed the project; L.L.R. and R.K. designed experiments; R.K., A.K., M.L.O. and C.B. performed laboratory experiments; R.K., K.B., G.M.M., J.M.-F., F.L.-V., M.L.O., C.B. and L.L.R. contributed data analysis; and L.L.R., R.K., K.B., G.M.M., F.L.-V. and J.M.-F. prepared and edited the manuscript and figures.
-
SUPPLEMENTARY MATERIAL
Supplementary material for this article is available at http://dmm.biologists.org/lookup/suppl/doi:10.1242/dmm.008243/-/DC1
- Received January 12, 2012.
- Accepted July 31, 2012.
- © 2012. Published by The Company of Biologists Ltd
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial Share Alike License (http://creativecommons.org/licenses/by-nc-sa/3.0), which permits unrestricted non-commercial use, distribution and reproduction in any medium provided that the original work is properly cited and all further distributions of the work or adaptation are subject to the same Creative Commons License terms.
REFERENCES
Article navigation
Related articles
Cited by...
More in this TOC section
Similar articles
Other journals from The Company of Biologists
Editor’s choice – DJ-1 plays a positive nonmitochondrial role in endocytosis in Dictyostelium cells
The Dictyostelium homologue of the Parkinson's disease-associated protein DJ-1 is located in the cytosol. In this Editor's choice article, Paul R. Fisher and colleagues reveal that loss of DJ-1 causes cytopathological defects in endocytic and autophagic cell death pathways, but stimulates respiration by functionally normal mitochondrial respiratory complexes.
Review – Central and peripheral circadian clocks and their role in Alzheimer's disease
Damian C. Crowther and colleagues summarise in a new Review how studies in Drosophila and mouse models accelerate our understanding of the links between circadian rhythm dysfunctions and age-related neurodegenerative disorders, especially Alzheimer's disease.
This article features in a special collection on circadian rhythm research from The Company of Biologists, in honour of the 2017 Nobel Prize in Physiology or Medicine.
Featured article – Microtubules protection through acetylation
In a new Research article, Yong Q. Zhang and co-workers show that increased acetylation of microtubules rescues human tau overexpression-induced toxicity in Drosophila muscles and neurons.
Organoid special collection from The Company of Biologists
Organoids, artificially grown masses of cells or tissue that resembles an organ, sit at the interface of developmental, cell and translational biology. To reflect the cross-disciplinary importance of the field, The Company of Biologists presents this special cross-journal collection of articles on organoids published in DMM and our sister journals Development, Journal of Cell Science and Biology Open.