Structure, function and disease relevance of Wnt inhibitory factor 1, a secreted protein controlling the Wnt and hedgehog pathways.

Wnts and Hedgehogs (Hh) are large, lipid-modified extracellular morphogens that play key roles in embryonic development and stem cell proliferation of Metazoa. Both morphogens signal through heptahelical Frizzled-type receptors of the G-Protein Coupled Receptor family and there are several other similarities that suggest a common evolutionary origin of the Hh and Wnt pathways. There is evidence that the secreted protein, Wnt inhibitory factor 1 (WIF1) modulates the activity of both Wnts and Hhs and may thus contribute to the intertwining of these pathways. In this article, we review the structure, evolution, molecular interactions and functions of WIF1 with major emphasis on its role in carcinogenesis.

Arginines of the CGN codon family are Achilles' heels of cancer genes.

Recent studies have revealed that arginine is the most favorable target of amino acid alteration in most cancer types and it has been suggested that the high preference for arginine mutations reflects the critical roles of this amino acid in the function of proteins. High rates of mutations of arginine residues in cancer, however, might also be due to increased mutability of arginine codons of the CGN family as the CpG dinucleotides of these codons may be methylated. In the present work we have analyzed spectra of single base substitutions of cancer genes (oncogenes, tumor suppressor genes) and passenger genes in cancer tissues to assess the contributions of CpG hypermutability and selection to arginine mutations. Our studies have shown that arginines encoded by the CGN codon family display higher rates of mutation in both cancer genes and passenger genes than arginine codons AGA and AGG that are devoid of CpG dinucleotide, suggesting that the predominance of arginine mutations in cancer is primarily due to CpG hypermutability, rather than selection for arginine replacement. Nevertheless, our results also suggest that CGN codons for arginines may serve as Achilles' heels of cancer genes. CpG hypermutability of key arginines of proto-oncogenes, leading to high rates of recurrence of driver mutations, contributes significantly to carcinogenesis. Similarly, our results indicate that hypermutability of the CpG dinucleotide of CGA codons (converting them to TGA stop codons) contributes significantly to recurrent truncation and inactivation of tumor suppressor genes.

Evolution of termination codons of proteins and the TAG-TGA paradox.

In most eukaryotes and prokaryotes TGA is used at a significantly higher frequency than TAG as termination codon of protein-coding genes. Although this phenomenon has been recognized several years ago, there is no generally accepted explanation for the TAG-TGA paradox. Our analyses of human mutation data revealed that out of the eighteen sense codons that can give rise to a nonsense codon by single base substitution, the CGA codon is exceptional: it gives rise to the TGA stop codon at an order of magnitude higher rate than the other codons. Here we propose that the TAG-TGA paradox is due to methylation and hypermutabilty of CpG dinucleotides. In harmony with this explanation, we show that the coding genomes of organisms with strong CpG methylation have a significant bias for TGA whereas those from organisms that lack CpG methylation use TGA and TAG termination codons with similar probability.

Use of Publication Dynamics to Distinguish Cancer Genes and Bystander Genes.

de Magalhães has shown recently that most human genes have several papers in PubMed mentioning cancer, leading the author to suggest that every gene is associated with cancer, a conclusion that contradicts the widely held view that cancer is driven by a limited number of cancer genes, whereas the majority of genes are just bystanders in carcinogenesis. We have analyzed PubMed to decide whether publication metrics supports the distinction of bystander genes and cancer genes. The dynamics of publications on known cancer genes followed a similar pattern: seminal discoveries triggered a burst of cancer-related publications that validated and expanded the discovery, resulting in a rise both in the number and proportion of cancer-related publications on that gene. The dynamics of publications on bystander genes was markedly different. Although there is a slow but continuous time-dependent rise in the proportion of papers mentioning cancer, this phenomenon just reflects the increasing publication bias that favors cancer research. Despite this bias, the proportion of cancer papers on bystander genes remains low. Here, we show that the distinctive publication dynamics of cancer genes and bystander genes may be used for the identification of cancer genes.

Wnt Inhibitory Factor 1 Binds to and Inhibits the Activity of Sonic Hedgehog.

The hedgehog (Hh) and Wnt pathways, crucial for the embryonic development and stem cell proliferation of Metazoa, have long been known to have similarities that argue for their common evolutionary origin. A surprising additional similarity of the two pathways came with the discovery that WIF1 proteins are involved in the regulation of both the Wnt and Hh pathways. Originally, WIF1 (Wnt Inhibitory Factor 1) was identified as a Wnt antagonist of vertebrates, but subsequent studies have shown that in Drosophila, the WIF1 ortholog serves primarily to control the distribution of Hh. In the present, work we have characterized the interaction of the human WIF1 protein with human sonic hedgehog (Shh) using Surface Plasmon Resonance spectroscopy and reporter assays monitoring the signaling activity of human Shh. Our studies have shown that human WIF1 protein binds human Shh with high affinity and inhibits its signaling activity efficiently. Our observation that the human WIF1 protein is a potent antagonist of human Shh suggests that the known tumor suppressor activity of WIF1 may not be ascribed only to its role as a Wnt inhibitor.

Use of signals of positive and negative selection to distinguish cancer genes and passenger genes.

A major goal of cancer genomics is to identify all genes that play critical roles in carcinogenesis. Most approaches focused on genes positively selected for mutations that drive carcinogenesis and neglected the role of negative selection. Some studies have actually concluded that negative selection has no role in cancer evolution. We have re-examined the role of negative selection in tumor evolution through the analysis of the patterns of somatic mutations affecting the coding sequences of human genes. Our analyses have confirmed that tumor suppressor genes are positively selected for inactivating mutations, oncogenes, however, were found to display signals of both negative selection for inactivating mutations and positive selection for activating mutations. Significantly, we have identified numerous human genes that show signs of strong negative selection during tumor evolution, suggesting that their functional integrity is essential for the growth and survival of tumor cells.

The DNA in the cells of the human body is usually copied correctly when a cell divides. However, errors (mutations) are sometimes introduced during the copying process. Although the majority of mutations have no major impact on cells, many mutations are harmful: they decrease the ability of cells to survive. There are, however, mutations that can lead to cells dividing more frequently or gaining the ability to spread, which can lead to cancer. These mutations are known as 'driver mutations' because they drive the growth of tumors. Since such 'driver mutations' provide a growth advantage to tumor cells, they are subject to positive selection, this is, their frequency in the tumor increases over time. Because of their selective advantage, driver mutations accumulate at significantly higher rates than the neutral 'passenger mutations' that do not play a role in tumor growth. Genes that carry driver mutations are called driver genes, while genes that carry only passenger mutations are known as passenger genes. Certain genes, however, do not fit into either category. For example, some genes that are essential for tumor growth must get rid of harmful mutations to maintain activity. Mutations of such 'tumor essential genes' are thus subject to 'negative' or 'purifying selection'. A major goal of cancer research is to identify genes that play critical roles in tumor growth. Earlier studies have identified numerous driver genes positively selected for driver mutations, exploiting the fact that driver genes show significantly higher mutation rates than passenger genes. Identification of tumor essential genes, however, is inherently more difficult since the paucity of mutations of negatively selected genes hinders the analysis of the mutation data. The failure to provide convincing evidence for negative selection in tumors has led to suggestions that it has no role in cancer evolution. Bányai et al. used a novel approach to address the question of whether negative selection occurs in cancer. Based on characteristic differences in the patterns of mutations in cancer they distinguished clusters of passenger genes, driver genes and tumor essential genes. The group of tumor essential genes includes genes that serve to satisfy the increased demand of rapidly dividing tumor cells for nutrients' and genes that are essential for cell migration and metastasis (the spread of cancer cells to other areas of the body). The tumor essential genes that Bányai et al. identified may prove to be valuable targets for cancer therapy, illustrating the importance of genome sequencing in cancer research. Identification of additional tumor essential genes is, however, hindered by the fact that they are likely to have low levels of mutations, which can exclude them from meaningful analyses. Progress with genomic sequencing of tumors is expected to overcome this limitation and help identify additional genes that are essential for cancer growth.

Absence of the Tks4 Scaffold Protein Induces Epithelial-Mesenchymal Transition-Like Changes in Human Colon Cancer Cells.

Epithelial to mesenchymal transition (EMT) is a multipurpose process involved in wound healing, development, and certain pathological processes, such as metastasis formation. The Tks4 scaffold protein has been implicated in cancer progression; however, its role in oncogenesis is not well defined. In this study, the function of Tks4 was investigated in HCT116 colon cancer cells by knocking the protein out using the CRISPR/Cas9 system. Surprisingly, the absence of Tks4 induced significant changes in cell morphology, motility, adhesion and expression, and localization of E-cadherin, which are all considered as hallmarks of EMT. In agreement with these findings, the marked appearance of fibronectin, a marker of the mesenchymal phenotype, was also observed in Tks4-KO cells. Analysis of the expression of well-known EMT transcription factors revealed that Snail2 was strongly overexpressed in cells lacking Tks4. Tks4-KO cells showed increased motility and decreased cell-cell attachment. Collagen matrix invasion assays demonstrated the abundance of invasive solitary cells. Finally, the reintroduction of Tks4 protein in the Tks4-KO cells restored the expression levels of relevant key transcription factors, suggesting that the Tks4 scaffold protein has a specific and novel role in EMT regulation and cancer progression.

Identification and correction of abnormal, incomplete and mispredicted proteins in public databases.

BACKGROUND: Despite significant improvements in computational annotation of genomes, sequences of abnormal, incomplete or incorrectly predicted genes and proteins remain abundant in public databases. Since the majority of incomplete, abnormal or mispredicted entries are not annotated as such, these errors seriously affect the reliability of these databases. Here we describe the MisPred approach that may provide an efficient means for the quality control of databases. The current version of the MisPred approach uses five distinct routines for identifying abnormal, incomplete or mispredicted entries based on the principle that a sequence is likely to be incorrect if some of its features conflict with our current knowledge about protein-coding genes and proteins: (i) conflict between the predicted subcellular localization of proteins and the absence of the corresponding sequence signals; (ii) presence of extracellular and cytoplasmic domains and the absence of transmembrane segments; (iii) co-occurrence of extracellular and nuclear domains; (iv) violation of domain integrity; (v) chimeras encoded by two or more genes located on different chromosomes. RESULTS: Analyses of predicted EnsEMBL protein sequences of nine deuterostome (Homo sapiens, Mus musculus, Rattus norvegicus, Monodelphis domestica, Gallus gallus, Xenopus tropicalis, Fugu rubripes, Danio rerio and Ciona intestinalis) and two protostome species (Caenorhabditis elegans and Drosophila melanogaster) have revealed that the absence of expected signal peptides and violation of domain integrity account for the majority of mispredictions. Analyses of sequences predicted by NCBI's GNOMON annotation pipeline show that the rates of mispredictions are comparable to those of EnsEMBL. Interestingly, even the manually curated UniProtKB/Swiss-Prot dataset is contaminated with mispredicted or abnormal proteins, although to a much lesser extent than UniProtKB/TrEMBL or the EnsEMBL or GNOMON-predicted entries. CONCLUSION: MisPred works efficiently in identifying errors in predictions generated by the most reliable gene prediction tools such as the EnsEMBL and NCBI's GNOMON pipelines and also guides the correction of errors. We suggest that application of the MisPred approach will significantly improve the quality of gene predictions and the associated databases.

Collagen-binding domains of gelatinase A and thrombospondin-derived peptides impede endocytic clearance of active gelatinase A and promote HT1080 fibrosarcoma cell invasion.

Gelatinase A (matrix metalloproteinase-2, MMP-2) binds to several proteins through its collagen-binding domains (CBDs). Surface plasmon resonance analysis revealed a strong interaction between CBD123 and thrombospondin-1 (TSP-1), with a K(D) value of 2x10(-9) M. CBD123, as well as individual domains, behave as competitive inhibitors of the TSP-1-directed endocytic clearance of active MMP-2, but not of its latent form, by HT1080 fibrosarcoma cells. Enhanced level of active MMP-2 in conditioned medium was associated to increased matrigel invasion. Similarly, GGWSHWSPWSS and GGWSHW peptides, as tryptophan-rich peptides within properdin-repeat motifs (TSRs) of TSP-1, promoted MMP-2 accumulation and cell invasiveness. Our data document the importance of TSP-1 in promoting MMP-2-mediated cancer cell invasion through interaction between CBDs of the enzyme and TSRs motifs of TSP-1.

Morphological Stasis and Proteome Innovation in Cephalochordates.

Lancelets, extant representatives of basal chordates, are prototypic examples of evolutionary stasis; they preserved a morphology and body-plan most similar to the fossil chordates from the early Cambrian. Such a low level of morphological evolution is in harmony with a low rate of amino acid substitution; cephalochordate proteins were shown to evolve slower than those of the slowest evolving vertebrate, the elephant shark. Surprisingly, a study comparing the predicted proteomes of Chinese amphioxus, Branchiostoma belcheri and the Florida amphioxus, Branchiostoma floridae has led to the conclusion that the rate of creation of novel domain combinations is orders of magnitude greater in lancelets than in any other Metazoa, a finding that contradicts the notion that high rates of protein innovation are usually associated with major evolutionary innovations. Our earlier studies on a representative sample of proteins have provided evidence suggesting that the differences in the domain architectures of predicted proteins of these two lancelet species reflect annotation errors, rather than true innovations. In the present work, we have extended these studies to include a larger sample of genes and two additional lancelet species, Asymmetron lucayanum and Branchiostoma lanceolatum. These analyses have confirmed that the domain architecture differences of orthologous proteins of the four lancelet species are because of errors of gene prediction, the error rate in the given species being inversely related to the quality of the transcriptome dataset that was used to aid gene prediction.

NMR structure of the WIF domain of the human Wnt-inhibitory factor-1.

The human Wnt-binding protein Wnt-inhibitory factor-1 (WIF-1) comprises an N-terminal WIF module followed by five EGF-like repeats. Here we report the three-dimensional structure of the WIF domain of WIF-1 determined by NMR spectroscopy. The fold consists of an eight-stranded beta-sandwich reminiscent of the immunoglobulin fold. Residual detergent (Brij-35) used in the refolding protocol was found to bind tightly to the WIF domain. The binding site was identified by intermolecular nuclear Overhauser effects observed between the WIF domain and the alkyl chain of the detergent. The results point to a possible role of WIF domains as a recognition motif of Wnt and Drosophila Hedgehog proteins that are activated by palmitoylation.

Putative extremely high rate of proteome innovation in lancelets might be explained by high rate of gene prediction errors.

A recent analysis of the genomes of Chinese and Florida lancelets has concluded that the rate of creation of novel protein domain combinations is orders of magnitude greater in lancelets than in other metazoa and it was suggested that continuous activity of transposable elements in lancelets is responsible for this increased rate of protein innovation. Since morphologically Chinese and Florida lancelets are highly conserved, this finding would contradict the observation that high rates of protein innovation are usually associated with major evolutionary innovations. Here we show that the conclusion that the rate of proteome innovation is exceptionally high in lancelets may be unjustified: the differences observed in domain architectures of orthologous proteins of different amphioxus species probably reflect high rates of gene prediction errors rather than true innovation.

Modules, multidomain proteins and organismic complexity.

Originally the term 'protein module' was coined to distinguish mobile domains that frequently occur as building blocks of diverse multidomain proteins from 'static' domains that usually exist only as stand-alone units of single-domain proteins. Despite the widespread use of the term 'mobile domain', the distinction between static and mobile domains is rather vague as it is not easy to quantify the mobility of domains. In the present work we show that the most appropriate measure of the mobility of domains is the number of types of local environments in which a given domain is present. Ranking of domains with respect to this parameter in different evolutionary lineages highlighted marked differences in the propensity of domains to form multidomain proteins. Our analyses have also shown that there is a correlation between domain size and domain mobility: smaller domains are more likely to be used in the construction of multidomain proteins, whereas larger domains are more likely to be static, stand-alone domains. It is also shown that shuffling of a limited set of modules was facilitated by intronic recombination in the metazoan lineage and this has contributed significantly to the emergence of novel complex multidomain proteins, novel functions and increased organismic complexity of metazoa.

Hydroxamate-based peptide inhibitors of matrix metalloprotease 2.

There is major interest in designing inhibitors for matrix metalloproteinase 2 (MMP-2, gelatinase A) since this enzyme is known to be involved in pathological processes such as tumor invasion or rheumatoid arthritis. The majority of MMP-2 inhibitor candidate drugs block the active site of MMP-2 by binding to its catalytic Zn2+ ion through a chelating (hydroxamate, sulphonate etc.) group. Despite the general interest in designing MMP-2 inhibitors, the results with many of the drug candidates were disappointing, their failure was usually explained by cross-reactions with other MMPs. One way to enhance MMP-2 selectivity is to design inhibitors that interact with both the active site and exosites such as the fibronectin type II (FN2) domains of the enzyme. In the present work, we have examined the inhibitory potential and MMP-2 selectivity of hydroxamates of three groups of peptides known to bind to the collagen-binding FN2 domains of MMP-2. The first type of peptides consisted of collagen-like (Pro-Pro-Gly)(n) repeats, peptides of the second group were identified from a random 15-mer phage display library based on their binding to immobilized FN2 domains of MMP-2. A hydroxamate of peptide p33-42, known to bind to the third FN2 domain of MMP-2 has also been tested. Our studies have shown that these compounds inhibited MMP-2 with IC50 values of 10-100 microM. The fact that their inhibitory potential was nearly identical for MMP-2del, a recombinant version of MMP-2 that lacks the FN2 domains, suggests that inhibition is not mediated by their binding to FN2 domains. It seems likely that the failure to exploit interaction with the FN2 domains is due to the fact that the FN2 domains and the catalytic domain of MMP-2 tumble independently, therefore only a tiny fraction of the conformational isomers can bind peptide hydroxamates via both the active site and the FN2 domain(s).

Wnts grasp the WIF domain of Wnt Inhibitory Factor 1 at two distinct binding sites.

Wnts have a structure resembling a hand with "thumb" and "index" fingers that grasp the cysteine rich domains of Frizzled receptors at two distinct binding sites. In the present work we show that the WIF domain of Wnt Inhibitory Factor 1 is also bound by Wnts at two sites. Using C-terminal domains of Wnt5a and Wnt7a and arginine-scanning mutagenesis of the WIF domain we demonstrate that, whereas the N-terminal, lipid-modified "thumb" of Wnts interacts with the alkyl-binding site of the WIF domain, the C-terminal domain of Wnts (Wnt-CTD) binds to a surface on the opposite side of the WIF domain.

Evidence that human genes of modular proteins have retained significantly more ancestral introns than their fly or worm orthologues.

Comparison of the exon-intron structures of human, fly and worm orthologues of mosaic genes assembled from class 1-1 modules by exon-shuffling has revealed that human genes retained significantly more of the original inter-module introns than their protostome orthologues. It is suggested that the much higher rate of intron loss in the worm- and insect lineages than in the chordate lineage reflects their greater tendency for genome compaction.

Characterization of a Wnt-binding site of the WIF-domain of Wnt inhibitory factor-1.

A Wnt-binding site of the WIF-domain of Wnt inhibitory factor-1 was localized by structure-guided arginine-scanning mutagenesis in combination with surface plasmon resonance assays. Our observation that substitution of some residues of WIF resulted in an increased affinity for Wnt5a, but decreased affinity for Wnt3a, suggests that these residues may define the specificity spectrum of WIF for Wnts. These results hold promise for a more specific targeting of Wnt family members with WIF variants in various forms of cancer.

Reassessing domain architecture evolution of metazoan proteins: major impact of errors caused by confusing paralogs and epaktologs.

In the accompanying paper (Nagy, Szláma, Szarka, Trexler, Bányai, Patthy, Reassessing Domain Architecture Evolution of Metazoan Proteins: Major Impact of Gene Prediction Errors) we showed that in the case of UniProtKB/TrEMBL, RefSeq, EnsEMBL and NCBI's GNOMON predicted protein sequences of Metazoan species the contribution of erroneous (incomplete, abnormal, mispredicted) sequences to domain architecture (DA) differences of orthologous proteins might be greater than those of true gene rearrangements. Based on these findings, we suggest that earlier genome-scale studies based on comparison of predicted (frequently mispredicted) protein sequences may have led to some erroneous conclusions about the evolution of novel domain architectures of multidomain proteins. In this manuscript we examine the impact of confusing paralogous and epaktologous multidomain proteins (i.e., those that are related only through the independent acquisition of the same domain types) on conclusions drawn about DA evolution of multidomain proteins in Metazoa. To estimate the contribution of this type of error we have used as reference UniProtKB/Swiss-Prot sequences from protein families with well-characterized evolutionary histories. We have used two types of paralogy-group construction procedures and monitored the impact of various parameters on the separation of true paralogs from epaktologs on correctly annotated Swiss-Prot entries of multidomain proteins. Our studies have shown that, although public protein family databases are contaminated with epaktologs, analysis of the structure of sequence similarity networks of multidomain proteins provides an efficient means for the separation of epaktologs and paralogs. We have also demonstrated that contamination of protein families with epaktologs increases the apparent rate of DA change and introduces a bias in DA differences in as much as it increases the proportion of terminal over internal DA differences.We have shown that confusing paralogous and epaktologous multidomain proteins significantly increases the apparent rate of DA change in Metazoa and introduces a positional bias in favor of terminal over internal DA changes. Our findings caution that earlier studies based on analysis of datasets of protein families that were contaminated with epaktologs may have led to some erroneous conclusions about the evolution of novel domain architectures of multidomain proteins. A reassessment of the DA evolution of multidomain proteins is presented in an accompanying paper [1].

Correction: Nagy, A., et al. Reassessing Domain Architecture Evolution of Metazoan Proteins: Major Impact of Gene Prediction Errors. Genes 2011, 2, 449-501.

We found some errors in the published versions of Figure S2, Figure S3 and Figure S8 of our paper [1]. The correct Figures are presented below. [...].

Reassessing domain architecture evolution of metazoan proteins: major impact of gene prediction errors.

In view of the fact that appearance of novel protein domain architectures (DA) is closely associated with biological innovations, there is a growing interest in the genome-scale reconstruction of the evolutionary history of the domain architectures of multidomain proteins. In such analyses, however, it is usually ignored that a significant proportion of Metazoan sequences analyzed is mispredicted and that this may seriously affect the validity of the conclusions. To estimate the contribution of errors in gene prediction to differences in DA of predicted proteins, we have used the high quality manually curated UniProtKB/Swiss-Prot database as a reference. For genome-scale analysis of domain architectures of predicted proteins we focused on RefSeq, EnsEMBL and NCBI's GNOMON predicted sequences of Metazoan species with completely sequenced genomes. Comparison of the DA of UniProtKB/Swiss-Prot sequences of worm, fly, zebrafish, frog, chick, mouse, rat and orangutan with those of human Swiss-Prot entries have identified relatively few cases where orthologs had different DA, although the percentage with different DA increased with evolutionary distance. In contrast with this, comparison of the DA of human, orangutan, rat, mouse, chicken, frog, zebrafish, worm and fly RefSeq, EnsEMBL and NCBI's GNOMON predicted protein sequences with those of the corresponding/orthologous human Swiss-Prot entries identified a significantly higher proportion of domain architecture differences than in the case of the comparison of Swiss-Prot entries. Analysis of RefSeq, EnsEMBL and NCBI's GNOMON predicted protein sequences with DAs different from those of their Swiss-Prot orthologs confirmed that the higher rate of domain architecture differences is due to errors in gene prediction, the majority of which could be corrected with our FixPred protocol. We have also demonstrated that contamination of databases with incomplete, abnormal or mispredicted sequences introduces a bias in DA differences in as much as it increases the proportion of terminal over internal DA differences. Here we have shown that in the case of RefSeq, EnsEMBL and NCBI's GNOMON predicted protein sequences of Metazoan species, the contribution of gene prediction errors to domain architecture differences of orthologs is comparable to or greater than those due to true gene rearrangements. We have also demonstrated that domain architecture comparison may serve as a useful tool for the quality control of gene predictions and may thus guide the correction of sequence errors. Our findings caution that earlier genome-scale studies based on comparison of predicted (frequently mispredicted) protein sequences may have led to some erroneous conclusions about the evolution of novel domain architectures of multidomain proteins. A reassessment of the DA evolution of orthologous and paralogous proteins is presented in an accompanying paper [1].