Discovery and Characterization of a Small Molecule Inhibitor of the PDZ Domain of Dishevelled*

Dishevelled (Dvl) is an essential protein in the Wnt signaling pathways; it uses its PDZ domain to transduce the Wnt signals from the membrane receptor Frizzled to downstream components. Here, we report identifying a drug-like small molecule compound through structure-based ligand screening and NMR spectroscopy and show the compound to interact at low micromolar affinity with the PDZ domain of Dvl. In a Xenopus testing system, the compound could permeate the cell membrane and block the Wnt signaling pathways. In addition, the compound inhibited Wnt signaling and reduced the levels of apoptosis in the hyaloid vessels of eye. Moreover, this compound also suppressed the growth of prostate cancer PC-3 cells. These biological effects suggest that by blocking the PDZ domain of Dvl, the compound identified in our studies effectively inhibits the Wnt signaling and thus provides a useful tool for studies dissecting the Wnt signaling pathways.

Dishevelled (Dvl) is an essential protein in the Wnt signaling pathways; it uses its PDZ domain to transduce the Wnt signals from the membrane receptor Frizzled to downstream components. Here, we report identifying a drug-like small molecule compound through structure-based ligand screening and NMR spectroscopy and show the compound to interact at low micromolar affinity with the PDZ domain of Dvl. In a Xenopus testing system, the compound could permeate the cell membrane and block the Wnt signaling pathways. In addition, the compound inhibited Wnt signaling and reduced the levels of apoptosis in the hyaloid vessels of eye. Moreover, this compound also suppressed the growth of prostate cancer PC-3 cells. These biological effects suggest that by blocking the PDZ domain of Dvl, the compound identified in our studies effectively inhibits the Wnt signaling and thus provides a useful tool for studies dissecting the Wnt signaling pathways.
The Wnt signaling pathways are regulated by a family of secreted Wnt glycoproteins. The canonical Wnt pathway, which is highly conserved, is best understood. In this pathway, Wnt molecules interact with the seven-transmembrane Frizzled (Fz) 2 proteins (1) by binding to an N-terminal cysteinerich-domain (2). The signal is then transduced into the cell through an internal sequence of Fz, C-terminal to the seventh transmembrane domain, which binds directly to the PDZ (postsynaptic density-95/discs large/zonula occludens-1) domain of the cytoplasmic protein Dishevelled (Dvl) (3). Dvl then transduces the Wnt signals to downstream components (4). Three Dvl homologs (Dvl-1, -2, and -3) have been identified in humans; all are expressed in both embryonic and adult tis-sues, including brain, heart, lung, kidney, skeletal muscle, and others (4). Up-regulation and overexpression of Dvl proteins have been reported in many cancers, including those of breast, colon, prostate, mesothelium, and lung (non-small cell) (5)(6)(7)(8).
The Dvl protein is made up of three conserved domains: an N-terminal DIX domain, a central PDZ domain, and a C-terminal DEP domain (9). The central PDZ domain is of particular interest because of its interaction with Fz and other Wnt pathway proteins (3,10). The direct interaction between the PDZ domain and Fz peptides is relatively weak, and other factors may play a role to ensure the communication between the two molecules (3). For example, several studies suggest that the DEP domain of Dvl has a membrane-targeting function that may facilitate PDZ-Fz interaction (11)(12)(13)(14). However, the weak PDZ-Fz interaction provides an opportunity to block Wnt signaling at the Dvl level by using a small molecule inhibitor. An earlier study in our laboratories used an NMR-assisted virtual ligand screening approach to identify a peptide mimic that can bind to the Dvl PDZ domain (15). We have now used an improved algorithm to conduct an additional structure-based virtual screen of the PDZ domain of Dvl and have discovered a group of drug-like compounds that bind to the PDZ domain with moderate to low micromolar affinity. One of these compounds effectively blocked Wnt signaling in vivo and reduced the growth rate of a prostate cancer cell line.

EXPERIMENTAL PROCEDURES
Ligand and Receptor Preparation-All compound data bases were obtained from the NCI (National Institutes of Health), Chemical Diversity Inc. (ChemDiv, San Diego, CA), or Sigma-Aldrich. Three-dimensional coordinates for all compounds were generated by using the Optive Research Concord program (Tripos Inc., St. Louis, MO) and stored in the Sybyl mol2 format. The Unity module of the Sybyl software package (Tripos Inc.) was used to select the compounds in the three-dimensional small molecule data base that matched the known ligand of the Dvl PDZ domain. The first Unity query was based on three-dimensional distance constraints determined by analyzing the structure of the Dapper (Dpr) peptide-PDZ complex (10), which was obtained from the Protein Data Bank (PDB code: 1L6O). The atoms of the bound Dpr peptide within hydrogen-bonding distance of suitable H-bond acceptors and donors on the backbone of the ␤B sheet of the PDZ domain were selected.
Docking and Scoring-The compounds selected from the Unity queries were docked to the Dpr peptide binding site of Dvl PDZ domain by using the FlexX module of the Sybyl software package. For the docking analysis, the receptor site was first defined as all of the residues in the PDZ domain within 7.0 Å of the Dpr peptide, and the core site (i.e. where the core fragment is to be placed during docking) was defined as all residues within 10.0 Å of the Thr (Ϫ2) in the Dpr peptide. In later stages of screening, the core site was redefined as a much smaller area made up of residues Ile-267, Ser-268, Ile-269, Leu-324, Arg-325, and Val-328, which make up the hydrophobic groove. This refined definition helped to prevent docked conformations that were likely to be incorrect. For example, when the larger core site was used, some docked compounds were almost entirely exposed to the solvent, although most of the molecules are hydrophobic. As many as 30 docked conformations were generated for each compound by FlexX. Docked conformations were scored by using the five consensus score functions included in the FlexX software. The Dpr peptide was used to calibrate docking parameters and to determine high (negative) values for each individual scoring function. The backbone root mean square deviation between the crystal structure conformation and the docked conformation of Dpr was 2.86 Å.
Protein Preparation-Different forms of the Dvl PDZ domain protein were synthesized by using an Escherichia coli system as described (3,15,16). To prepare TMR-labeled PDZ domain, a cysteine located outside of the binding site was mutated to alanine to increase solubility of the PDZ domain. A construct without a T328C mutation was made for use in fluorescence spectroscopy experiments. The fluorescent label 2-((5(6)-tetramethylrhodamine)carboxyamino)ethyl methanethiosulfonate (TMR) was covalently bound to this Cys residue, which is the only Cys in this PDZ domain. The PDZ solution was dialyzed overnight at 4°C against 100 mM potassium phosphate buffer at pH 7.5 to remove dithiothreitol, which had been added to prevent disulfide bond formation. A 10-fold excess of TMR dissolved in DMSO was added dropwise to the PDZ solution while stirring. After 2 h at room temperature, the unbound TMR was removed from the solution by dialysis against 100 mM potassium phosphate buffer, pH 7.5, at 4°C. NMR Spectroscopy-All NMR studies used either a 600-MHz Varian INOVA spectrometer or a 600-MHz Bruker Avance spectrometer. 15 N-labeled PDZ samples were prepared at a concentration of 0.3 mM in 100 mM KHPO 4 , 0.5 mM EDTA, and 10% D 2 O at pH 7.5. The compounds were dissolved either in the same buffer as the PDZ domain or in DMSO, depending on the aqueous solubility of the individual compound. Titration was carried out by adding small amounts of the compound to the PDZ domain and taking 15 N-HSQC spectra of the mixture. Compound concentrations varied from 0.3 to 6.0 mM during the course of these titrations. All NMR spectra were processed with NMRPipe (17) software and analyzed with the Sparky program (18).
Fluorescence Spectroscopy-All fluorescence measurements were obtained by using a Jobin-Yvon Fluorolog-3 spectroflu-orometer (HORIBA Jobin-Yvon Inc., Edison, NJ) (15). Fluorescence anisotropy measurements of TMR-labeled PDZ domain were obtained for binding affinity calculations. Titrations of compounds to the solution of TMR-labeled PDZ domain were performed at 25°C in 100 mM KHPO 4 , 0.5 mM EDTA buffer, at pH 7.5. The excitation wavelength for TMR-labeled PDZ domain was 551 nm with an entrance slit width of 5 nm. The maximum fluorescence emission wavelength was 578 nm with an exit slit width of 5 nm. The compounds were prepared to a concentration of ϳ1.0 -20.0 mM in the same buffer as that used for the protein. During titration, the range of concentrations of the compound was ϳ100 nM-1.0 mM, depending on binding affinity for the particular compound. The anisotropy data were analyzed by fitting the data to the standard ligand binding curve in the program Prism (GraphPad Software Inc., San Diego, CA). Best-fit curves were obtained by using a global, non-linear regression model that assumed that the law of mass action was followed. Although changes of fluorescence anisotropy due to the compounds binding were small, we were able to obtain binding affinity values for the small molecule compounds. Under the same conditions, fluorescence polarization measurements of the ROX-labeled Dpr peptide (ROX-N-butyric-SGSLKMTTV-COOH) were also performed. The excitation wavelength for the ROX-labeled Dpr peptide was 578 nm with an entrance slit width of 5 nm. The maximum fluorescence emission wavelength was 605 nm with an exit slit width of 5 nm. By titrating the PDZ domain into 50 nM ROX-Dpr peptide solution in the absence and presence of 6 M compound 3289-8625, respectively, we obtained the binding affinity (K d ) of Dpr and the apparent binding affinity (K d app) of the peptide in the presence of 3289-8625. The competition binding constant (K i ) between PDZ and 3289-8625 was calculated by the equation Xenopus Embryo Studies-Xenopus eggs were obtained from females that had been injected with 500 IU of human chorionic gonadotropin (Sigma-Aldrich) and had been artificially fertilized. Synthesis and microinjection of mRNAs were carried out as described previously (15,19). Briefly, for the luciferase assay, the siamois promoter-driven reporter DNA construct (3, 15) (400 pg) was injected alone or co-injected with Wnt3A mRNA (1 pg) into the animal pole region at the two-cell stage. Injected embryos were cultured in the absence or presence of the compound 3289-8625 at different concentrations, and animal cap explants were dissected at late blastula stage. For secondary axis assay, Wnt3A mRNA (1 pg) was injected in the ventro-vegetal region at the four-cell stage, and injected embryos were cultured in the absence or presence of the compound 3289-8625 until larval stage.
Cell Proliferation Assay-The effect of the compound 3289-8625 on the PC-3 cells was examined by using the MTT cell proliferation assay (Promega, Madison, WI). In the experiments, PC-3 cells were plated in 24-well culture plates at a density of 1 ϫ 10 4 cells/well and cultured overnight in Dulbecco's modified Eagle's medium containing 10% fetal bovine serum plus penicillin and streptomycin at 37°C. On the second day, compound 3289-8625 was added to the cells at different concentrations. Before being added to the cells, the compound was dissolved in DMSO at 200 mM. While adding the compound to the cells, the same amounts of DMSO were added to the control cells. After 68 -69 h of treatment, 10% volume of MTT stock solution (5 mg/ml) was added to the cell culture. The cells were incubated at 37°C until the total treatment time reached 72 h. The converted dye was then solubilized, and the absorbance was measured at 570 nm. Three sets of independent experiments were performed, and each data point was normalized against the control cells.
The ␤-Catenin Level in PC-3-PC-3 cells were seeded in 100-mm tissue culture dish, cultured in completed Dulbecco's modified Eagle's medium overnight, and then treated with 3289-8625 compound (final concentration 80 M) or DMSO vehicle for 72 h. For the extraction of membrane and cytosolic proteins (20,21), cells were collected in TES suspension buffer and homogenized on ice. The lysate was centrifuged for 10 min at 500 ϫ g. The crude supernatant was then fractionated at 100,000 ϫ g for 90 min at 4°C to generate a supernatant or cytosolic fraction and a membrane-rich pellet fraction. The membrane-enriched pellet was dissolved in phosphate-buffered saline buffer containing 1% Triton X-100 and 1% Nonidet P-40. Equal volumes of Laemmli buffer were added to protein solutions, and the samples were boiled for 5 min. The proteins were then separated by 10% SDS-PAGE electrophoresis under reducing condition, transferred to nitrocellulose membranes, blocked with 5% nonfat dry milk in phosphate-buffered saline with Tween, probed with antibody against ␤-catenin (Santa Cruz Biotechnology Inc., Santa Cruz, CA) and horseradish peroxidase-conjugated anti-rabbit IgG secondary antibody (Sigma-Aldrich), and visualized by a chemiluminescence detection system (PerkinElmer Life Sciences). Membranes were then stripped and reprobed with antibody against integrin ␣3 protein (22) as equal loading control.

RESULTS
Screening of Dvl PDZ Binding Compounds-To identify additional scaffolds of Dvl PDZ domain inhibitors, we carried out several rounds of new computational screening on the basis of our earlier studies (15,23,24). To identify possible PDZ binding compounds, we first used the UNITY module in the Sybyl software (Tripos Inc.) to screen data bases of drug-like compounds from the NCI (National Institutes of Health), ChemDiv, and Sigma-Aldrich. Hits returned from these searches were docked to the protein receptor site, and the conformations of the complexes were scored by using the FlexX module in Sybyl. The overall procedure of the computational screening was similar to that used in our earlier work (15), but many details were further refined in this study.
The first Unity query was based on three-dimensional distance constraints determined by analyzing the structure of the Dpr peptide-PDZ complex (10); we chose the ligand-based query approach, which was not used in our earlier study (15). The atoms of the bound Dpr peptide within hydrogen-bonding distance of suitable H-bond acceptors and donors on the backbone of the ␤B sheet of the PDZ domain were selected. Distance and angle constraints between those atoms were used to run a three-dimensional flexible screen of the NCI (National Institutes of Health) and ChemDiv data bases; in addition, the Sigma-Aldrich data bases were searched by using the Sigma-Al-drich online two-dimensional search utility. The initial Unity query generated a list of several thousand compounds as potential hits.
The program FlexX (25) was then used in the docking studies. This program uses five scoring functions to evaluate docking results: the standard FlexX scoring function F_score (26 -28), the Chemscore function (29), the knowledge-based proton motive force (PMF) score (based entirely on protein-ligand atom pairs and their distances (30)), the G_score function (calculated from ligand-protein atom pair interactions based on values from the Tripos force field (31)), and the D_score (based on the atom charges and the Van der Waals interactions between the ligand and protein (32)). Data from our earlier studies (15) and from the additional experiments in this study revealed that the F_score and Chemscore functions were the most reliable in predicting binding to the PDZ domain. Therefore, in this study, we used only the F_score and Chemscore rather than the consensus score of all five scoring functions. The G_score function was used as an additional reference because in our experience, G_score can help to identify internal steric hindrance in a docked conformation. Approximately 50 high scoring compounds were then obtained from the NCI (National Institutes of Health) and from Sigma-Aldrich and further screened by using NMR spectroscopy. Before the NMR studies, the FlexX docked complexes of all selected compounds were visually inspected to confirm that the ligands were in the peptide binding groove of the PDZ domain and that there was no internal steric hindrance.
To validate docking results, we performed 1 H-15 N correlated NMR spectroscopy. We obtained the 15 N-HSQC spectra (15) by titrating various concentrations of the small molecules into samples of 15 N-labeled mDvl1 PDZ domain. Examination of the spectra for chemical shift perturbations (33) revealed several small molecules that bound to the conventional C-terminal peptide binding groove of the PDZ domain; many of those NMR-confirmed hits showed a butyric acid substructure at one end. On the basis of this butyric acid group, we used the program Unity to perform a substructure search of the NCI (National Institutes of Health), Sigma-Aldrich, and ChemDiv compound data bases. The program FlexX was then used to dock and score the hits returned from the Unity screen. Several of the highest scoring compounds in the FlexX docking had a similar core structure that was predicted by FlexX to bind in the traditional binding groove of the PDZ domain. Those compounds were obtained from ChemDiv, Inc. and were tested by NMR spectroscopy.
Characterization of Compound 3289-8625, a PDZ Domain Inhibitor-In the NMR experiments, compound 3289-5066 and compound 3289-8625 (Fig. 1) displayed the most significant chemical shift perturbations when titrated into the solution of 15 N-labled Dvl PDZ domain. Their chemical shift patterns were generally similar and closely resembled that seen with Dpr and Fz binding (3). Compound 3289-8625 bound more strongly to the Dvl PDZ domain, as judged by the chemical shift perturbations due to ligand binding (Fig. 2). When the weighted sums of the chemical shift data were plotted onto a tube representation of the mDvl1 backbone, compound 3289-the ␣B helix (Fig. 2). To further assess the binding characteristics of the molecules identified in the virtual screen, we used fluorescence anisotropy to measure the binding affinities of the identified inhibitors to the Dvl PDZ domain. In these experiments, each small molecule inhibitor was titrated into a solution of TMR-labeled PDZ domain, and the anisotropy change due to ligand binding was used to determine the binding affinity of the inhibitor to the Dvl PDZ domain. To verify the accuracy of this method, we first determined the binding affinity between the PDZ domain and the Dpr peptide; this value (11.0 Ϯ 1.4 M) was consistent with the value obtained by other methods (3,15). The fluorescence method yielded a measured binding affinity (K d ) value of 10.6 Ϯ 1.7 M between the PDZ domain and compound 3289-8625; the fluorescence data are plotted in Fig. 3. We also measured the binding of other, similar compounds identified in this study to the Dvl PDZ domain. As suggested by the NMR studies, all of these compounds had weaker binding affinities than did compound 3289-8625; for example, measured by the fluorescence anisotropy method, the binding affinity between compound 2372-2393 and the Dvl PDZ domain was 18.9 Ϯ 2.1 M.
The binding affinity between compound 3289-8625 and the Dvl PDZ domain is comparable with the binding affinity of the Dpr peptide to the Dvl PDZ domain. To further demonstrate that compound 3289-8625 can compete with the Dpr peptide, we also measured the binding affinity between the Dpr peptide and the PDZ domain by monitoring the change of fluorescence polarization during titration of unlabeled Dvl PDZ into a solution of the fluorescent ROX-labeled Dpr peptide. This method yielded a measured binding affinity of 6.1 Ϯ 0.4 M., which is in good agreement with the value obtained by using TMR-labeled PDZ domain as described above (the slight difference may due to the effects of fluorescence tag labeling or experimental errors). By using the same assay, we showed that compound 3289-8625 inhibited the interaction between the Dpr peptide and the PDZ domain in the manner of classical competitive inhibition (K I values, 4.9 Ϯ 1.7 M, Fig. 3B). The data clearly show that compound 3289-8625 and the Dpr peptide competed for the same site on the surface of Dvl PDZ domain. Indeed, the conformation of PDZ-bound compound 3289-8625 as calculated by FlexX lies in this same region and closely resembles the crystal structure conformation of the MTTV motif of the Dpr   peptide (10). When overlaid, the conformations of Dpr and this compound are seen to be very similar. The benzene ring on one end of the compound mimics the Val side chain of Dpr, and the benzene ring on the other end of the compound lies in the same region as the Met side chain. The C-terminal carboxylate ion of Dpr and that of the compound lie in the same region, where they can interact similarly with the PDZ domain (Fig. 4).
Inhibition of Wnt Signaling by Compound 3289-8625-The PDZ domain of Dvl interacts directly with the conserved sequence C-terminal to the seventh transmembrane helix of the Wnt receptor Fz (3), and Wnt signaling can be inhibited by blocking this interaction (3,15). To test whether compound 3289-8625 inhibits Wnt signaling in the cell, we evaluated the ability of compound 3289-8625 to block the canonical Wnt pathway using a 293 cell line stably transfected with a luciferase reporter under the control of lymphoid enhancer factor/T-cell factor binding sites (SuperTopflash) (34). To activate the Wnt pathway, we applied recombinant Wnt3a to the cells at either 10 ng/ml or 50 ng/ml (Fig. 5). When used at 3 M, 3289-8625 was effective at reducing luciferase activity by about 2-fold (Fig. 5A). However, at 50 ng/ml Wnt3a 3 M, compound 3289-8625 was not effective (Fig. 5B). These observations clearly indicated that the compound 3289-8625 competitively inhibits the Wnt signaling.
To further assess the ability of the compound to inhibit the Wnt signaling pathway in vivo, we used a previously established Xenopus system (15,19). Xenopus Wnt target gene siamois promoter-driven luciferase reporter construct (35) was injected alone or co-injected with synthetic Wnt3A mRNA (1 pg) or ␤-catenin RNA (100 pg) in the animal pole region at the twocell stage. Groups of 10 injected embryos were cultured in the absence or presence of different concentrations of the compound 3289-8625. Ectodermal explants were dissected at the late blastula stage and assayed for luciferase activity. Because the molecular mass of compound 3289-8625 is 374 daltons and its calculated octanol-water partition (logP) value and polar surface are 3.15 and 116 Å 2 , respectively, it is likely that the compound can penetrate the cell membrane. Therefore, instead of co-injecting the compound with Wnt3A mRNA as in earlier studies (15), we added the compounds to the culture medium at different concentrations. As shown in Fig. 6A, compound 3289-8625 was able to penetrate the keratin-rich outer embryonic cell membrane and to block Wnt signaling in a concentration-dependent fashion. At 10 M, it inhibited Wnt3A-induced luciferase activity by 2-fold. However, the compound had no effect if the Wnt signal was activated by ␤-catenin, indicating that the compound inhibited Wnt signaling upstream of ␤-catenin. Furthermore, ventro-vegetal injection of Wnt3A mRNA at the four-cell stage induced a complete secondary axis. When these Wnt3A-injected embryos were incubated with compound 3289-8625 at 25 M, the formation of a complete secondary axis was significantly reduced (Fig. 6, B and C). This further suggests that compound 3289-8625 was efficient at inhibiting Wnt signaling in vivo. Compound 3289-8625 Suppresses Wnt Signaling in the Hyaloid Vessel System-In the mouse eye, the canonical Wnt pathway is critical for scheduled vascular regression. In this system, resident macrophages produce Wnt7b and, through close contact with vascular endothelial cells of the temporary hyaloid vessels, Wnt7b activates the Wnt pathway to promote programmed cell death (36). To further examine the effects of the compound 3289-8625 on the Wnt signaling pathway in vivo, we took advantage of the earlier studies (36) and examined whether the compound was capable to interfere with the canonical Wnt pathway and reduce programmed cell death and vascular regression in the hyaloid vessel system. In the experiments, 120 nl of 12 M compound 3289-8625 was  injected into the vitreous of the eye, where the hyaloid vessels reside. Based on an estimated vitreous volume of 780 nl, we anticipated an immediate 6.5-fold dilution to give an initial concentration no less than 2.5 M. After 24 h, hyaloid vessel networks were removed by dissection, and a TUNEL analysis (37) was performed for detection of apoptotic cells (Fig. 7). Quantification revealed that compound 3289-8625 could suppress apoptosis consistent with inhibition of the Wnt pathway.
Compound 3289-8625 Slows the Growth of Prostate Cancer PC-3 Cells-Amplification, up-regulation, and overexpression of Dvl have been reported in many cancers, including prostate cancer (6 -8, 38, 39). Therefore, compounds that inhibit Wnt signaling by blocking the Fz-Dvl interaction may lead to the development of anticancer therapeutic agents (40). To test this premise, we experimentally assessed the effect of compound 3289-8625 on PC-3 prostate cancer cells, whose growth can be suppressed by inhibition of Wnt signaling (20).

DISCUSSION
Because of advances in technology, it is now possible to identify small molecules that can block the active site of a specific protein by screening an extensive small molecule data base in silico (41,42). However, the current technology is far from perfect. Frequently, a large percentage of the potential ligands identified by in silico screening yields false positive results, and additional methods are needed to evaluate these screening results. In our studies, NMR spectroscopy was used to perform such a final validation to identify small molecules that yield true positive results. This approach has proved to be very effective in our study. The aim of the current study was to identify and evaluate small molecule inhibitors of the PDZ domain of Dvl, and using the approach of NMR-assisted virtual screening, we identified several such inhibitors.
Wnt signaling plays an important role in embryonic development and the regulation of cell growth, and its inappropriate activation has been implicated in cancer and other human diseases such as ocular vascular disease (36,(43)(44)(45). For this reason, different elements within the Wnt signaling pathways have emerged as potentially useful targets for the development of therapeutic reagents (40). The cytoplasmic protein Dvl relays the Wnt signal from the membrane-bound Wnt receptor to downstream partners through the interaction between its PDZ domain and Wnt receptor Fz. Therefore, the ligand binding site on the surface of the PDZ domain provides a perfect site for formulating small molecule inhibitors to block the Wnt signal- ing pathway, and such small molecules may be a powerful tool to study Wnt signaling. In addition, it may also be helpful in formulating rational approaches to the development of novel pharmaceutical agents that can interfere with specific Wnt signal events that contribute to human diseases.
The best Dvl PDZ domain inhibitor identified in our studies, compound 3289-8625, has a binding affinity to Dvl, which is similar to that of Fz, a native Dvl PDZ domain binding partner in the Wnt signaling pathway (3). Therefore, this compound may be useful for blocking the Wnt signaling pathway at the Fz-Dvl interaction point. Indeed, different lines of experiments in our studies demonstrated that this compound does have the capability of blocking Wnt signaling in vivo.
For a successful small molecule inhibitor that can be used to target an intracellular protein-protein interaction event, the compound also has to be able to penetrate through the cell membrane. The molecular mass of compound 3289-8625 is 374 daltons. In addition, the calculated logP coefficient of the compound is 3.15, and the calculated polar surface area of the molecules is about 116.3 Å. All these parameters indicate that the compound is membrane-permeable, although we have not experimentally assessed its membrane permeability. Indeed, our studies with Xenopus embryos and cell culture showed that the compound was able to penetrate through the membranes of Xenopus embryos, 293 cells, and PC-3 cells. Furthermore, in vivo studies also showed that the compound could penetrate the membrane of vascular endothelial cells.
In summary, by screening data bases of small organic molecules that were readily synthesized, we were able to identify several PDZ binding compounds without the need for in-house design and synthesis of the ligand molecules. In addition, we also showed that one of the best inhibitors could effectively inhibit the Wnt signals in vivo. Therefore, the identified small molecule inhibitors not only can be used as templates for further chemical optimization(s) but also provide us a tool for dissection of the Wnt signaling pathways.  Cells were examined with the Promega MTT cell proliferation assay, and at the data points, cell proliferation was normalized against the control cells. Three sets of independent experiments were performed, and each data point represents the measurements from the three experiments. Error bars indicate S.E. B, the levels of ␤-catenin protein in the PC-3 cells treated with 80 M compound 3289-8625 for 72 h (close bar). The ␤-catenin proteins isolated from cell membrane and cytosolic fractions were analyzed in Western blots as described under "Experimental Procedures." The ␤-catenin protein levels were calibrated with integrin ␣3, a housekeeping protein that served as loading control, and normalized with the ␤-catenin levels in the cells that were treated with DMSO, i.e. 100% (open bar). The densitometries of protein bands presented here were quantified using the Scion Image software (Scion Corp., Frederick, MD). Each bar represents the mean of three Western blot results. Error bars indicate S.E.