Identification Of Small Molecule TRABID Deubiquitinase Inhibitors By Computation-Based Virtual Screen
© Shi et al.; licensee BioMed Central Ltd. 2012
Received: 5 October 2011
Accepted: 19 April 2012
Published: 14 May 2012
Wnt/β-catenin-mediated gene transcription plays important roles in a wide range of biological and pathophysiological processes including tumorigenesis where β-catenin-mediated transcription activity frequently elevates. TRABID, a deubiquitinase, was shown to have a positive Wnt/β-catenin-mediated gene transcription and hence holds a promise as a putative anti-cancer target.
In this study, we used a combination of structure based virtual screening and an in vitro deubiquitinase (DUB) assay to identify several small molecules that inhibit TRABID DUB activity. However, these inhibitors failed to show inhibitory effects on β-catenin-mediated gene transcription. In addition, expression of TRABID shRNAs, wildtype TRABID, or the DUB activity-deficient mutant showed little effects on β-catenin-mediated gene transcription.
TRABID may not be a critical component in canonical Wnt/β-catenin signal transduction or that a minute amount of this protein is sufficient for its role in regulating Wnt activity.
The evolutionarily conserved Wnt signaling pathway regulates various developmental processes during embryogenesis and plays an important role for tissue homeostasis in adults. Wnt signaling pathway also plays an important role in tumorigenesis, particularly the formation of inherited and sporadic colorectal cancer as result of adenomatous polyposis coli (APC) mutation that leads to β-catenin accumulation in the nucleus [1–4]. Nuclear beta-catenin binds to and functions as a cofactor of lymphoid enhancer-binding factor (LEF-1)  and T cell factors (TCF)  to stimulate the transcription of Wnt target genes . Wnt-β-catenin signaling is essential in sustaining the cancer stem cell phenotype and is also involved in the transformation into malignant human squamous cell carcinomas . The loss of Wnt signaling pathway component APC in stem cells results in progressively growing neoplasia . Activation of Wnt/TCF pathway is also a determinant of lung adenocarcinoma metastasis to brain and bone. The phenotype of those metastatic derivatives of adenocarcinoma resembles the bronchioalveolar stem cells . Thus, it holds promises to prevent and/or treat cancers by targeting Wnt signaling pathway components.
TRAF-binding domain (TRABID), one of deubiquitination enzymes, was recently reported to specifically and positively regulate Wnt signaling pathway . TRABID is composed of a TRAF-binding domain in the C-terminus and three Zinc-finger (ZnF) motifs at the N-terminus . TRABID preferentially binds to lysine 63-linked polyubiquitin chains (K63 chains), but not K48-linked polyubiquitin chains (K48 chains), and specifically cleaves K63 chains . K63-linked ubiquitination was suggested to regulate substrate activity rather than protein stability [13–16]. The catalytic residues of TRABID reside in the OTU (ovarian tumor) domain within the TRAF-binding domain .
The OTU domain is conserved within the members of the OTU family of DUBs that include A20 and possess the cysteine protease activity . TRABID was also shown to bind to APC and may be responsible for its deubiquitination. However, the direct interaction between APC and TRABID was not detected . Knockdown of TRABID with RNAi resulted in downregulated expression of canonical Wnt target genes and decreased Wnt transcription activity, whereas its knockdown does not affect TNF-2 pathway . Epistasis analysis suggested that TRABID might act downstream of beta-catenin stabilization and affect the interaction of beta-catenin with LEF1 . Furthermore, TRABID heterozygosity suppressed the rough eye phenotype caused by ectopic Wingless expression in fly, but did not affect the phenotypes caused by the inhibition of Notch signaling or inhibition of EGF signaling. These results suggest that TRABID might be a potential drug target for controlling Wnt pathway activation .
In this study, we screened for small molecule inhibitors of TRABID by employing a combination of structure based virtual screening and an in vitro DUB assay. We searched for compounds in a chemical library from the National Cancer Institute that potentially bind to TRABID catalytic site based on the crystal structure of A20 catalytic domain . We identified several compounds that are able to inhibit the DUB activity of TRABID. However, these inhibitors failed to show inhibitory effects on Wnt activity. Furthermore, neither shRNAs that silenced TRABID efficiently nor overexpression of wild type (WT) TRABID or its DUB activity-deficient mutant showed inhibitory effects on canonical Wnt signaling activity.
Results and discussion
Searching on the surface of the equilibrated structure of TRABID OTU domain modeled from the crystal structure of A20 OTU revealed that there was a potential pocket adjacent to the catalytic center of OTU (Figure 1E). This pocket is part of a distal ubiquitin binding site at the catalytic site in A20 . Because of its spatial proximity to the catalytic site and its role in ubiquitin binding, it is reasonable to postulate that targeting this site may block the entrance of the substrates and hence inhibit the OTU DUB activity. Of note, during the molecular dynamics refinement, considerable conformational changes occurred in this site, making it more expanded comparing with the pocket size in the A20 crystal structure (Figure 1F).
Wnt signaling plays a crucial role in cancer cell growth and cancer stem cell maintenance, which are believe to be responsible for cancer metastasis and poor prognosis of chemotherapy. Thereby, regulation of Wnt signaling holds promise in chemotherapy of cancer. TRABID, one of K63-specific DUB enzymes, became an attractive candidate to achieve this task since it has been reported to regulate Wnt signal positively , as inhibition of TRABID presumably decreases Wnt signaling activity. However, in this study, we failed to confirm the notion that TRABID positively regulates canonical Wnt signaling based on several lines of evidence. Firstly, the compounds, which blocked TRABID DUB activity in vitro, failed to inhibit Wnt transcription activity in two tumor cell lines whereβcatenin-mediated transcription are known to be elevated. Secondly, the TRABID shRNAs failed to downregulate Wnt/β-catenin-mediated transcription activity or endogenous Wnt target gene Axin2 mRNA levels in these tumor cells. Thirdly, overexpression of WT and C443S mutant TRABID failed to show any effect on Wnt target transcription activity. We did observe that the siRNA described in the previous report  was able to inhibit Wnt reporter gene activity (Figure 4C-E). One possibility is that the effect of the siRNA reported in  is off target. A less likely possibility is that a minute amount of TRABID is sufficient for its regulation of Wnt signaling in these cells, and our shRNAs or inhibitors were not efficient enough to alter Wnt activity. Nevertheless, the fly data  have already suggested that TRABID might not be a core component of Wnt signaling. The observation that the TRABID inhibitors were able to inhibit cell growth suggests that TRABID may have a role in cell growth regulation. However, the mechanisms for this role of TRABID need to be further investigated.
We have used computation-based virtual screen to identify chemical inhibitors for TRABID, but failed to confirm that TRABID has a significant role in Wnt signaling. Nevertheless, inhibition of TRABID may inhibit cell growth.
Construction of TRABID expression vectors and shRNA expression vectors Vectors expressing GST-fused full length (FL) and N-terminal fragment (NT) (amino acides: 1–354) of human TRABID were constructed using pET-GST vector (EMD Chemicals Inc., Gibbstown, NJ). Proteins were expressed in bacteria and purified using glutathione sepharose 4B (Roche Diagnostics, Indianapolis, IN). FLAG-expression vectors for expression of the full-length and DUB-deficient (C443S) TRABID in mammalian cells were gifts from Dr. Hoanh Tran. FLAG-tagged proteins were immunoprecipitated using the M2 agarose affinity gel (Sigma-Aldrich, St. Louis, MO). TRABID DNA sequences were confirmed by DNA sequencing. TRABID shRNA targeting sequences were selected using the RNAi Codex algorithm (http://katahdin.cshl.org:9331_portal_scripts_main2.pl)  and the shRNAs were expressed using Migr-CMV-YFP-miR30 vector as previously reported . The shRNA target sequences are: shTrabd1, TGTCTCAACAAGCAGCAAAGT; shTrabd2, AGGAGCTAGGTAATGAGGAAC; shTrabd3, TGTCAGAACGTGGAATTAAGT; shTrabd4, TGATCATCCCAGACCTAATAA; shTrabd5, CTGGCACATATTCTTAGACGA; shTrabd6, TTGGAAAAGTCCGATTGCTCT. β-catenin siRNA  and TRABID siRNA  target sequence were described previously.
DUB enzymatic activity assay
For determining TRABID DUB enzymatic activity, hexa-K63 ubiquitin (0.01 ug/ul, Boston Biochem Inc., Cambridge, MA) and the compounds were incubated with TRABID immunoprecipitated from transfected HEK 293 T cells or prepared from bacteria for 3 hrs at 37°C in a DUB reaction buffer (50 mM HEPES at pH 7.4, 150 mM KCl, 10 mM DTT, 5% glycerol, 0.01% Triton X-100) . Samples were applied to SDS-PAGE under non-reducing conditions and then transferred onto members. Ubiquitin chains were detected using a rabbit anti-ubiquitin antibody (Bethyl Laboratories Inc., Montgomery, TX).For determining A20 DUB enzymatic activity, Penta-K48 ubiquitin (0.04 ug/ul, Boston Biochem Inc., Cambridge, MA), 1.5 uM of A20 catalytic domain (Boston Biochem Inc., Cambridge, MA), and various concentrations of compounds were incubated for 3 hrs at 37 C in an A20 DUB reaction buffer (50 mM HEPES at pH 8.0, 3 mM DTT, 0.01% Brij-35) .
Homology modeling and potential pocket identification There are two crystal structures of A20 OTU domains deposited in PDB with nearly identical sequences and structures. Because The OTU domain of A20, shared a high degree of amino acid homology with that of TRABID, its sturcutre is the ideal template for the homology modeling of TRABID OTU structure. We used one of the A20 OTU structures (from PDB ID 2VFJ ) as the template for homology modeling.. The homology model of TRABID was built by the Prime module of the Schrodinger molecular modeling package [22, 23]. A 3 ns molecular dynamics was performed by AMBER10  to equilibrate the initial model since there is considerable gap (15%) between the sequences of TRABID and A20 (Figure 1A). The final equilibrated structure is used for the further modeling. The ICM PocketFinder module based on an algorithm of utilizing a Van de Waals grid potential map and a carbon probe  was used to detect potential relevant pockets around the catalytic cysteine residue of the TRABID OTU domain.
Since there is no any knowledge about the characteristics and properties of the potential ligands, it is necessary to search all the compounds in the library to find potential hits. We thus carried out virtual screens of the National Cancer Institute (NCI) database (http://188.8.131.52/ncidb2) for chemical compounds that could potentially bind to the detected cavity. This database includes the coordinates of 250,251 small (M.W. < 1,000 Da) chemical compounds. In order to eliminate the large amount of irrelevant compounds we used Glide module [26, 27] in the Schrodinger package to perform hierarchical searches. The scoring function of High-throughput in Glide, which is less accurate but very fast, was used to perform the first round screening. Subsequently, the top 5000 compounds were selected and further tested with the normal scoring function in Glide. Finally, the most precise, but slower, scoring function in Glide  was used to perform the final screening with the top 500 compounds, and the top 200 compounds were requested from NCI and evaluated by the experimental assay.
Where ∆Gvw is the Van der Waals energy, ∆Ghb is the htdrogen bonding energy, ∆Gto is torsional energy, ∆Gel is electrostastic energy; and ∆sf is the surface energy.
TOPFLASH reporter assays
For testing the effect of compounds, cells were seeded in 48-well culture plates at a density of 0.2-0.3 million/ml, 250 ul/well and transfected with TOPFLASH or FOPFLASH and Renilla luciferase using the Lipofectamine Plus reagent (Invitrogen Corporation, Carlsbad, CA). The compounds were added 3 hrs later. For shRNA knockdown experiments, cells were transfected with an shRNA expressing plasmid and re-transfected with the Wnt reporter gene plasmids and the shRNA two days later. Wnt reporter gene activity was determined 24 hrs after the re-transfection and presented after normalization against the Renilla luciferase activity.
Cells were plated into 6 well plates at a density of 0.2 million/ml for 12–16 hrs. For compound inhibition experiments, cells were incubated with 20 μM of a compound for 24 hrs. For shRNA inhibition experiments, cells were transfected with the shRNAs and incubated for 72 hrs. Total RNAs were extracted using the RNAeasy mini kit (Qiagen, Valencia, CA) and transcripted into cDNAs. The levels of Wnt target gene mRNAs were determined using a MyiQ Single color real-time PCR detection system (Bio-Rad, Hercules, CA) with beta-actin as a control.
Cell growth assay
Cells were seeded at a density of 0.05 million/ml in 24-well culture plates and added with 20 μM compound. On indicated date, the numbers of living cells were counted using Guava Easycyte Mini Flow Cytometry System (Millipore Corporation, Billerica, MA).
We thank Xiaofeng Li and Yong Zhang for technical help. This work is supported by NIH grant to D.W. (CA132317).
- Nishisho I, Nakamura Y, Miyoshi Y, Miki Y, Ando H, Horii A, Koyama K, Utsunomiya J, Baba S, Hedge P: Mutations of chromosome 5q21 genes in FAP and colorectal cancer patients. Science. 1991, 253 (5020): 665-669. 10.1126/science.1651563.PubMedView ArticleGoogle Scholar
- Kinzler KW, Nilbert MC, Su LK, Vogelstein B, Bryan TM, Levy DB, Smith KJ, Preisinger AC, Hedge P, McKechnie D: Identification of FAP locus genes from chromosome 5q21. Science. 1991, 253 (5020): 661-665. 10.1126/science.1651562.PubMedView ArticleGoogle Scholar
- Korinek V, Barker N, Morin PJ, van Wichen D, de Weger R, Kinzler KW, Vogelstein B, Clevers H: Constitutive transcriptional activation by a beta-catenin-Tcf complex in APC−/− colon carcinoma. Science. 1997, 275 (5307): 1784-1787. 10.1126/science.275.5307.1784.PubMedView ArticleGoogle Scholar
- Groden J, Thliveris A, Samowitz W, Carlson M, Gelbert L, Albertsen H, Joslyn G, Stevens J, Spirio L, Robertson M: Identification and characterization of the familial adenomatous polyposis coli gene. Cell. 1991, 66 (3): 589-600. 10.1016/0092-8674(81)90021-0.PubMedView ArticleGoogle Scholar
- Behrens J, von Kries JP, Kuhl M, Bruhn L, Wedlich D, Grosschedl R, Birchmeier W: Functional interaction of beta-catenin with the transcription factor LEF-1. Nature. 1996, 382 (6592): 638-642. 10.1038/382638a0.PubMedView ArticleGoogle Scholar
- Molenaar M, van de Wetering M, Oosterwegel M, Peterson-Maduro J, Godsave S, Korinek V, Roose J, Destree O, Clevers H: XTcf-3 transcription factor mediates beta-catenin-induced axis formation in Xenopus embryos. Cell. 1996, 86 (3): 391-399. 10.1016/S0092-8674(00)80112-9.PubMedView ArticleGoogle Scholar
- Galceran J, Farinas I, Depew MJ, Clevers H, Grosschedl R: Wnt3a−/−−like phenotype and limb deficiency in Lef1(−/−)Tcf1(−/−) mice. Genes Dev. 1999, 13 (6): 709-717. 10.1101/gad.13.6.709.PubMedPubMed CentralView ArticleGoogle Scholar
- Malanchi I, Peinado H, Kassen D, Hussenet T, Metzger D, Chambon P, Huber M, Hohl D, Cano A, Birchmeier W: Cutaneous cancer stem cell maintenance is dependent on beta-catenin signalling. Nature. 2008, 452 (7187): 650-653. 10.1038/nature06835.PubMedView ArticleGoogle Scholar
- Barker N, Ridgway RA, van Es JH, van de Wetering M, Begthel H, van den Born M, Danenberg E, Clarke AR, Sansom OJ, Clevers H: Crypt stem cells as the cells-of-origin of intestinal cancer. Nature. 2009, 457 (7229): 608-611. 10.1038/nature07602.PubMedView ArticleGoogle Scholar
- Nguyen DX, Chiang AC, Zhang XH, Kim JY, Kris MG, Ladanyi M, Gerald WL, Massague J: WNT/TCF signaling through LEF1 and HOXB9 mediates lung adenocarcinoma metastasis. Cell. 2009, 138 (1): 51-62. 10.1016/j.cell.2009.04.030.PubMedPubMed CentralView ArticleGoogle Scholar
- Tran H, Hamada F, Schwarz-Romond T, Bienz M: Trabid, a new positive regulator of Wnt-induced transcription with preference for binding and cleaving K63-linked ubiquitin chains. Genes Dev. 2008, 22 (4): 528-542. 10.1101/gad.463208.PubMedPubMed CentralView ArticleGoogle Scholar
- Evans PC, Taylor ER, Coadwell J, Heyninck K, Beyaert R, Kilshaw PJ: Isolation and characterization of two novel A20-like proteins. Biochem J. 2001, 357 (Pt 3): 617-623.PubMedPubMed CentralView ArticleGoogle Scholar
- Deng L, Wang C, Spencer E, Yang L, Braun A, You J, Slaughter C, Pickart C, Chen ZJ: Activation of the IkappaB kinase complex by TRAF6 requires a dimeric ubiquitin-conjugating enzyme complex and a unique polyubiquitin chain. Cell. 2000, 103 (2): 351-361. 10.1016/S0092-8674(00)00126-4.PubMedView ArticleGoogle Scholar
- Wang C, Deng L, Hong M, Akkaraju GR, Inoue J, Chen ZJ: TAK1 is a ubiquitin-dependent kinase of MKK and IKK. Nature. 2001, 412 (6844): 346-351. 10.1038/35085597.PubMedView ArticleGoogle Scholar
- Zhou H, Wertz I, O'Rourke K, Ultsch M, Seshagiri S, Eby M, Xiao W, Dixit VM: Bcl10 activates the NF-kappaB pathway through ubiquitination of NEMO. Nature. 2004, 427 (6970): 167-171. 10.1038/nature02273.PubMedView ArticleGoogle Scholar
- Xia ZP, Sun L, Chen X, Pineda G, Jiang X, Adhikari A, Zeng W, Chen ZJ: Direct activation of protein kinases by unanchored polyubiquitin chains. Nature. 2009, 461 (7260): 114-119. 10.1038/nature08247.PubMedPubMed CentralView ArticleGoogle Scholar
- Komander D, Barford D: Structure of the A20 OTU domain and mechanistic insights into deubiquitination. Biochem J. 2008, 409 (1): 77-85. 10.1042/BJ20071399.PubMedView ArticleGoogle Scholar
- Shin KJ, Wall EA, Zavzavadjian JR, Santat LA, Liu J, Hwang JI, Rebres R, Roach T, Seaman W, Simon MI: A single lentiviral vector platform for microRNA-based conditional RNA interference and coordinated transgene expression. Proc Natl Acad Sci U S A. 2006, 103 (37): 13759-13764. 10.1073/pnas.0606179103.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang Y, Tang W, Jones MC, Xu W, Halene S, Wu D: Different roles of G protein subunits beta1 and beta2 in neutrophil function revealed by gene expression silencing in primary mouse neutrophils. J Biol Chem. 2010, 285 (32): 24805-24814. 10.1074/jbc.M110.142885.PubMedPubMed CentralView ArticleGoogle Scholar
- Pan W, Choi SC, Wang H, Qin Y, Volpicelli-Daley L, Swan L, Lucast L, Khoo C, Zhang X, Li L: Wnt3a-mediated formation of phosphatidylinositol 4,5-bisphosphate regulates LRP6 phosphorylation. Science. 2008, 321 (5894): 1350-1353. 10.1126/science.1160741.PubMedPubMed CentralView ArticleGoogle Scholar
- Wertz IE, O'Rourke KM, Zhou H, Eby M, Aravind L, Seshagiri S, Wu P, Wiesmann C, Baker R, Boone DL: De-ubiquitination and ubiquitin ligase domains of A20 downregulate NF-kappaB signalling. Nature. 2004, 430 (7000): 694-699. 10.1038/nature02794.PubMedView ArticleGoogle Scholar
- Jacobson MP, Friesner RA, Xiang ZX, Honig B: On the role of the crystal environment in determining protein side-chain conformations. J Mol Biol. 2002, 320 (3): 597-608. 10.1016/S0022-2836(02)00470-9.PubMedView ArticleGoogle Scholar
- Jacobson MP, Pincus DL, Rapp CS, Day TJF, Honig B, Shaw DE, Friesner RA: A hierarchical approach to all-atom protein loop prediction. Proteins. 2004, 55 (2): 351-367. 10.1002/prot.10613.PubMedView ArticleGoogle Scholar
- Case DA, Darden TA, Cheatham ITE, Simmerling CL, Wang J, Duke RE, Luo R, Crowley M, Walker RC, Zhang W: AMBER10. In. San Francisco, CA. University of California, San Francisco. 2008Google Scholar
- Abagyan R, An JH, Totrov M: Pocketome via comprehensive identification and classification of ligand binding envelopes. Mol Cell Proteomics. 2005, 4 (6): 752-761. 10.1074/mcp.M400159-MCP200.PubMedView ArticleGoogle Scholar
- Friesner RA, Banks JL, Murphy RB, Halgren TA, Klicic JJ, Mainz DT, Repasky MP, Knoll EH, Shelley M, Perry JK: Glide: A new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. Journal of Medicinal Chemistry. 2004, 47 (7): 1739-1749. 10.1021/jm0306430.PubMedView ArticleGoogle Scholar
- Halgren TA, Murphy RB, Friesner RA, Beard HS, Frye LL, Pollard WT, Banks JL: Glide: A new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. Journal of Medicinal Chemistry. 2004, 47 (7): 1750-1759. 10.1021/jm030644s.PubMedView ArticleGoogle Scholar
- Friesner RA, Murphy RB, Repasky MP, Frye LL, Greenwood JR, Halgren TA, Sanschagrin PC, Mainz DT: Extra precision glide: Docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes. Journal of Medicinal Chemistry. 2006, 49 (21): 6177-6196. 10.1021/jm051256o.PubMedView ArticleGoogle Scholar
- Schapira M, Totrov M, Abagyan R: Prediction of the binding energy for small molecules, peptides and proteins. J Mol Recognit. 1999, 12 (3): 177-190. 10.1002/(SICI)1099-1352(199905/06)12:3<177::AID-JMR451>3.0.CO;2-Z.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.