Abstract
Wei Zhang1, Wenyan Lu2, Subramaniam Ananthan1, Mark J. Suto1 and Yonghe Li2
1Department of Chemistry, Drug Discovery Division, Southern Research Institute, Birmingham, Alabama 35205, United States
2Department of Oncology, Drug Discovery Division, Southern Research Institute, Birmingham, Alabama 35205, United States
Correspondence to:
Wei Zhang, email: [email protected]
Yonghe Li, email: [email protected]
Keywords: frizzled, Wnt/β-catenin, inhibitor, virtual screening, cancer
Received: March 29, 2017 Accepted: June 28, 2017 Published: September 06, 2017
ABSTRACT
Frizzled (Fzd) proteins are seven transmembrane receptors that belong to a novel and separated family of G-protein-coupled receptors (GPCRs). The Fzd receptors can respond to Wnt proteins to activate the canonical β-catenin pathway which is important for both initiation and progression of cancers. Disruption of the Wnt/β-catenin signal thus represents an opportunity for rational cancer prevention and therapy. Of the 10 members of the Fzd family, Fzd7 is the most important member involved in cancer development and progression. In the present studies, we applied structure-based virtual screening targeting the transmembrane domain (TMD) of Fzd7 to select compounds that could potentially bind to the Fzd7-TMD and block the Wnt/Fzd7 signaling and further evaluated them in biological assays. Six small molecule compounds were confirmed as Fzd7 inhibitors. The best hit, SRI37892, significantly blocked the Wnt/Fzd7 signaling with IC50 values in the sub-micromolar range and inhibited cancer cell proliferation with IC50 values around 2 μM. Our results provide the first proof of concept of targeting Fzd-TMD for the development of Wnt/Fzd modulators. The identified small molecular Fzd7 inhibitors can serve as a useful tool for studying the regulation mechanism(s) of Wnt/Fzd7 signaling as well as a starting point for the development of cancer therapeutic agents.
INTRODUCTION
Frizzled (Fzd) proteins are a family of seven transmembrane receptors (7-TM) for secreted Wnt glycolipoproteins. The Fzd receptors can respond to Wnt proteins in the presence of the Wnt co-receptor low density lipoprotein receptor-related protein 5 (LRP5) or LRP6 to activate the canonical β-catenin pathway. Activation of Wnt/β-catenin signaling has been found to be important for both initiation and progression of cancers of different tissues [1, 2]. Thus, disruption of the Wnt/β-catenin signal represents an opportunity for rational cancer chemoprevention and therapy [3, 4, 5].
Of the 10 members of the Fzd family, Fzd7 is the most important member regulating cancer development and progression. Fzd7 is frequently up-regulated in a variety of cancers including colon cancer [6, 7], breast cancer [8], hepatocellular carcinoma [9, 10], lymphoblastic leukemia [11], Wilm’s tumors [12, 13], glioblastoma [14], ovarian cancer [15] melanoma [16], chronic myeloid leukemia [17] and gastric cancer [18]. Recent studies have demonstrated that the anti-Fzd7 antibody [12, 19, 20, 21], siRNA/shRNA knockdown of Fzd7, [6, 7, 12, 15, 16, 17], soluble Fzd7 peptides [22] and small interfering peptides directed against Fzd7 [23] antagonized Wnt/β-catenin signaling in Fzd7-expressing cancer cells, and suppressed tumor progression and metastasis. Therefore, targeted inhibition of Fzd7 represents a rational and promising new approach for cancer therapy [24, 25].
The Fzd family (Fzd1-10) and the closely related Smoothened (SMO) have been listed by the International Union of Pharmacology (IUPHAR) as a novel and separated family of G-protein-coupled receptors (GPCRs) : the class-F (Frizzled) [26]. In addition to the traditional seven transmembrane helices domain (7TMD), Fzds also process a large N-terminal extracellular cystein rich domain (CRD) and a C-terminal intracellular PDZ (Psd-95/disc large/ZO-1 homologs)-binding domain. The CRD is fully exposed and interacts with extracellular Fzd binding proteins, including the receptor’s cognate ligands Wnt proteins [27, 28]. The PDZ-binding domain interact with the intracellular Fzd binding proteins, including Disheveled (Dvl), an essential and functionally necessary phosphoprotein for Fzd signaling pathways [29, 30]. While it has been demonstrated that interruption of the extracellular Wnt/CRD [19, 31], or the intracellular Fzd/Dvl [32, 33] interactions could effectively inhibit Wnt/β-catenin signaling, studies targeting the TMD of Fzd have not been reported to date.
The TMD of Fzd shares the same topology as classic GPCRs, including the conventional 7TM core, the intracellular C-terminal helix-8, the conserved cysteine residues in the extracellular loops as well as other common residues in the TMD of GPCRs [34, 35]. The genuine GPCR features of Fzd, including both G-protein and β-arrestin dependent signaling pathways, have also been observed in recent years [36, 37, 38, 39]. As major targets for drug development, GPCRs have been studied for many years and it is now well accepted that GPCR proteins exist in a conformational equilibrium between active and inactive biophysical states [40, 41, 42; 43], and ligands that shift the equilibrium toward the active/inactive conformational states of GPCRs have had applications in many different therapeutic areas. Therefore, it would be interesting to find out whether targeting Fzd-TMD could be an alternative strategy to modulate the Wnt/Fzd signaling.
The recently reported crystal structure of the TMD of human SMO receptor bound to a small-molecule antagonist [44] further supports the notion of blocking Wnt signaling using small molecule inhibitors. It also provides an excellent structural template for molecular modeling study of Fzd-TMD. In the present study, we constructed a homology model of the TMD of Fzd7 and applied structure-based virtual screening (SBVS) to select putative small molecular Fzd7 inhibitors. Our study led to the identification of an Fzd7 inhibitor, SRI37892, with potent activities against Wnt/β-catenin signaling and cancer cell viability.
RESULTS
Homology modeling and structural analysis of Fzd7-TMD
The SMO receptor has high sequence similarity to the frizzled receptors. The TMDs of Fzd7 and SMO share 28% identical and 47% homologous residues. Based on the SMO crystal structural and multi-sequence alignment of the GPCR class-F family (see Supplementary Tables 1 and 2 in the supporting information), we constructed a homology model for the Fzd7-TMD (Figure 1A). The SMO antagonist, LY2940680, which was co-crystallized in the SMO crystal structure was kept in the Fzd7-TMD model as a reference molecule to hold the putative ligand binding site. The top part of this putative binding site of the Fzd7 model is exposed and surrounded by the extracellular loops (ELs) as well as the N-terminal linker region (NLR), the rest part of the binding site is buried in the TMD and formed by residues from transmembrane helices (THs) I, III, VI and VII. This binding site mainly consists of hydrophobic residues which form two hydrophobic cores at the top and the bottom regions of the pocket (Figure 1B and 1C). The top hydrophobic core is formed by residues from NLR (L242, Y244, F245), EL-2 (L407, L415) and EL-3 (W499, W503, P525), the bottom core includes residues from TH-I (A2531.35, Y2571.39), TH-III (L3433.32), TH-VI (Y4896.52), TH-VII (F5307.38, F5347.42) as well as a residue Y412 from EL-2. Results of multiple sequence alignment of the class-Frizzled showed that the properties of most of the binding site residues are well conserved (see Supplementary Table 1 in the supporting information), especially residues that form the two hydrophobic cores, which suggests that this putative ligand binding site may exist commonly in the Frizzled family. These results encouraged us to target this binding site of Fzd7-TMD for small molecule inhibitors.
Figure 1: Structural presentation of the homology model of Fzd7-TMD. (A) Ribbon presentation of the Fzd-TMD model with LY2940680 (shown in solid sticks) buried in the putative binding site. (B) Close-up view of the binding site residues that form the top and the bottom hydrophobic cores. (C) 2D interaction diagram showing all binding site residues. Residues are colored by properties: hydrophobic (green), polar (cyan), positively charged (blue), negatively charged (red).
Identification of Fzd7 inhibitors through SBVS
To identify compounds that can bind to the Fzd7-TMD and block Wnt/Fzd7 signaling, we conducted structure based virtual screening (SBVS) targeting the TMD binding site of the constructed Fzd7 model. 500,000 structurally diverse compounds were screened through a three-step docking/scoring process using the Glide program. The docked results of 738 top-scored compounds at the final docking stage were further visually inspected to select potential Fzd7 inhibitors that showed both shape and electrostatic complementarities to the receptor binding site. We hypothesized that occupying the two hydrophobic cores could be important for ligand binding and gave certain priorities to compounds that can fit into both hydrophobic regions, especially to those that fitted well to the bottom hydrophobic core. By taking into consideration structural diversity in our final selection, 67 structurally representative compounds were purchased and tested for their inhibitory effect on Wnt/β-catenin signaling in Wnt3-expressing HEK293 cells. At 10 μM concentrations, 3 of 67 compounds, SRI35959, SRI35961 and SRI35958, significantly ( > 50%) blocked Wnt3A induced Wnt/β-catenin signaling and their inhibitory effects were further confirmed in concentration dependent assays with IC50 values of 2.9, 8.6 and 11.9 μM, respectively (Figures 2 and 3A). SRI35959, the most potent of these three hits, was further tested in LRP6-expressing HEK293 cells, and was found to display a similar inhibition potency with an IC50 value of 4.3 μM (Figure 3B).
Figure 2: Structures of the identified small molecule Fzd7 inhibitors. The IC50 values for inhibition of Wnt/β-catenin signaling were measured in Wnt3A-expressing HEK293 cells as described in Figure 3A and 3C.
Figure 3: SRI35959 and SRI37892 block Wnt/β-catenin signaling in HEK293 cells induced by Wnt3A and LRP6. HEK293 cells in 24-well plates were transiently transfected with the Wnt3A plasmid (A, C) or the LRP6 plasmid (B, D) along with the Super8XTOPFlash luciferase construct and the β-galactosidase-expressing vector in each well. After 24 h incubation, cells were treated with SRI35959 (A, B) or SRI37892 (C, D) at the indicated concentrations for 24 h. The luciferase activity was then measured with normalization to the activity of the β-galactosidase. Values are the average of triple determinations with the SD indicated by error bars. *P < 0.05, **P < 0.01 versus corresponding control value.
While the 67 selected compounds included structurally diverse compounds (see Supplementary Table 3 in the supporting information), the three hits identified clearly share some structural similarities, including the presence of a common phenylbenzimidazole unit. The most potent compound, SRI35959, however, possesses a potentially reactive α,β-unsaturated amide group and a 1,3-benzdioxole moiety that is prone for metabolic conversion to toxic metabolites [45]. Therefore, to identify additional compounds with improved potency and devoid of structural liabilities, we conducted analog searching and docking-screening to select a second set of compounds. By using different analog search methods, including similarity, substructure and topomer searches, we assembled 5000 “analogs” and docked them into the binding site of our Fzd7-TMD model. Following the same compound selection procedures described above, 35 analogs (Supplementary Table 4) were finally purchased from the top-scored results and tested in the Wnt/β-catenin assay. Three of 35 compounds, SRI37892, SRI37893 and SRI34284, were confirmed as actives with IC50 values of 0.66, 9.9 and 11.9 μM, respectively (Figures 2 and 3C). The best hit, SRI37892, also displayed potent activity against Wnt/β-catenin signaling in LRP6-expressing HEK293 cells with an IC50 value of 0.78 μM (Figure 3D).
All the 6 hits from the two sets of tested compounds share the same phenylbenzimidazole fragment, suggesting this group may be important for receptor binding. 4 of the 6 hits (SRI35959, SRI35961, SRI37893 and SRI34284) possess a potentially metabolically oxidizable p-phenylenediamine-like unit. The most potent compound, SR37892, however, does not contain unfavorable chemical features and was therefore selected for further studies described below.
Inhibitory effects of SRI37892 on Wnt/β-catenin signaling in cancer cells
Aberrant Wnt/β-catenin signaling is associated with a poorer prognosis in breast cancer patients [46], and is predominantly found in triple negative breast cancer (TNBC) which is distinguished by negative immunohistochemical assays for expression of the estrogen and progesterone receptors (ER/PR) and human epidermal growth factor receptor-2 (HER2) [47, 48, 49]. It has been reported that FZD7 is upregulated in TNBC, and that FZD7 plays an important role on Wnt/β-catenin signaling in TNBC cells and cancer cell proliferation [50]. Therefore, we tested SRI37892 in TNBC HS578T and BT549 cells to confirm its inhibitory effect on Wnt/β-catenin signaling. As expected, treatment of SRI37892 at 1 or 2 μM resulted in suppression of LRP6 phosphorylation, down-regulation of cytosolic free β-catenin level, and inhibition of expression of specific Wnt targets axin2 and survivin in HS578T and BT549 cells (Figure 4). Moreover, SRI37892 displayed potent activity against HS578T and BT549 cell proliferation with IC50 values of 2.2 and 1.9 μM, respectively (Figure 5A). Finally, SRI37892 at 0.5-2 μM significantly suppressed colony formation in HS578T and BT549 cells (Figure 5B).
Figure 4: Effects of SRI37892 on Wnt/β-catenin signaling in breast cancer HS578T and BT549 cells. Cancer cells in 6-well plates were treated with SRI37892 at the indicated concentrations for 24 h. The levels of cytosolic free β-catenin, total cellular β-catenin, Fzd7, LRP6, phospho-LRP6, axin2 and survivin were examined by Western blotting. All the samples were also probed with anti-human actin antibody to verify equal loading.
Figure 5: SRI37892 inhibits breast cancer cell viability and colony formation. (A) Breast cancer HS578T and BT549 cells in 96-well plates were treated with SRI37892 at the indicated concentrations for 96 h. Cell viability was measured by the CellTiter-Glo assay. (B) Breast cancer HS578T and BT549 cells were treated with SR37892 at the indicated concentrations for 10–12 days. The medium was changed every 3 days. Colonies were fixed with formaldehyde and stained with Crystal Violet. All values are averages of triplicate determinations with the S.D. indicated by error bars. **P <0.01 compared with corresponding control value.
Predicted binding mode of the active compounds
To gain structural insights into Fzd7-inhibitor binding, we further conducted docking studies following a more sophisticated induced-fit docking (IFD) protocol which allows binding site residues to be flexible. The two most potent hits, SRI35959 and SRI37892, were docked separately into the putative binding site of the Fzd7-TMD model. The docked models suggested that SRI37892 bound relatively tighter than SRI35959 with a docking score (which mimics the binding affinity) of -12.0 kcal/mol as compared to -10.3 kcal/mol of SRI35959, which is consistent with the observed experimental results that SRI37892 is a more potent inhibitor. Both compounds fit well into the Fzd7 binding site. In the docked models (Figure 6), the phenylbenzimidazole moiety of these two compounds docked very similarly: both occupied the lower hydrophobic pocket of the binding site and formed multiple Pi-Pi interactions with the aromatic sidechains of the surrounding residues (Trp534, Phe530 and His303). While the other end of the two compounds occupied the top hydrophobic core of the binding site, there are clear differences among the protein-inhibitor interactions. SRI35959 has a relatively linear conformation and thus its benzodioxole end can only occupy part of the top hydrophobic core. In contrast, the kinked conformation of SRI37892 enables its dihydroquinolinone end to interact with residues at the other side of the top hydrophobic core, including Leu242 and Met243 from the N-terminal linker region. While the middle fragment of SRI35959 did not form any close interactions with the receptor, the amide group of SRI37982 formed a hydrogen bond with Glu492 and its middle thiazole ring fitted tightly into a hydrophobic region formed by Pro525 and Trp499. Taken together, occupying both the top and the bottom hydrophobic regions of the binding site is likely important for ligand binding, and forming close hydrophobic interactions and additional hydrogen bond(s) with nearby residues can further contribute to favorable Fzd7-inhibitor interactions.
Figure 6: Predicted binding modes of identified Fzd7 inhibitors. The docked models of SRI35953 (A) and SRI37892 (B). Inhibitor molecules are represented in green-colored solid sticks. Protein residues are shown in gray lines. 2D ligand interaction diagram of SRI35953 (C) and SRI37892 (D). Pi-Pi stacking is marked with solid green lines. Hydrogen bond is marked with dashed purple lines. Residues are colored by properties: hydrophobic (green), polar (cyan), positively charged (blue), negatively charged (red).
DISCUSSION
In addition to the canonical β-catenin dependent signaling, there are also other noncannonical Wnt/Fzd signaling pathways that are independent of β-catenin, including Wnt/Ca2+ and Wnt/PCP signaling [51]. The mechanism(s) that controls the activation of a Wnt toward a certain signaling path are still unclear. Meanwhile, while structure biology studies have demonstrated that Wnts bind directly to the extracellular CRD of Fzd [52], other studies suggested that the CRD might be dispensable for the Wnt/Fzd signaling [53]. A possible explanation for such a conundrum is that CRD may act as a ligand fishing rod that brings Wnt close to the receptor thus facilitating further interactions with other regions of the core TMD to activate the Wnt/Fzd signaling [54]. Nonetheless, the underlying mechanism(s) of Fzds modulated signaling is still unclear partly because of the lack of proper research tools. To the best of our knowledge, this is the first report on small molecule Wnt inhibitors that target the Fzd-TMD. The identified Fzd7 inhibitors thus provide a useful tool for the study of Wnt/Fzd signaling mechanism.
Wnt/β-catenin signaling is activated aberrantly in many human cancers, and FZD7 is frequently up-regulated in a variety of cancers including breast and colorectal cancer [24, 25]. While genetic mutations of certain components of the Wnt/β-catenin pathway are not the predominate mechanism associated with many types of cancer including breast cancer [4], mutations of APC and CTNNB1 (β-catenin encoding gene) are two major factors for Wnt/β-catenin signaling overactivation in colorectal cancers [55, 56]. However recent studies also indicate that the levels of Wnt/β-catenin signaling in colorectal cancer cells can be modulated on the cell surface [57, 58, 59]. Fzd7 can activate Wnt/β-catenin signaling in colorectal cancer cells despite the presence of APC or CTNNB1 mutations [6, 7]. In the present study, we demonstrated that SRI37892 has a potent activity against Wnt/β-catenin signaling in TNBC cells. It would be interesting to see whether SRI37892 can indeed effectively inhibit Wnt/β-catenin signaling in colorectal cancer cells harboring APC and CTNNB1 mutations.
Wnt/β-catenin signaling is a critical regulator of stem cells [2]. Recent studies have demonstrated that Fzd7 is required for maintenance of the pluripotent state in human embryonic stem cells [60], and functions as a Wnt receptor for robust Wnt-dependent processes in Lgr5(+) intestinal stem cells [61]. Moreover, Fzd7 plays an important role in activation of the non-canonical Wnt signaling pathway during human embryonic hematopoiesis [62]. Further studies are required to examine whether the identified compounds inhibit non-canonical Wnt signaling and modulate stem cell fate and function.
While increased Wnts expression has been reported in different forms of cancers, reduced expression of Wnts is associated with other disorders such as obesity [63]. Recent studies also demonstrated that the Wnt/β-catenin pathway is a critical regulator of osteoclastogenesis, and that Wnt activation attenuates osteoclast differentiation and bone resorption, which suggests that activating Wnt/β-catenin signaling is an attractive strategy for the treatment of bone loss disorders [64, 65, 66]. Targeting the TMD of GPCRs is a clinically validated strategy for drug development. Ligands that bind to the TMD and modulate the conformations of GPCRs have been successfully applied to develop different drug candidates, including agonists, antagonist as well as biased modulators that regulate specific signaling pathways [67, 68, 69, 70, 71]. The GPCR characterizations of Fzds have been established in recent years. Therefore, targeting Fzd-TMD could be an attractive therapeutic strategy to modulate Wnt/Fzd signaling. Nonetheless, so far there were no studies demonstrating the feasibility of such a strategy. In the present study, by specifically targeting a putative binding site at the Fzd7-TMD, we successfully identified small molecule compounds that significantly blocked the Wnt/Fzd signaling and inhibited the proliferation of cancer cells, thus providing a strong support to the concept of targeting Fzd-TMD for the development of chemotherapeutic candidates. Sequence analysis shows that most of the SRI37892 binding site residues of Fzd7-TMD are well conserved (see Supplementary Tables 1 and 2 in the supporting information). Therefore, the compounds identified here may also bind to other FZD receptors. Nevertheless, the identified SRI37892, a potent Fzd7 inhibitor, represents a potential lead for further drug development.
MATERIALS AND METHODS
Molecular modeling
Structural model generation and molecular docking studies were conducted using the programs of the Schrödinger Suite 2015 (Schrödinger, LLC, New York, NY, 2015).
Homology modeling was conducted using the Prime program. The crystal structure of SMO-LY2940680 (a small-molecule SMO antagonist) complex (PDB ID: 4JVK) was used as the structural template. The sequences of SMO and frizzled receptors (Fzds1-10) was first aligned using ClusterX. The alignment results were then input into Prime to build the structural model of Fzd7-TMD. We also used the Ballesteros and Weinstein numbering system [72] on the basis of structural superposition. Specifically, for each of the THs, the following residues are assigned number 50: T2681.50, Y2962.50, W3613.50, W3864.50, Y4325.50, A4876.50, I5427.50.
Molecular docking studies were performed using the Glide program with the default parameters. The 3D structures of ligands were prepared using the LigPrep program. The induced-Fit-Docking (IFD) protocol [73], which is capable of sampling dramatic side-chain conformational changes as well as minor changes in protein backbone structure, was applied to explore the potential binding modes of identified active compounds. Residues within 5 Å of the docked ligands were allowed to be flexible and the docked results were finally scored using the extra-precision (XP) mode of Glide.
Structure based virtual screening
For the purpose of applying structure-based virtual screening, we assembled a virtual library consisting structurally diverse compounds selected from different commercial sources. Specifically, structures of a total of approximately eight million compounds were downloaded directly from the websites of twelve vendors (Asinex, Chembridge, ChemDiv, Enamine, FCH group, InterBioScreen, Life chemicals, Maybridge, Princeton, TIMTEC, SPECS and Vitas-M). From the eight million commercially available compounds, we first identified the most diverse set of 100,000 structurally representative compounds using the clustering and diversity analysis protocols of Pipeline Pilot; [74] then for each of the 100,000 compounds, we further selected four structurally most similar analogs from the rest of the eight million compounds based on the Tanimoto coefficients calculated from the 2D structural fingerprints. Such an assembled library of 500,000 compounds covers a large portion of chemical space, contains diverse structures, and could provide initial structure-activity relationship (SAR) information.
Structure-based virtual screening was performed using the virtual screening workflow of maestro, which adopts a three-step docking/scoring protocol implemented in Glide. 3D conformations of the 500,000 assembled compounds were prepared using LigPrep, and a total of 923,408 conformers were generated. All of the 923,408 conformers were first docked into the ligand binding site of the Fzd7-TMD model using the high-throughput virtual screening (HTVS) mode; the top 5% best-scored conformers were then re-docked and scored using the standard-precision (SP) mode; the top 10% of the SP docked conformers were further docked and scored using the extra-precision (XP) mode. Finally, the top 25% of the best XP-scored conformers were output for visual examination.
Analog search
Commercial analogs of confirmed active compounds were identified through substructure search and similarity search protocols of Pipeline Pilot and the TOPOMER search module of SYBYLX program. The identified analogs were then docked to the binding site of the Fzd7-TMD model. The top-scored docked models were further visually examined for structural complementary in both shape and polarity to the receptor to select a final list of compounds for biological evaluations.
Cell culture
All cell lines were obtained from ATCC and grown under standard cell culture conditions at 37°C in a humidified atmosphere with 5% CO2. The cells were cultured in DMEM medium containing 10% of FBS, 2 mM of L-glutamine, 100 units/ml of penicillin, and 100 μg/ml of streptomycin.
Luciferase Wnt reporter assay
Plasmid pCS-Myc-hLRP6 containing the full-length human LRP6 cDNA was provided by Dr. Christof Niehrs (Deutsches Krebsforschungszentrum, Heidelberg, Germany), and the Super8XTOPFlash luciferase construct was provided by Dr. Randall T. Moon (University of Washington, Seattle). Plasmid pcDNA3-Wnt3-HA was constructed as described previously [75]. HEK293 cells were plated into 24-well plates. After overnight culture, the cells were transiently transfected with the Super8XTOPFlash luciferase construct and β-galactosidase-expressing vector (Promega) along with or without Wnt3A or LRP6 plasmid. After 24 h incubation, cells were treated with the test compounds at the indicated concentrations. Cells were then lysed 24 h later and both luciferase and β-galactosidase activities were determined. The luciferase activity was normalized to the β-galactosidase activity.
Western blotting
Polyclonal anti-LRP6 and monoclonal anti-survivin were from Santa Cruz Biotechnology. Polyclonal anti-Fzd7 was from Novus Biologicals. Monoclonal anti-phospho-LRP6 and anti-axin2 were purchased from Cell Signaling Technology. Monoclonal anti-β-catenin was from BD Biosciences. Monoclonal anti-actin was from Sigma. Cells in 6-well plates were lysed in 0.5 ml of lysis buffer (phosphate-buffered saline containing 1% Triton X-100 and 1 mM phenylmethylsulfonyl fluoride) at 4°C for 30 min. Equal quantities of protein were subjected to SDS-PAGE under reducing conditions. Following transfer to immobilon-P transfer membrane, successive incubations with a primary antibody, and a horseradish peroxidase-conjugated secondary antibody were carried out for 60-120 min at room temperature. The immunoreactive proteins were then detected using the ECL system. Films showing immunoreactive bands were scanned by HP Scanjet 5590 (Hewlett Packard, Palo Alto, CA).
Cytosolic free β-catenin analysis
Plasmid pGST-E-cadherin was provided by Dr. Gail Johnson (University of Rochester). The level of cytosolic free β-catenin was analyzed by the GST-E-cadherin binding assay as previously described [75]. Uncomplexed cytosolic free β-catenin present in 100 μg of total cell lysate was subjected to SDS-PAGE and detected using the monoclonal antibody to β-catenin (BD Biosciences).
Cell viability assay
Cells were seeded into 96-well tissue culture treated microtiter plates at a density of 4000 cells/well. After overnight incubation, the cells were treated with each individual compound for 96 h. Cell viability was measured by the CellTiter-Glo Assay (Promega).
Colony formation assay
Cancer cells were seeded at a density of 800 cells/well into six-well plates. Sixteen hours after the plates had been set up, SR37892 was added, and media were replenished every 3 days. After being incubated for 12 days, colonies were fixed with 4% formaldehyde, stained with 0.5 mg/ml crystal violet, and imaged on a FluorChem HD2 Imager System (Alpha Innotech).
Author contributions
W.Z. coordinated the study, performed the computational studies and wrote the paper. W.L. and Y.L. were responsible for the design and performance of biological evaluations and the related writings. S.A. and M.S. contributed their medicinal chemistry expertise to the study and the manuscript.
CONFLICTS OF INTEREST
The authors declare no competing financial interest.
FUNDING
This work was partially supported by Southern Research Internal Fund and the National Institutes of Health grant R21CA182056.
REFERENCES
1. Niehrs C. The complex world of WNT receptor signalling. Nat Rev Mol Cell Biol. 2012; 13:767-779.
2. Clevers H, Nusse R. Wnt/beta-catenin signaling and disease. Cell. 2012; 149:1192-1205.
3. Anastas JN, Moon RT. WNT signalling pathways as therapeutic targets in cancer. Nat Rev Cancer. 2013; 13:11-26.
4. King TD, Suto MJ, Li Y. The Wnt/beta-catenin signaling pathway: a potential therapeutic target in the treatment of triple negative breast cancer. J Cell Biochem. 2012; 113:13-18.
5. Polakis P. Drugging Wnt signalling in cancer. EMBO J. 2012; 31:2737-2746.
6. Ueno K, Hazama S, Mitomori S, Nishioka M, Suehiro Y, Hirata H, Oka M, Imai K, Dahiya R, Hinoda Y. Down-regulation of frizzled-7 expression decreases survival, invasion and metastatic capabilities of colon cancer cells. Br J Cancer. 2009; 101:1374-1381.
7. Ueno K, Hiura M, Suehiro Y, Hazama S, Hirata H, Oka M, Imai K, Dahiya R, Hinoda Y. Frizzled-7 as a potential therapeutic target in colorectal cancer. Neoplasia. 2008; 10:697-705.
8. Chakrabarti R, Wei Y, Hwang J, Hang X, Andres Blanco M, Choudhury A, Tiede B, Romano RA, DeCoste C, Mercatali L, Ibrahim T, Amadori D, Kannan N, et al. ΔNp63 promotes stem cell activity in mammary gland development and basal-like breast cancer by enhancing Fzd7 expression and Wnt signaling. Nat Cell Biol. 2014; 16:1004-15.
9. Bengochea A, de Souza MM, Lefrancois L, Le Roux E, Galy O, Chemin I, Kim M, Wands JR, Trepo C, Hainaut P, Scoazec JY, Vitvitski L, Merle P. Common dysregulation of Wnt/Frizzled receptor elements in human hepatocellular carcinoma. Br J Cancer. 2008; 99:143-150.
10. Merle P, Kim M, Herrmann M, Gupte A, Lefrancois L, Califano S, Trepo C, Tanaka S, Vitvitski L, de la Monte S, Wands JR. Oncogenic role of the frizzled-7/beta-catenin pathway in hepatocellular carcinoma. J Hepatol. 2005; 43:854-862.
11. Khan NI, Bradstock KF, Bendall LJ. Activation of Wnt/beta-catenin pathway mediates growth and survival in B-cell progenitor acute lymphoblastic leukaemia. Br J Haematol. 2007; 138:338-348.
12. Pode-Shakked N, Harari-Steinberg O, Haberman-Ziv Y, Rom-Gross E, Bahar S, Omer D, Metsuyanim S, Buzhor E, Jacob-Hirsch J, Goldstein RS, Mark-Danieli M, Dekel B. Resistance or sensitivity of Wilms’ tumor to anti-FZD7 antibody highlights the Wnt pathway as a possible therapeutic target. Oncogene. 2011; 30:1664-1680.
13. Li CM, Guo M, Borczuk A, Powell CA, Wei M, Thaker HM, Friedman R, Klein U, Tycko B. Gene expression in Wilms’ tumor mimics the earliest committed stage in the metanephric mesenchymal-epithelial transition. Am J Pathol. 2002; 160:2181-2190.
14. Schiffgens S, Wilkens L, Brandes AA, Meier T, Franceschi E, Ermani M, Hartmann C, Sandalcioglu IE, Dumitru CA. Sex-specific clinicopathological significance of novel (Frizzled-7) and established (MGMT, IDH1) biomarkers in glioblastoma. Oncotarget. 2016; 7:55169-55180. https://doi.org/10.18632/oncotarget.10465.
15. Asad M, Wong MK, Tan TZ, Choolani M, Low J, Mori S, Virshup D, Thiery JP, Huang RY. FZD7 drives in vitro aggressiveness in Stem-A subtype of ovarian cancer via regulation of non-canonical Wnt/PCP pathway. Cell Death Dis. 2014; 5:e1346.
16. Tiwary S, Xu L. FRIZZLED7 is required for tumor inititation and metastatic growth of melanoma cells. PLoS One. 2016; 11:e0147638.
17. Liu N, Zang S, Liu Y, Wang Y, Li W, Liu Q, Ji M, Ma D, Ji C. FZD7 regulates BMSCs-mediated protection of CML cells. Oncotarget. 2016; 7:6175-6187. https://doi.org/10.18632/oncotarget.6742.
18. Schmuck R, Warneke V, Behrens HM, Simon E, Weichert W, Rocken C. Genotypic and phenotypic characterization of side population of gastric cancer cell lines. Am J Pathol. 2011; 178:1792-1804.
19. Gurney A, Axelrod F, Bond CJ, Cain J, Chartier C, Donigan L, Fischer M, Chaudhari A, Ji M, Kapoun AM, Lam A, Lazetic S, Ma S, et al. Wnt pathway inhibition via the targeting of Frizzled receptors results in decreased growth and tumorigenicity of human tumors. Proc Natl Acad Sci U S A. 2012; 109:11717-11722.
20. Riley RS, Day ES. Frizzled7 antibody-functionalized nanoshells enable multivalent binding for Wnt signaling inhibition in triple negative breast cancer cells. Small. 2017; 13.
21. Nickho H, Younesi V, Aghebati-Maleki L, Motallebnezhad M, Majidi Zolbanin J, Movassagh Pour A, Yousefi M. Developing and characterization of single chain variable fragment (scFv) antibody against frizzled 7 (Fzd7) receptor. Bioengineered. 2016; 1-10.
22. Wei W, Chua MS, Grepper S, So SK. Soluble Frizzled-7 receptor inhibits Wnt signaling and sensitizes hepatocellular carcinoma cells towards doxorubicin. Mol Cancer. 2011; 10:16.
23. Nambotin SB, Lefrancois L, Sainsily X, Berthillon P, Kim M, Wands JR, Chevallier M, Jalinot P, Scoazec JY, Trepo C, Zoulim F, Merle P. Pharmacological inhibition of Frizzled-7 displays anti-tumor properties in hepatocellular carcinoma. J Hepatol. 2011; 54:288-299.
24. King TD, Zhang W, Suto MJ, Li Y. Frizzled7 as an emerging target for cancer therapy. Cell Signal. 2012; 24: 846-851.
25. Phesse T, Flanagan D, Vincan E. Frizzled7: a promising Achilles’ heel for targeting the Wnt receptor complex to treat cancer. Cancers (Basel). 2016; 8.
26. Foord SM, Bonner TI, Neubig RR, Rosser EM, Pin JP, Davenport AP, Spedding M, Harmar AJ. International Union of Pharmacology. XLVI. G protein-coupled receptor list. Pharmacol Rev. 2005; 57:279-288.
27. Nusse R. Wnts and Hedgehogs: lipid-modified proteins and similarities in signaling mechanisms at the cell surface. Development. 2003; 130:5297-5305.
28. Schulte G, Bryja V. The Frizzled family of unconventional G-protein-coupled receptors. Trends Pharmacol Sci. 2007; 28:518-525.
29. Cong F, Schweizer L, Varmus H. Wnt signals across the plasma membrane to activate the beta-catenin pathway by forming oligomers containing its receptors, Frizzled and LRP. Development. 2004; 131:5103-5115.
30. Wallingford JB, Habas R. The developmental biology of Dishevelled: an enigmatic protein governing cell fate and cell polarity. Development. 2005; 132:4421-4436.
31. Lee HJ, Bao J, Miller A, Zhang C, Wu J, Baday YC, Guibao C, Li L, Wu D, Zheng JJ. Structure-based discovery of novel small molecule Wnt signaling inhibitors by targeting the cysteine-rich domain of frizzled. J Biol Chem. 2015; 290:30596-30606.
32. Shan JF, Shi DL, Wang JM, Zheng J. Identification of a specific inhibitor of the dishevelled PDZ domain. Biochemistry. 2005; 44:15495-15503.
33. Grandy D, Shan J, Zhang X, Rao S, Akunuru S, Li H, Zhang Y, Alpatov I, Zhang XA, Lang RA, Shi DL, Zheng JJ. Discovery and characterization of a small molecule inhibitor of the PDZ domain of dishevelled. J Biol Chem. 2009; 284:16256-16263.
34. Schulte G. International union of basic and clinical pharmacology. LXXX. The class frizzled receptors. Pharmacol Rev. 2010; 62:632-967.
35. Lagerstrom MC, Schioth HB. Structural diversity of G protein-coupled receptors and significance for drug discovery. Nat Rev Drug Discov. 2008; 7:339-357.
36. Angers S, Moon RT. Proximal events in Wnt signal transduction. Nat Rev Mol Cell Biol. 2009; 10:468-477.
37. Egger-Adam D, Katanaev VL. Trimeric G protein-dependent signaling by Frizzled receptors in animal development. Front Biosci. 2008; 13:4740-4755.
38. Koval A, Katanaev VL. Wnt3a stimulation elicits G-protein-coupled receptor properties of mammalian Frizzled proteins. Biochem J. 2011; 433L:435-440.
39. Schulte G, Schambony A, Bryja V. beta-Arrestins - scaffolds and signalling elements essential for WNT/Frizzled signalling pathways? Br J Pharmacol. 2010; 159:1051-1058.
40. Yeagle PL, Albert AD. A conformational trigger for activation of a G protein by a G protein-coupled receptor. Biochemistry. 2003; 42:1365-1368.
41. Kobilka BK, Deupi X. Conformational complexity of G-protein-coupled receptors. Trends Pharmacol Sci. 2007; 28:397-406.
42. Hamm HE. How activated receptors couple to G proteins. Proc Natl Acad Sci U S A. 2001; 98:4819-4821.
43. Fanelli F, De Benedetti PG, Raimondi F, Seeber M. Computational modeling of intramolecular and intermolecular communication in GPCRs. Curr Protein Pept Sci. 2009; 10:173-185.
44. Wang C, Wu H, Katritch V, Han GW, Huang XP, Liu W, Siu FY, Roth BL, Cherezov V, Stevens RC. Structure of the human smoothened receptor bound to an antitumour agent. Nature. 2013; 497:338-343.
45. Stepan AF, Walker DP, Bauman J, Price DA, Baillie TA, Kalgutkar AS, Aleo MD. Structural alert/reactive metabolite concept as applied in medicinal chemistry to mitigate the risk of idiosyncratic drug toxicity: a perspective based on the critical examination of trends in the top 200 drugs marketed in the United States. Chem Res Toxicol. 2011; 24:1345-1410.
46. Lin SY, Xia W, Wang JC, Kwong KY, Spohn B, Wen Y, Pestell RG, Hung MC. Beta-catenin, a novel prognostic marker for breast cancer: its roles in cyclin D1 expression and cancer progression. Proc Natl Acad Sci U S A. 2000; 97:4262-4266.
47. Khramtsov AI, Khramtsova GF, Tretiakova M, Huo D, Olopade OI, Goss KH. Wnt/beta-catenin pathway activation is enriched in basal-like breast cancers and predicts poor outcome. Am J Pathol. 2010; 176:2911-2920.
48. Geyer FC, Lacroix-Triki M, Savage K, Arnedos M, Lambros MB, MacKay A, Natrajan R, Reis-Filho JS. beta-Catenin pathway activation in breast cancer is associated with triple-negative phenotype but not with CTNNB1 mutation. Mod Pathol. 2011; 24: 209-231.
49. Dey N, Barwick BG, Moreno CS, Ordanic-Kodani M, Chen Z, Oprea-Ilies G, Tang W, Catzavelos C, Kerstann KF, Sledge GW Jr, Abramovitz M, Bouzyk M, De P, Leyland-Jones BR. Wnt signaling in triple negative breast cancer is associated with metastasis. BMC Cancer. 2013; 13:537.
50. Yang L, Wu X, Wang Y, Zhang K, Wu J, Yuan YC, Deng X, Chen L, Kim CC, Lau S, Somlo G, Yen Y. FZD7 has a critical role in cell proliferation in triple negative breast cancer. Oncogene. 2011; 30:4437-4446.
51. Schulte G. International union of basic and clinical pharmacology. LXXX. The class Frizzled receptors. Pharmacol Rev. 2010; 62:632-667.
52. Janda CY, Waghray D, Levin AM, Thomas C, Garcia KC. Structural basis of Wnt recognition by Frizzled. Science. 2012; 337:59-64.
53. Chen CM, Strapps W, Tomlinson A, Struhl G. Evidence that the cysteine-rich domain of Drosophila Frizzled family receptors is dispensable for transducing Wingless. Proc Natl Acad Sci U S A. 2004; 101:15961-15966.
54. Povelones M, Nusse R. The role of the cysteine-rich domain of Frizzled in Wingless-Armadillo signaling. EMBO J. 2005; 24:3493-3503.
55. Kinzler KW, Vogelstein B. Lessons from hereditary colorectal cancer. Cell. 1996; 87:159-170.
56. Morin PJ, Sparks AB, Korinek V, Barker N, Clevers H, Vogelstein B, Kinzler KW. Activation of beta-catenin-Tcf signaling in colon cancer by mutations in beta-catenin or APC. Science. 1997; 275:1787-1790.
57. Voloshanenko O, Erdmann G, Dubash TD, Augustin I, Metzig M, Moffa G, Hundsrucker C, Kerr G, Sandmann T, Anchang B, Demir K, Boehm C, Leible S, et al. Wnt secretion is required to maintain high levels of Wnt activity in colon cancer cells. Nat Commun. 2013; 4:2610.
58. Giannakis M, Hodis E, Jasmine Mu X, Yamauchi M, Rosenbluh J, Cibulskis K, Saksena G, Lawrence MS, Qian ZR, Nishihara R, Van Allen EM, Hahn WC, Gabriel SB, et al. RNF43 is frequently mutated in colorectal and endometrial cancers. Nat Genet. 2014; 46:1264-1266.
59. Bond CE, McKeone DM, Kalimutho M, Bettington ML, Pearson SA, Dumenil TD, Wockner LF, Burge M, Leggett BA, Whitehall VL. RNF43 and ZNRF3 are commonly altered in serrated pathway colorectal tumorigenesis. Oncotarget. 2016; 7:70589-70600. https://doi.org/10.18632/oncotarget.12130.
60. Fernandez A, Huggins IJ, Perna L, Brafman D, Lu D, Yao S, Gaasterland T, Carson DA, Willert K. The WNT receptor FZD7 is required for maintenance of the pluripotent state in human embryonic stem cells. Proc Natl Acad Sci U S A. 2014; 111:1409-1414.
61. Flanagan DJ, Phesse TJ, Barker N, Schwab RH, Amin N, Malaterre J, Stange DE, Nowell CJ, Currie SA, Saw JT, Beuchert E, Ramsay RG, Sansom OJ, et al. Frizzled7 functions as a Wnt receptor in intestinal epithelial Lgr5(+) stem cells. Stem Cell Reports. 2015; 4:759-767.
62. Vijayaragavan K, Szabo E, Bosse M, Ramos-Mejia V, Moon RT, Bhatia M. Noncanonical Wnt signaling orchestrates early developmental events toward hematopoietic cell fate from human embryonic stem cells. Cell Stem Cell. 2009; 4:248-262.
63. Christodoulides C, Scarda A, Granzotto M, Milan G, Dalla Nora E, Keogh J, De Pergola G, Stirling H, Pannacciulli N, Sethi JK, Federspil G, Vidal-Puig A, Farooqi IS, et al. WNT10B mutations in human obesity. Diabetologia. 2006; 49:678-684.
64. Albers J, Keller J, Baranowsky A, Beil FT, Catala-Lehnen P, Schulze J, Amling M, Schinke T. Canonical Wnt signaling inhibits osteoclastogenesis independent of osteoprotegerin. J Cell Biol. 2013; 200:537-549.
65. Otero K, Shinohara M, Zhao H, Cella M, Gilfillan S, Colucci A, Faccio R, Ross FP, Teitelbaum SL, Takayanagi H, Colonna M. TREM2 and beta-catenin regulate bone homeostasis by controlling the rate of osteoclastogenesis. J Immunol. 2012; 188:2612-2621.
66. Wei W, Zeve D, Suh JM, Wang X, Du Y, Zerwekh JE, Dechow PC, Graff JM, Wan Y. Biphasic and dosage-dependent regulation of osteoclastogenesis by beta-catenin. Mol Cell Biol. 2011; 31:4706-4719.
67. Tikhonova IG, Sum CS, Neumann S, Engel S, Raaka BM, Costanzi S, Gershengorn MC. Discovery of novel agonists and antagonists of the free fatty acid receptor 1 (FFAR1) using virtual screening. J Med Chem. 2008; 51:625-633.
68. Engel S, Skoumbourdis AP, Childress J, Neumann S, Deschamps JR, Thomas CJ, Colson AO, Costanzi S, Gershengorn MC. A virtual screen for diverse ligands: discovery of selective G protein-coupled receptor antagonists. J Am Chem Soc. 2008; 130:5115-5123.
69. Pupo AS, Duarte DA, Lima V, Teixeira LB, Parreiras ES, Costa-Neto CM. Recent updates on GPCR biased agonism. Pharmacol Res. 2016; 112:49-57.
70. Moller D, Gmeiner P. Arrestin-bound rhodopsin: a molecular structure and its impact on the development of biased GPCR ligands. Angew Chem Int Ed Engl. 2015; 54:13166-13168.
71. Chen JZ, Wang J, Xie XQ. GPCR structure-based virtual screening approach for CB2 antagonist search. J Chem Inf Model. 2007; 47:1626-1637.
72. Ballesteros JA, Weinstein H. Integrated methods for the construction of three-dimensional models and computational probing of structure-function relations in G protein-coupled receptors. Methods Neurosci. 1995; 25:366-428.
73. Sherman W, Day T, Jacobson MP, Friesner RA, Farid R. Novel procedure for modeling ligand/receptor induced fit effects. J Med Chem. 2006; 49:534-553.
74. Baurin N, Baker R, Richardson C, Chen I, Foloppe N, Potter A, Jordan A, Roughley S, Parratt M, Greaney P, Morley D, Hubbard RE. Drug-like annotation and duplicate analysis of a 23-supplier chemical database totalling 2.7 million compounds. J Chem Inf Comput Sci. 2004; 44:643-651.
75. Lin C, Lu W, Zhai L, Bethea T, Berry K, Qu Z, Waud WR, Li Y. Mesd is a general inhibitor of different Wnt ligands in Wnt/LRP signaling and inhibits PC-3 tumor growth in vivo. FEBS Lett. 2011; 585:3120-3125.
