- Research article
- Open Access
Identification Of Small Molecule TRABID Deubiquitinase Inhibitors By Computation-Based Virtual Screen
BMC Chemical Biologyvolume 12, Article number: 4 (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
To screen for compounds that can inhibit TRABID DUB activity in vitro, we established an assay to measure the DUB activity of TRABID. Consistent with the previous report by Tran et al. , we observed that recombinant TRABID proteins purified from an E coli expression system or pulled down from HEK293 cells overexpressing Flag-TRABID by an Flag antibody were able to specifically cleave hexa-K63 (Figure 1B,C), but not penta-K48 (data not shown), ubiquitin chains. We also tested a TRABID DUB-deficient mutant containing a substitution mutation (C443A), a residue located in the OTU domain and critical for its catalytic activity . The mutation abrogated the ability of TRABID to cleave the Hexa-K63 ubiquitin substrate (Figure 1C), confirming the importance of this residue for the DUB activity. Homology modeling of TRABID OTU domain based on the crystal structure of A20 OTU domain  revealed that the common cysteine catalytic site was well conserved in TRABID. This site is characterized by a cysteine residue that forms an electrostatic network with a histidine and a residue with an acidic side chain. In TRABID, these 3 residues are C443, H596 and D410, respectively, based on the homology model (Figure 1D).
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).
We performed hierarchical virtual screening towards this pocket as described in the “Materials and Methods”, and requested 200 compounds that were ranked highest in the virtual screen and were provided with 125 by NCI. These 125 compounds were screened for their abilities to inhibit the cleavage of Hexa-K63 ubiquitin by TRABID. Seventeen of these compounds showed more than 50 percent (denote with ++ in Additional file 1. Table S1) and five showed between 25–50 percent (denote with + in Additional file 1. Table S1) inhibition of TRABID DUB activity. Figure 2A shows a representative Western blot used in the TRABID DUB activity assay. All of the positive compounds were retested in a dose-dependent experiment. Figure 2B shows the data for the two of the strongest inhibitors, which showed the IC50 values of about 3 μM. Neither compound showed significant inhibition of A20 up to 30 μM, though both inhibited A20 at 100 μM (Figure 2C-D). Compound NSC60650 as well as a number of other compounds shown in Additional file 1. Table S1, which share similar chemical structures to NSC112200 and NSC267309 (Figure 2E), did not inhibit TRABID DUB activity (Figure 2B), suggesting that the two hydroxyl groups and their locations on the benzene ring may be important for the inhibitory activity. The molecular modeling of the bindings of NSC112200, NSC267309 and NSC60650 to TRABID OTU domain suggests that while all three compounds can form electrostatic interactions with the side chain of E522 and the backbone NH of T556, both NSC112200 and NSC267309 can form an additional hydrogen bond with the side chain of S491 (Figure 2F-H). Moreover, NSC112200 can form another hydrogen bond with the side chain of S520, and its one of methyl groups provides extra hydrophobic interaction with the carbon atoms located at the side chain of R557 (Figure 2F). The rescoring of binding free energies by ICM shows that the binding free energies of NSC112200 and NSC267309 with TRABID OTU are quite similar (−20.28 kcal/mol and −20.03 kcal/mol, respectively), while the binding free energy of NSC60650 with TRABID OTU is −16.89 kcal/mol, which is 3.39 kcal/mol higher than NSC112200. On the other hand, the calculated binding free energies of NSC112200 and NSC267309 with A20 OTU are −15.23 kcal/mol and −13.43 kcal/mol, respectively. As the high binding free energy means a lower binding affinity, our modeling predictions are consistent with our experimental data. TRABID was shown to play a positive role in canonical Wnt signaling, and its DUB activity appeared to be essential for this function . We hence tested the TRABID DUB inhibiting compounds for their effects on canonical Wnt signaling. SW480 or HCT116 colorectal cancer cell lines, which were used in the previous study for TRABID’s role in Wnt signaling , were transfected with a Wnt reporter gene construct TOPFLASH and incubated with compounds for 24 hrs. Because these two cell lines contain mutations in the APC and β-catenin genes, respectively, there was high reporter gene activity in these cells (Figure 3A-B). However, neither NSC112200 nor NSC267309 inhibited the reporter gene activity (Figure 3A-B). As controls knockdown β-catenin resulted in significant reduction in Wnt reporter gene activity (Figure 3A,B). We also tested other TRABID inhibitors in Additional file 1. Table S1, and none of them inhibited the reporter gene activity (data not shown). We then examined two endogenous Wnt target gene expression. Neither Axin2 nor c-Myc expression in SW480 cells was affected by NSC112200 or NSC267309 treatment (Figure 3C). These results together suggest that TRABID might not be involved in Wnt$β-catenin-mediated gene transcription regulation. To further evaluate the role of TRABID in the regulation of Wnt/β-catenin-mediated gene transcription, we generated 6 constructs that produce six independent shRNAs for TRABID. When these shRNAs were coexpressed with TRABID cDNA in HEK293 cells, they showed varying efficiencies in suppressing TRABID expression with shTrbd5 and 1 being the most efficiency and shTrbd2 and 4 showing little effects (Figure 4A). We also evaluated their knockdown efficiencies using quantitative RT-PCR analysis of RNAs from YFP-positive HCT116 transfectants (shRNA and YFP were in the same transcript) (Figure 4B). However, expression of these shRNAs had no inhibitory effect on Wnt reporter gene expression in either SW480 or HCT-116 cells (Figure 4C,D). These shRNAs also did not show inhibitory effect on the expression of endogenous Wnt target gene Axin2 in SW480 cells (Figure 4E). Furthermore, overexpression of wildtype TRABID or its DUB-deficient mutant had little effects on Wnt reporter gene expression in either SW480 or HCT-116 cells (Figure 4F,G) or in HEK293 cells (Figure 4H). Thus, we failed to demonstrate that TRABID has an important role in Wnt target gene expression regulation in these two cancer cells.
Nevertheless, both TRABID inhibitors NSC112200 and NSC267309, but not the non-active analog NSC60650, inhibited the growth of not only the colorectal tumor cell lines HCT-116 and SW480, but also two para-normal cell lines NIH3T3 and HEK293 at 20 μM (Figure 5). These compounds showed no effects on the growth of these cells at 5 μM (data not shown). Because of the lack of effect of Compound NSC60650 on cell growth, we believe that the effects of Compounds NSC112200 and NSC267309 are likely due to their inhibition of the deubiquitinase activity.
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://184.108.40.206/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.
In order to further interpret the inhibitory mechanism, the binding free energies of two strong TRABID deubiquitinase inhibitors, NSC112200 and NSC267309, identified in the screen were reevaluated by using the ICM empirical function, which is fast and fairly accurate in predicting relative binding free energies . The binding free energy of NSC60650, which shares similar chemical structure with NSC112200 and NSC267309 but shows no inhibition, was also recalculated for comparison. The overall complexes were optimized and the binding free were calculated using Equation:
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).
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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. 2008
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.
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.
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.
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.
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.
We thank Xiaofeng Li and Yong Zhang for technical help. This work is supported by NIH grant to D.W. (CA132317).
T. S. carried out biological experiments and drafted the manuscript. J.B., N.X.W. and J.Z. performed the virtual screen, modeling and manuscript preparation. D.W. participated in experimental design and drafted the manuscript. All authors read and approved the final manuscript
Tong Shi, Ju Bao, Nick X Wang contributed equally to this work.