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.

Analysis of Tks4 Knockout Mice Suggests a Role for Tks4 in Adipose Tissue Homeostasis in the Context of Beigeing.

Obesity and adipocyte malfunction are related to and arise as consequences of disturbances in signaling pathways. Tyrosine kinase substrate with four Src homology 3 domains (Tks4) is a scaffold protein that establishes a platform for signaling cascade molecules during podosome formation and epidermal growth factor receptor (EGFR) signaling. Several lines of evidence have also suggested that Tks4 has a role in adipocyte biology; however, its roles in the various types of adipocytes at the cellular level and in transcriptional regulation have not been studied. Therefore, we hypothesized that Tks4 functions as an organizing molecule in signaling networks that regulate adipocyte homeostasis. Our aims were to study the white and brown adipose depots of Tks4 knockout (KO) mice using immunohistology and western blotting and to analyze gene expression changes regulated by the white, brown, and beige adipocyte-related transcription factors via a PCR array. Based on morphological differences in the Tks4-KO adipocytes and increased uncoupling protein 1 (UCP1) expression in the white adipose tissue (WAT) of Tks4-KO mice, we concluded that the beigeing process was more robust in the WAT of Tks4-KO mice compared to the wild-type animals. Furthermore, in the Tks4-KO WAT, the expression profile of peroxisome proliferator-activated receptor gamma (PPARγ)-regulated adipogenesis-related genes was shifted in favor of the appearance of beige-like cells. These results suggest that Tks4 and its downstream signaling partners are novel regulators of adipocyte functions and PPARγ-directed white to beige adipose tissue conversion.

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.

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.

Influence of WFIKKN1 on BMP1-mediated activation of latent myostatin.

The NTR domain of WFIKKN1 protein has been shown to have significant affinity for the prodomain regions of promyostatin and latent myostatin but the biological significance of these interactions remained unclear. In view of its role as a myostatin antagonist, we tested the assumption that WFIKKN1 inhibits the release of myostatin from promyostatin and/or latent myostatin. WFIKKN1 was found to have no effect on processing of promyostatin by furin, the rate of cleavage of latent myostatin by BMP1, however, was significantly enhanced in the presence of WFIKKN1 and this enhancer activity was superstimulated by heparin. Unexpectedly, WFIKKN1 was also cleaved by BMP1 and our studies have shown that the KKN1 fragment generated by BMP1-cleavage of WFIKKN1 contributes most significantly to the observed enhancer activity. Analysis of a pro-TGF-ß -based homology model of homodimeric latent myostatin revealed that the BMP1-cleavage sites are buried and not readily accessible to BMP1. In view of this observation, the most plausible explanation for the BMP1-enhancer activity of the KKN1 fragment is that it shifts a conformational equilibrium of latent myostatin from the closed circular structure of the homodimer to a more open form, making the cleavage sites more accessible to BMP1. On the other hand, the observation that the enhancer activity of KKN1 is superstimulated in the presence of heparin is explained by the fact KKN1, latent myostatin, and BMP1 have affinity for heparin and these interactions with heparin increase the local concentrations of the reactants thereby facilitating the action of BMP1. ENZYMES: Furin: EC 3.4.21.75; BMP1, bone morphogentic protein 1 or procollagen C-endopeptidase: EC 3.4.24.19.

K153R polymorphism in myostatin gene increases the rate of promyostatin activation by furin.

Recent studies demonstrated an association between the K153R polymorphism in the myostatin gene with extreme longevity, lower muscle strength and obesity but the molecular basis of these associations has not been clarified. Here, we show that the K153R mutation significantly increases the rate of proteolysis of promyostatin by furin, but has no effect on the activity of the latent complex or the cleavage of the latent complex by bone morphogenetic protein 1 (BMP-1). The increased rate of activation of K153R mutant promyostatin may explain why this polymorphism is associated with obesity, lower muscle strength and extension of lifespan.

Both LCCL-domains of human CRISPLD2 have high affinity for lipid A.

The LCCL-domain is a recently defined protein module present in diverse extracellular multidomain proteins. Practically nothing is known about the molecular function of these domains; based on functional features of proteins harboring LCCL-domains it has been suggested that these domains might function as lipopolysaccharide-binding domains. Here we show that the two LCCL-domains of human CRISPLD2 protein, a lipopolysaccharide-binding serum protein involved in defense against endotoxin shock, have higher affinity for the lipid A, the toxic moiety of lipopolysaccharides than for ipopolysaccharide. Our observation that the LCCL-domains of CRISPLD2 are specific for the toxic lipid A moiety of the endotoxin suggests that it may block the interaction between endotoxins and the host endotoxin receptors without interfering with the development of antibacterial immunity against the polysaccharide moiety of LPS. We suggest that the anti-inflammatory function of CRISPLD2 protein may account for its role in various pathological and developmental processes.

Latent myostatin has significant activity and this activity is controlled more efficiently by WFIKKN1 than by WFIKKN2.

Myostatin, a negative regulator of skeletal muscle growth, is produced from myostatin precursor by multiple steps of proteolytic processing. After cleavage by a furin-type protease, the propeptide and growth factor domains remain associated, forming a noncovalent complex, the latent myostatin complex. Mature myostatin is liberated from latent myostatin by bone morphogenetic protein 1/tolloid proteases. Here, we show that, in reporter assays, latent myostatin preparations have significant myostatin activity, as the noncovalent complex dissociates at an appreciable rate, and both mature and semilatent myostatin (a complex in which the dimeric growth factor domain interacts with only one molecule of myostatin propeptide) bind to myostatin receptor. The interaction of myostatin receptor with semilatent myostatin is efficiently blocked by WAP, Kazal, immunoglobulin, Kunitz and NTR domain-containing protein 1 or growth and differentiation factor-associated serum protein 2 (WFIKKN1), a large extracellular multidomain protein that binds both mature myostatin and myostatin propeptide [Kondás et al. (2008) J Biol Chem 283, 23677-23684]. Interestingly, the paralogous protein WAP, Kazal, immunoglobulin, Kunitz and NTR domain-containing protein 2 or growth and differentiation factor-associated serum protein 1 (WFIKKN2) was less efficient than WFIKKN1 as an antagonist of the interactions of myostatin receptor with semilatent myostatin. Our studies have shown that this difference is attributable to the fact that only WFIKKN1 has affinity for the propeptide domain, and this interaction increases its potency in suppressing the receptor-binding activity of semilatent myostatin. As the interaction of WFIKKN1 with various forms of myostatin permits tighter control of myostatin activity until myostatin is liberated from latent myostatin by bone morphogenetic protein 1/tolloid proteases, WFIKKN1 may have greater potential as an antimyostatic agent than WFIKKN2.

A novel COCH mutation associated with autosomal dominant nonsyndromic hearing loss disrupts the structural stability of the vWFA2 domain.

Mutations in COCH have been associated with autosomal dominant nonsyndromic hearing loss (DFNA9) and are frequently accompanied by vestibular hypofunction. Here, we report identification of a novel missense mutation, p.F527C, located in the vWFA2 domain in members of a Korean family with late-onset and progressive hearing loss. To assess the molecular characteristics of this cochlin mutant, we constructed both wild-type and mutant cochlin constructs and transfected these into mammalian cell lines. Results of immunocytochemistry analysis demonstrated localization of the cochlin mutant in the endoplasmic reticulum/Golgi complex, whereas western blot analyses of cell lysates revealed that the mutant cochlin tends to form covalent complexes that are retained in the cell. Biochemical analyses of recombinant vWFA2 domain of cochlin carrying the p.F527C mutation revealed that the mutation increases propensity of the protein to form covalent disulfide-bonded dimers and affects the structural stability but not the collagen-affinity of the vWFA2 domain. We suggest that the instability of mutant cochlin is the major driving force for cochlin aggregation in the inner ear in DFNA9 patients carrying the COCH p.F527C mutation.

Biological functions of the WAP domain-containing multidomain proteins WFIKKN1 and WFIKKN2.

WFIKKN1 and WFIKKN2 are two closely related multidomain proteins consisting of a WAP (whey acidic protein)-, a follistatin-, an immunoglobulin-, two Kunitz-type protease inhibitor-domains and an NTR domain (netrin domain). Recent experiments have shown that both WFIKKN1 and WFIKKN2 bind myostatin and GDF11 (growth and differentiation factor 11) with high affinity and are potent antagonists of these growth factors. Structure-function studies on WFIKKN proteins have revealed that their interactions with GDF8 and GDF11 are mediated primarily by the follistatin and NTR domains.

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].

WFIKKN1 and WFIKKN2 bind growth factors TGFß1, BMP2 and BMP4 but do not inhibit their signalling activity.

WFIKKN1 and WFIKKN2 are large extracellular multidomain proteins consisting of a WAP domain, a follistatin domain, an immunoglobulin domain, two Kunitz-type protease inhibitor domains and an NTR domain. Recent experiments have shown that both proteins have high affinity for growth and differentiation factor (GDF)8 and GDF11. Here we study the interaction of WFIKKN proteins with several additional representatives of the transforming growth factor (TGF)ß family using SPR measurements. Analyses of SPR sensorgrams suggested that, in addition to GDF8 and GDF11, both WFIKKN proteins bind TGFß1, bone morphogenetic protein (BMP)2 and BMP4 with relatively high affinity (K(d) ∼ 10(-6) m). To assess the biological significance of these interactions we studied the effect of WFIKKN proteins on the activity of GDF8, GDF11, TGFß1, BMP2 and BMP4 using reporter assays. These studies revealed that WFIKKN1 and WFIKKN2 inhibited the biological activity of GDF8 and GDF11 in the nanomolar range, whereas they did not inhibit the activities of TGFß1, BMP2 and BMP4 even in the micromolar range. Our data indicate that WFIKKN proteins are antagonists of GDF8 and GDF11, but in the case of TGFß1, BMP2 and BMP4 they function as growth factor binding proteins. It is suggested that the physical association of WFIKKN proteins with these growth factors may localize their action and thus help to establish growth factor gradients in the extracellular space.

The second von Willebrand type A domain of cochlin has high affinity for type I, type II and type IV collagens.

Cochlin is colocalized with type II collagen in the extracellular matrix of cochlea and has been suggested to interact with this collagen. Here we show that the second von Willebrand type A domain of cochlin has affinity for type II collagen, as well as type I and type IV collagens whereas the LCCL-domain of cochlin has no affinity for these proteins. The implications of these findings for the mechanism whereby cochlin mutations cause the dominant negative DFNA9-type hearing loss are discussed.

NMR solution structure of the neurotrypsin Kringle domain.

Neurotrypsin is a multidomain protein that serves as a brain-specific serine protease. Here we report the NMR structure of its kringle domain, NT/K. The data analysis was performed with the BACUS (Bayesian analysis of coupled unassigned spins) algorithm. This study presents the first application of BACUS to the structure determination of a 13C unenriched protein for which no prior experimental 3D structure was available. NT/K adopts the kringle fold, consisting of an antiparallel beta-sheet bridged by an overlapping pair of disulfides. The structure reveals the presence of a surface-exposed left-handed polyproline II helix that is closely packed to the core beta-structure. This feature distinguishes NT/K from other members of the kringle fold and points toward a novel functional role for a kringle domain. Functional divergence among kringle domains is discussed on the basis of their surface and electrostatic characteristics.

Both WFIKKN1 and WFIKKN2 have high affinity for growth and differentiation factors 8 and 11.

WFIKKN1 and WFIKKN2 are large extracellular multidomain proteins consisting of a WAP, a follistatin, an immunoglobulin, two Kunitz-type protease inhibitor domains, and an NTR domain. Recent experiments have shown that WFIKKN2 protein binds mature GDF8/myostatin and myostatin propeptide and inhibits the biological activity of myostatin (Hill, J. J., Qiu, Y., Hewick, R. M., and Wolfman, N. M. (2003) Mol. Endocrinol. 17, 1144-1154). Here we show that the paralogue of this protein, WFIKKN1, also binds to both myostatin and myostatin propeptide and that both WFIKKN1 and WFIKKN2 bind GDF11, the growth and differentiation factor most closely related to myostatin, with high affinity. Structure-function studies on WFIKKN1 have revealed that the follistatin domain is primarily responsible for the binding of mature growth factor, whereas the NTR domain contributes most significantly to the interaction with myostatin propeptide. Analysis of the evolutionary histories of WFIKKN1/WFIKKN2 and GDF8/GDF11 proteins indicates that the functional association of an ancestral WFIKKN protein with an ancestor of GDF8/11 may date back to cephalochordates/urochordates. Although duplication of the corresponding genes gave rise to WFIKKN1/WFIKKN2 and GDF8/GDF11 in early vertebrates, the data presented here suggest that there is significant functional overlap of the paralogous proteins.