DOG 1.0: illustrator of protein domain structures

Structural disorder provides increased adaptability for vesicle trafficking pathways

Vesicle trafficking systems play essential roles in the communication between the organelles of eukaryotic cells and also between cells and their environment. Endocytosis and the late secretory route are mediated by clathrin-coated vesicles, while the COat Protein I and II (COPI and COPII) routes stand for the bidirectional traffic between the ER and the Golgi apparatus. Despite similar fundamental organizations, the molecular machinery, functions, and evolutionary characteristics of the three systems are very different. In this work, we compiled the basic functional protein groups of the three main routes for human and yeast and analyzed them from the structural disorder perspective. We found similar overall disorder content in yeast and human proteins, confirming the well-conserved nature of these systems. Most functional groups contain highly disordered proteins, supporting the general importance of structural disorder in these routes, although some of them seem to heavily rely on disorder, while others do not. Interestingly, the clathrin system is significantly more disordered (∼23%) than the other two, COPI (∼9%) and COPII (∼8%). We show that this structural phenomenon enhances the inherent plasticity and increased evolutionary adaptability of the clathrin system, which distinguishes it from the other two routes. Since multi-functionality (moonlighting) is indicative of both plasticity and adaptability, we studied its prevalence in vesicle trafficking proteins and correlated it with structural disorder. Clathrin adaptors have the highest capability for moonlighting while also comprising the most highly disordered members. The ability to acquire tissue specific functions was also used to approach adaptability: clathrin route genes have the most tissue specific exons encoding for protein segments enriched in structural disorder and interaction sites. Overall, our results confirm the general importance of structural disorder in vesicle trafficking and suggest major roles for this structural property in shaping the differences of evolutionary adaptability in the three routes.

Vesicle trafficking systems are fundamental among cellular transport mechanisms; various cargo molecules are transported via different coated vesicles to their specific destinations in every eukaryotic cell. Clathrin-coated vesicles mediate endocytosis and the late secretory route, while the COat Protein I and II (COPI and COPII) vesicle trafficking routes are responsible for the bidirectional traffic between the ER and the Golgi apparatus. Despite similar basic principles, regulatory mechanisms and structural features of the three systems, their molecular machinery, functions, and evolutionary characteristics vastly differ. We investigated and compared these three routes and their basic functional protein groups from the structural disorder point of view, since disordered protein regions could provide a broad variety of functional and evolutionary advantages for them. We found that structurally disordered protein segments are most abundant in the clathrin system, which might explain the observed inherent plasticity, increased adaptability and exceptional robustness of this route. We support our hypothesis by two analyses on protein multi-functionality and tissue specificity, both being indicative of evolutionary adaptability. Clathrin pathway proteins stand out in both measures, with their disordered regions being largely responsible for their outstanding capabilities.

UHRF2 mRNA expression is low in malignant glioma but silencing inhibits the growth of U251 glioma cells in vitro

UHRF2 is a member of the ubiquitin plant homeo domain RING finger family, which has been proven to be frequently up-regulated in colorectal cancer cells and play a role as an oncogene in breast cancer cells. However, the role of UHRF2 in glioma cells remains unclear. In this study, we performed real-time quantitative PCR on 32 pathologically confirmed glioma samples (grade I, 4 cases; grade II, 11 cases; grade III, 10 cases; and grade IV, 7 cases; according to the 2007 WHO classification system) and four glioma cell lines (A172, U251, U373, and U87). The expression of UHRF2 mRNA was significantly lower in the grade III and grade IV groups compared with the noncancerous brain tissue group, whereas its expression was high in A172, U251, and U373 glioma cell lines. An in vitro assay was performed to investigate the functions of UHRF2. Using a lentivirus-based RNA interference (RNAi) approach, we down-regulated UHRF2 expression in the U251 glioma cell line. This down- regulation led to the inhibition of cell proliferation, an increase in cell apoptosis, and a change of cell cycle distribution, in which S stage cells decreased and G2/M stage cells increased. Our results suggest that UHRF2 may be closely related to tumorigenesis and the development of gliomas.

Interaction of gamma-herpesvirus genome maintenance proteins with cellular chromatin

The capacity of gamma-herpesviruses to establish lifelong infections is dependent on the expression of genome maintenance proteins (GMPs) that tether the viral episomes to cellular chromatin and allow their persistence in latently infected proliferating cells. Here we have characterized the chromatin interaction of GMPs encoded by viruses belonging to the genera Lymphocryptovirus (LCV) and Rhadinovirus (RHV). We found that, in addition to a similar diffuse nuclear localization and comparable detergent resistant interaction with chromatin in transfected cells, all GMPs shared the capacity to promote the decondensation of heterochromatin in the A03-1 reporter cell line. They differed, however, in their mobility measured by fluorescence recovery after photobleaching (FRAP), and in the capacity to recruit accessory molecules required for the chromatin remodeling function. While the AT-hook containing GMPs of LCVs were highly mobile, a great variability was observed among GMPs encoded by RHV, ranging from virtually immobile to significantly reduced mobility compared to LCV GMPs. Only the RHV GMPs recruited the bromo- and extra terminal domain (BET) proteins BRD2 and BRD4 to the site of chromatin remodeling. These findings suggest that differences in the mode of interaction with cellular chromatin may underlie different strategies adopted by these viruses for reprogramming of the host cells during latency.

Analysis of a deep transcriptome from the mantle tissue of Patella vulgata Linnaeus (Mollusca: Gastropoda: Patellidae) reveals candidate biomineralising genes

The gastropod Patella vulgata is abundant on rocky shores in Northern Europe and a significant grazer of intertidal algae. Here we report the application of Illumina sequencing to develop a transcriptome from the adult mantle tissue of P. vulgata. We obtained 47,237,104 paired-end reads of 51 bp, trialled de novo assembly methods and settled on the additive multiple K method followed by redundancy removal as resulting in the most comprehensive assembly. This yielded 29,489 contigs of at least 500 bp in length. We then used three methods to search for candidate genes relevant to biomineralisation: searches via BLAST and Hidden Markov Models for homologues of biomineralising genes from other molluscs, searches for predicted proteins containing tandem repeats and searches for secreted proteins that lacked a transmembrane domain. From the results of these searches we selected 15 contigs for verification by RT-PCR, of which 14 were successfully amplified and cloned. These included homologues of Pif-177/BSMP, Perlustrin, SPARC, AP24, Follistatin-like and Carbonic anhydrase, as well as three containing extensive G-X-Y repeats as found in nacrein. We selected two for further verification by in situ hybridisation, demonstrating expression in the larval shell field. We conclude that de novo assembly of Illumina data offers a cheap and rapid route to a predicted transcriptome that can be used as a resource for further biological study.

Low-resolution structure of the soluble domain GPAA1 (yGPAA170-247) of the glycosylphosphatidylinositol transamidase subunit GPAA1 from Saccharomyces cerevisiae

The GPI (glycosylphosphatidylinositol) transamidase complex catalyses the attachment of GPI anchors to eukaryotic proteins in the lumen of ER (endoplasmic reticulum). The Saccharomyces cerevisiae GPI transamidase complex consists of the subunits yPIG-K (Gpi8p), yPIG-S (Gpi17p), yPIG-T (Gpi16p), yPIG-U (CDC91/GAB1) and yGPAA1. We present the production of the two recombinant proteins yGPAA1^70–247 and yGPAA1^70–339 of the luminal domain of S. cerevisiae GPAA1, covering the amino acids 70–247 and 70–339 respectively. The secondary structural content of the stable and monodisperse yGPAA1^70–247 has been determined to be 28% α-helix and 27% β-sheet. SAXS (small-angle X-ray scattering) data showed that yGPAA1^70–247 has an R_g (radius of gyration) of 2.72±0.025 nm and D_max (maximum dimension) of 9.14 nm. These data enabled the determination of the two domain low-resolution solution structure of yGPAA1^70–247. The large elliptical shape of yGPAA1^70–247 is connected via a short stalk to the smaller hook-like domain of 0.8 nm in length and 3.5 nm in width. The topological arrangement of yGPAA1^70–247 will be discussed together with the recently determined low-resolution structures of yPIG-K^24–337 and yPIG-S^38–467 from S. cerevisiae in the GPI transamidase complex.

Novel interaction of Rab13 and Rab8 with endospanins

Rab GTPases regulate vesicular traffic in eukaryotic cells by cycling between the active GTP-bound and inactive GDP-bound states. Their functions are modulated by the diverse selection of effector proteins that bind to specific Rabs in their activated state. We previously described the expression of Rab13 in bone cells. To search for novel Rab13 interaction partners, we screened a newborn rat bone marrow cDNA library for Rab13 effectors with a bacterial two-hybrid system. We found that Rab13 binds to the C-terminus of Endospanin-2, a small transmembrane protein. In addition to Rab13 also Rab8 bound to Endospanin-2, while no binding of Rab7, Rab10, Rab11 or Rab32 was observed. Rab13 and Rab8 also interacted with Endospanin-1, a close homolog of Endospanin-2. Rab13 and Endospanin-2 colocalised in perinuclear vesicular structures in Cos1 cells suggesting direct binding also in vivo. Endospanin-2 is implicated in the regulation of the cell surface growth hormone receptor (GHR), but the inhibition of Rab13 expression did not affect GHR cell surface expression. This suggests that the Rab13–Endospanin-2 interaction may have functions other than GHR regulation. In conclusion, we have identified a novel interaction for Rab13 and Rab8 with Endospanin-2 and Endospanin-1. The role of this interaction in cell physiology, however, remains to be elucidated.

Multiplex targeted sequencing identifies recurrently mutated genes in autism spectrum disorders

Exome sequencing studies of autism spectrum disorders (ASDs) have identified many de novo mutations but few recurrently disrupted genes. We therefore developed a modified molecular inversion probe method enabling ultra-low-cost candidate gene resequencing in very large cohorts. To demonstrate the power of this approach, we captured and sequenced 44 candidate genes in 2446 ASD probands. We discovered 27 de novo events in 16 genes, 59% of which are predicted to truncate proteins or disrupt splicing. We estimate that recurrent disruptive mutations in six genes-CHD8, DYRK1A, GRIN2B, TBR1, PTEN, and TBL1XR1-may contribute to 1% of sporadic ASDs. Our data support associations between specific genes and reciprocal subphenotypes (CHD8-macrocephaly and DYRK1A-microcephaly) and replicate the importance of a β-catenin-chromatin-remodeling network to ASD etiology.

G9a, a multipotent regulator of gene expression

Lysine methylation of histone and non-histone substrates by the methyltransferase G9a is mostly associated with transcriptional repression. Recent studies, however, have highlighted its role as an activator of gene expression through mechanisms that are independent of its methyltransferase activity. Here we review the growing repertoire of molecular mechanisms and substrates through which G9a regulates gene expression. We also discuss emerging evidence for its wide-ranging functions in development, pluripotency, cellular differentiation and cell cycle regulation that underscore the complexity of its functions. The deregulated expression of G9a in cancers and other human pathologies suggests that it may be a viable therapeutic target in various diseases.

Medulloblastoma exome sequencing uncovers subtype-specific somatic mutations

Medulloblastomas are the most common malignant brain tumors in children¹. Identifying and understanding the genetic events that drive these tumors is critical for the development of more effective diagnostic, prognostic and therapeutic strategies. Recently, our group and others described distinct molecular subtypes of medulloblastoma based on transcriptional and copy number profiles^2–5. Here, we utilized whole exome hybrid capture and deep sequencing to identify somatic mutations across the coding regions of 92 primary medulloblastoma/normal pairs. Overall, medulloblastomas exhibit low mutation rates consistent with other pediatric tumors, with a median of 0.35 non-silent mutations per megabase. We identified twelve genes mutated at statistically significant frequencies, including previously known mutated genes in medulloblastoma such as CTNNB1, PTCH1, MLL2, SMARCA4 and TP53. Recurrent somatic mutations were identified in an RNA helicase gene, DDX3X, often concurrent with CTNNB1 mutations, and in the nuclear co-repressor (N-CoR) complex genes GPS2, BCOR, and LDB1, novel findings in medulloblastoma. We show that mutant DDX3X potentiates transactivation of a TCF promoter and enhances cell viability in combination with mutant but not wild type beta-catenin. Together, our study reveals the alteration of Wnt, Hedgehog, histone methyltransferase and now N-CoR pathways across medulloblastomas and within specific subtypes of this disease, and nominates the RNA helicase DDX3X as a component of pathogenic beta-catenin signaling in medulloblastoma.

Specification, annotation, visualization and simulation of a large rule-based model for ERBB receptor signaling

BACKGROUND

Mathematical/computational models are needed to understand cell signaling networks, which are complex. Signaling proteins contain multiple functional components and multiple sites of post-translational modification. The multiplicity of components and sites of modification ensures that interactions among signaling proteins have the potential to generate myriad protein complexes and post-translational modification states. As a result, the number of chemical species that can be populated in a cell signaling network, and hence the number of equations in an ordinary differential equation model required to capture the dynamics of these species, is prohibitively large. To overcome this problem, the rule-based modeling approach has been developed for representing interactions within signaling networks efficiently and compactly through coarse-graining of the chemical kinetics of molecular interactions.

RESULTS

Here, we provide a demonstration that the rule-based modeling approach can be used to specify and simulate a large model for ERBB receptor signaling that accounts for site-specific details of protein-protein interactions. The model is considered large because it corresponds to a reaction network containing more reactions than can be practically enumerated. The model encompasses activation of ERK and Akt, and it can be simulated using a network-free simulator, such as NFsim, to generate time courses of phosphorylation for 55 individual serine, threonine, and tyrosine residues. The model is annotated and visualized in the form of an extended contact map.

CONCLUSIONS

With the development of software that implements novel computational methods for calculating the dynamics of large-scale rule-based representations of cellular signaling networks, it is now possible to build and analyze models that include a significant fraction of the protein interactions that comprise a signaling network, with incorporation of the site-specific details of the interactions. Modeling at this level of detail is important for understanding cellular signaling.

Histone acetyltransferases in rice (Oryza sativa L.): phylogenetic analysis, subcellular localization and expression

BACKGROUND

Histone acetyltransferases (HATs) play an important role in eukaryotic transcription. Eight HATs identified in rice (OsHATs) can be organized into four families, namely the CBP (OsHAC701, OsHAC703, and OsHAC704), TAFII250 (OsHAF701), GNAT (OsHAG702, OsHAG703, and OsHAG704), and MYST (OsHAM701) families. The biological functions of HATs in rice remain unknown, so a comprehensive protein sequence analysis of the HAT families was conducted to investigate their potential functions. In addition, the subcellular localization and expression patterns of the eight OsHATs were analyzed.

RESULTS

On the basis of a phylogenetic and domain analysis, monocotyledonous CBP family proteins can be subdivided into two groups, namely Group I and Group II. Similarly, dicotyledonous CBP family proteins can be divided into two groups, namely Group A and Group B. High similarities of protein sequences, conserved domains and three-dimensional models were identified among OsHATs and their homologs in Arabidopsis thaliana and maize. Subcellular localization predictions indicated that all OsHATs might localize in both the nucleus and cytosol. Transient expression in Arabidopsis protoplasts confirmed the nuclear and cytosolic localization of OsHAC701, OsHAG702, and OsHAG704. Real-time quantitative polymerase chain reaction analysis demonstrated that the eight OsHATs were expressed in all tissues examined with significant differences in transcript abundance, and their expression was modulated by abscisic acid and salicylic acid as well as abiotic factors such as salt, cold, and heat stresses.

CONCLUSIONS

Both monocotyledonous and dicotyledonous CBP family proteins can be divided into two distinct groups, which suggest the possibility of functional diversification. The high similarities of protein sequences, conserved domains and three-dimensional models among OsHATs and their homologs in Arabidopsis and maize suggested that OsHATs have multiple functions. OsHAC701, OsHAG702, and OsHAG704 were localized in both the nucleus and cytosol in transient expression analyses with Arabidopsis protoplasts. OsHATs were expressed constitutively in rice, and their expression was regulated by exogenous hormones and abiotic stresses, which suggested that OsHATs may play important roles in plant defense responses.

FERM-containing protein FRMD5 is a p120-catenin interacting protein that regulates tumor progression

FERM family proteins have been known to play an important role in tumor progression. FERM-domain containing protein 5 (FRMD5), a novel putative cytoskeletal protein, is an unknown function protein. Here, we reported that FRMD5 localized at the cell adherens junction and formed a molecular complex with p120-catenin through its C-terminal region. Functionally, we found that knockdown of endogenous FRMD5 promotes lung cancer cell migration and invasion in vitro as well as tumor growth in vivo, suggesting a tumor suppressive effect. These findings indicated that FRMD5 may play a role in p120-catenin-based cell-cell contact and is involved in the regulation of tumor progression.

Rho guanine nucleotide exchange factor is an NFL mRNA destabilizing factor that forms cytoplasmic inclusions in amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is an adult-onset progressive disorder of unknown etiology characterized by the selective degeneration of motor neurons. Recent evidence supports the hypothesis that alterations in RNA metabolism in motor neurons can explain the development of protein inclusions, including neurofilamentous aggregates, observed in this pathology. In mice, p190RhoGEF, a guanine nucleotide exchange factor, is involved in neurofilament protein aggregation in an RNA-triggered transgenic model of motor neuron disease. Here, we observed that rho guanine nucleotide exchange factor (RGNEF), the human homologue of p190RhoGEF, binds low molecular weight neurofilament mRNA and affects its stability via 3' untranslated region destabilization. We observed that the overexpression of RGNEF in a stable cell line significantly decreased the level of low molecular weight neurofilament protein. Furthermore, we observed RGNEF cytoplasmic inclusions in ALS spinal motor neurons that colocalized with ubiquitin, p62/sequestosome-1, and TAR (trans-active regulatory) DNA-binding protein 43 (TDP-43). Our results provide further evidence that RNA metabolism pathways are integral to ALS pathology. This is also the first described link between ALS and an RNA binding protein with aggregate formation that is also a central cell signaling pathway molecule.

Hax-1 is rapidly degraded by the proteasome dependent on its PEST sequence

BACKGROUND

HS-1-associated protein X-1 (Hax-1), is a multifunctional protein that has sequence homology to Bcl-2 family members. HAX-1 knockout animals reveal that it plays an essential protective role in the central nervous system against various stresses. Homozygous mutations in the HAX-1 gene are associated with autosomal recessive forms of severe congenital neutropenia along with neurological symptoms. The protein level of Hax-1 has been shown to be regulated by cellular protease cleavage or by transcriptional suppression upon stimulation.

RESULTS

Here, we report a novel post-translational mechanism for regulation of Hax-1 levels in mammalian cells. We identified that PEST sequence, a sequence rich in proline, glutamic acid, serine and threonine, is responsible for its poly-ubiquitination and rapid degradation. Hax-1 is conjugated by K48-linked ubiquitin chains and undergoes a fast turnover by the proteasome system. A deletion mutant of Hax-1 that lacks the PEST sequence is more resistant to the proteasomal degradation and exerts more protective effects against apoptotic stimuli than wild type Hax-1.

CONCLUSION

Our data indicate that Hax-1 is a short-lived protein and that its PEST sequence dependent fast degradation by the proteasome may contribute to the rapid cellular responses upon different stimulations.

Comparative Analysis of SWIRM Domain-Containing Proteins in Plants

Chromatin-remodeling complexes affect gene expression by using the energy of ATP hydrolysis to locally disrupt or alter the association of histones with DNA. SWIRM (Swi3p, Rsc8p, and Moira) domain is an alpha-helical domain of about 85 residues in chromosomal proteins. SWIRM domain-containing proteins make up large multisubunit complexes by interacting with other chromatin modification factors and may have an important function in plants. However, little is known about SWIRM domain-containing proteins in plants. In this study, 67 SWIRM domain-containing proteins from 6 plant species were identified and analyzed. Plant SWIRM domain proteins can be divided into three distinct types: Swi-type, LSD1-type, and Ada2-type. Generally, the SWIRM domain forms a helix-turn-helix motif commonly found in DNA-binding proteins. The genes encoding SWIRM domain proteins in Oryza sativa are widely expressed, especially in pistils. In addition, OsCHB701 and OsHDMA701 were downregulated by cold stress, whereas OsHDMA701 and OsHDMA702 were significantly induced by heat stress. These observations indicate that SWIRM domain proteins may play an essential role in plant development and plant responses to environmental stress.

Characterization of barley Prp1 gene and its expression during seed development and under abiotic stress

The pre-mRNA processing (Prp1) gene encodes a spliceosomal protein. It was firstly identified in fission yeast and plays a regular role during spliceosome activation and cell cycle. Plant Prp1 genes have only been identified from rice, Sorghum and Arabidopsis thaliana. In this study, we reported the identification and isolation of a novel Prp1 gene from barley, and further explored its expressional pattern by using real-time quantitative RTPCR, promoter prediction and analysis of microarray data. The putative barley Prp1 protein has a similar primary structure features to those of other known Prp1 protein in this family. The results of amino acid comparison indicated that Prp1 protein of barley and other plant species has a highly conserved 30 termnal region while their 50 sequences greatly varied. The results of expressional analysis revealed that the expression level of barley Prp1 gene is always stable in different vegetative tissues, except it is up-regulated at the mid- and late stages of seed development or under the condition of cold stress. This kind of expressional pattern for barley Prp1 is also supported by our results of comparison of microarray data from barley, rice and Arabidopsis. For the molecular mechanism of its expressional pattern, we conclude that the expression of Prp1 gene may be up-regulated by the increase of pre-mRNAs and not be constitutive or ubiquitous.

The evolving proteome of a complex extracellular matrix, the Oikopleura house

Extracellular matrices regulate biological processes at the level of cells, tissues, and in some cases, entire multicellular organisms. The subphylum Urochordata exemplifies the latter case, where animals are partially or completely enclosed in “houses” or “tunics”. Despite this common strategy, we show that the house proteome of the appendicularian, Oikopleura, has very little in common with the proteome of the sister class, ascidian, Ciona. Of 80 identified house proteins (oikosins), ∼half lack domain modules or similarity to known proteins, suggesting de novo appearance in appendicularians. Gene duplication has been important in generating almost 1/3 of the current oikosin complement, with serial duplications up to 8 paralogs in one family. Expression pattern analyses revealed that individual oikosins are produced from specific fields of cells within the secretory epithelium, but in some cases, migrate up to at least 20 cell diameters in extracellular space to combine in defined house structures. Interestingly, peroxidasin and secretory phospholipase A₂ domains, implicated in innate immune defence are secreted from the anlage associated with the food-concentrating filter, suggesting that this extra-organismal structure may play, in part, such a role in Oikopleura. We also show that sulfation of proteoglycans is required for the hydration and inflation of pre-house rudiments into functional houses. Though correct proportioning in the production of oikosins would seem important in repetitive assembly of the complex house structure, the genomic organization of oikosin loci appears incompatible with common enhancers or locus control regions exerting such a coordinate regulatory role. Thus, though all tunicates employ extracellular matrices based on a cellulose scaffold as a defining feature of the subphylum, they have evolved radically different protein compositions associated with this common underlying structural theme.

Evolutionary histories of expanded peptidase families in Schistosoma mansoni

Schistosoma mansoni is one of the three main causative agents of human schistosomiasis, a major health problem with a vast socio-economic impact. Recent advances in the proteomic analysis of schistosomes have revealed that peptidases are the main virulence factors involved in the pathogenesis of this disease. In this context, evolutionary studies can be applied to identify peptidase families that have been expanded in genomes over time in response to different selection pressures. Using a phylogenomic approach, we searched for expanded endopeptidase families in the S. mansoni predicted proteome with the aim of contributing to the knowledge of such enzymes as potential therapeutic targets. We found three endopeptidase families that comprise leishmanolysins (metallopeptidase M8 family), cercarial elastases (serine peptidase S1 family) and cathepsin D proteins (aspartic peptidase A1 family). Our results suggest that the Schistosoma members of these families originated from successive gene duplication events in the parasite lineage after its diversification from other metazoans. Overall, critical residues are conserved among the duplicated genes/proteins. Furthermore, each protein family displays a distinct evolutionary history. Altogether, this work provides an evolutionary view of three S. mansoni peptidase families, which allows for a deeper understanding of the genomic complexity and lineage-specific adaptations potentially related to the parasitic lifestyle.

RUNX1 mutations in cytogenetically normal acute myeloid leukemia are associated with a poor prognosis and up-regulation of lymphoid genes

BACKGROUND

The RUNX1 (AML1) gene is a frequent mutational target in myelodysplastic syndromes and acute myeloid leukemia. Previous studies suggested that RUNX1 mutations may have pathological and prognostic implications.

DESIGN AND METHODS

We screened 93 patients with cytogenetically normal acute myeloid leukemia for RUNX1 mutations by capillary sequencing of genomic DNA. Mutation status was then correlated with clinical data and gene expression profiles.

RESULTS

We found that 15 out of 93 (16.1%) patients with cytogenetically normal acute myeloid leukemia had RUNX1 mutations. Seventy-three patients were enrolled in the AMLCG-99 trial and carried ten RUNX1 mutations (13.7%). Among these 73 patients RUNX1 mutations were significantly associated with older age, male sex, absence of NPM1 mutations and presence of MLL-partial tandem duplications. Moreover, RUNX1-mutated patients had a lower complete remission rate (30% versus 73% P=0.01), lower relapse-free survival rate (3-year relapse-free survival 0% versus 30.4%; P=0.002) and lower overall survival rate (3-year overall survival 0% versus 34.4%; P<0.001) than patients with wild-type RUNX1. RUNX1 mutations remained associated with shorter overall survival in a multivariate model including age and the European Leukemia Net acute myeloid leukemia genetic classification as covariates. Patients with RUNX1 mutations showed a unique gene expression pattern with differential expression of 85 genes. The most prominently up-regulated genes in patients with RUNX1-mutated cytogenetically normal acute myeloid leukemia include lymphoid regulators such as HOP homeobox (HOPX), deoxynucleotidyltransferase (DNTT, terminal) and B-cell linker (BLNK), indicating lineage infidelity.

CONCLUSIONS

Our findings firmly establish that RUNX1 mutations are a marker of poor prognosis and provide insights into the pathogenesis of RUNX1 mutation-positive acute myeloid leukemia.

Interactions of an Arabidopsis RanBPM homologue with LisH-CTLH domain proteins revealed high conservation of CTLH complexes in eukaryotes

BACKGROUND

RanBPM (Ran-binding protein in the microtubule-organizing centre) was originally reported as a centrosome-associated protein in human cells. However, RanBPM protein containing highly conserved SPRY, LisH, CTLH and CRA domains is currently considered as a scaffolding protein with multiple cellular functions. A plant homologue of RanBPM has not yet been characterized.

RESULTS

Based on sequence similarity, we identified a homologue of the human RanBPM in Arabidopsis thaliana. AtRanBPM protein has highly conserved SPRY, LisH, CTLH and CRA domains. Cell fractionation showed that endogenous AtRanBPM or expressed GFP-AtRanBPM are mainly cytoplasmic proteins with only a minor portion detectable in microsomal fractions. AtRanBPM was identified predominantly in the form of soluble cytoplasmic complexes ~230-500 kDa in size. Immunopurification of AtRanBPM followed by mass spectrometric analysis identified proteins containing LisH and CRA domains; LisH, CRA, RING-U-box domains and a transducin/WD40 repeats in a complex with AtRanBPM. Homologues of identified proteins are known to be components of the C-terminal to the LisH motif (CTLH) complexes in humans and budding yeast. Microscopic analysis of GFP-AtRanBPM in vivo and immunofluorescence localization of endogenous AtRanBPM protein in cultured cells and seedlings of Arabidopsis showed mainly cytoplasmic and nuclear localization. Absence of colocalization with γ-tubulin was consistent with the biochemical data and suggests another than a centrosomal role of the AtRanBPM protein.

CONCLUSION

We showed that as yet uncharacterized Arabidopsis RanBPM protein physically interacts with LisH-CTLH domain-containing proteins. The newly identified high molecular weight cytoplasmic protein complexes of AtRanBPM showed homology with CTLH types of complexes described in mammals and budding yeast. Although the exact functions of the CTLH complexes in scaffolding of protein degradation, in protein interactions and in signalling from the periphery to the cell centre are not yet fully understood, structural conservation of the complexes across eukaryotes suggests their important biological role.

Striatal-enriched protein-tyrosine phosphatase (STEP) regulates Pyk2 kinase activity

Proline-rich tyrosine kinase 2 (Pyk2) is a member of the focal adhesion kinase family and is highly expressed in brain and hematopoietic cells. Pyk2 plays diverse functions in cells, including the regulation of cell adhesion, migration, and cytoskeletal reorganization. In the brain, it is involved in the induction of long term potentiation through regulation of N-methyl-d-aspartate receptor trafficking. This occurs through the phosphorylation and activation of Src family tyrosine kinase members, such as Fyn, that phosphorylate GluN2B at Tyr(1472). Phosphorylation at this site leads to exocytosis of GluN1-GluN2B receptors to synaptic membranes. Pyk2 activity is modulated by phosphorylation at several critical tyrosine sites, including Tyr(402). In this study, we report that Pyk2 is a substrate of striatal-enriched protein-tyrosine phosphatase (STEP). STEP binds to and dephosphorylates Pyk2 at Tyr(402). STEP KO mice showed enhanced phosphorylation of Pyk2 at Tyr(402) and of the Pyk2 substrates paxillin and ASAP1. Functional studies indicated that STEP opposes Pyk2 activation after KCl depolarization of cortical slices and blocks Pyk2 translocation to postsynaptic densities, a key step required for Pyk2 activation and function. This is the first study to identify Pyk2 as a substrate for STEP.

Phosphorylation of the Drosophila melanogaster RNA-binding protein HOW by MAPK/ERK enhances its dimerization and activity

Drosophila melanogaster Held Out Wings (HOW) is a conserved RNA-binding protein (RBP) belonging to the STAR family, whose closest mammalian ortholog Quaking (QKI) has been implicated in embryonic development and nervous system myelination. The HOW RBP modulates a variety of developmental processes by controlling mRNA levels and the splicing profile of multiple key regulatory genes; however, mechanisms regulating its activity in tissues have yet to be elucidated. Here, we link receptor tyrosine kinase (RTK) signaling to the regulation of QKI subfamily of STAR proteins, by showing that HOW undergoes phosphorylation by MAPK/ERK. Importantly, we show that this modification facilitates HOW dimerization and potentiates its ability to bind RNA and regulate its levels. Employing an antibody that specifically recognizes phosphorylated HOW, we show that HOW is phosphorylated in embryonic muscles and heart cardioblasts in vivo, thus documenting for the first time Serine/Threonine (Ser/Thr) phosphorylation of a STAR protein in the context of an intact organism. We also identify the sallimus/D-titin (sls) gene as a novel muscle target of HOW-mediated negative regulation and further show that this regulation is phosphorylation-dependent, underscoring the physiological relevance of this modification. Importantly, we demonstrate that HOW Thr phosphorylation is reduced following muscle-specific knock down of Drosophila MAPK rolled and that, correspondingly, Sls is elevated in these muscles, similarly to the HOW RNAi effect. Taken together, our results provide a coherent mechanism of differential HOW activation; MAPK/ERK-dependent phosphorylation of HOW promotes the formation of HOW dimers and thus enhances its activity in controlling mRNA levels of key muscle-specific genes. Hence, our findings bridge between MAPK/ERK signaling and RNA regulation in developing muscles.

Bisulfite sequencing reveals that Aspergillus flavus holds a hollow in DNA methylation

Aspergillus flavus first gained scientific attention for its production of aflatoxin. The underlying regulation of aflatoxin biosynthesis has been serving as a theoretical model for biosynthesis of other microbial secondary metabolites. Nevertheless, for several decades, the DNA methylation status, one of the important epigenomic modifications involved in gene regulation, in A. flavus remains to be controversial. Here, we applied bisulfite sequencing in conjunction with a biological replicate strategy to investigate the DNA methylation profiling of A. flavus genome. Both the bisulfite sequencing data and the methylome comparisons with other fungi confirm that the DNA methylation level of this fungus is negligible. Further investigation into the DNA methyltransferase of Aspergillus uncovers its close relationship with RID-like enzymes as well as its divergence with the methyltransferase of species with validated DNA methylation. The lack of repeat contents of the A. flavus' genome and the high RIP-index of the small amount of remanent repeat potentially support our speculation that DNA methylation may be absent in A. flavus or that it may possess de novo DNA methylation which occurs very transiently during the obscure sexual stage of this fungal species. This work contributes to our understanding on the DNA methylation status of A. flavus, as well as reinforces our views on the DNA methylation in fungal species. In addition, our strategy of applying bisulfite sequencing to DNA methylation detection in species with low DNA methylation may serve as a reference for later scientific investigations in other hypomethylated species.

Structural basis for the substrate specificity of a novel β-N-acetylhexosaminidase StrH protein from Streptococcus pneumoniae R6

The β-N-acetylhexosaminidase (EC 3.2.1.52) from glycoside hydrolase family 20 (GH20) catalyzes the hydrolysis of the β-N-acetylglucosamine (NAG) group from the nonreducing end of various glycoconjugates. The putative surface-exposed N-acetylhexosaminidase StrH/Spr0057 from Streptococcus pneumoniae R6 was proved to contribute to the virulence by removal of β(1,2)-linked NAG on host defense molecules following the cleavage of sialic acid and galactose by neuraminidase and β-galactosidase, respectively. StrH is the only reported GH20 enzyme that contains a tandem repeat of two 53% sequence-identical catalytic domains (designated as GH20-1 and GH20-2, respectively). Here, we present the 2.1 Å crystal structure of the N-terminal domain of StrH (residues Glu-175 to Lys-642) complexed with NAG. It adopts an overall structure similar to other GH20 enzymes: a (β/α)(8) TIM barrel with the active site residing at the center of the β-barrel convex side. The kinetic investigation using 4-nitrophenyl N-acetyl-β-d-glucosaminide as the substrate demonstrated that GH20-1 had an enzymatic activity (k(cat)/K(m)) of one-fourth compared with GH20-2. The lower activity of GH20-1 could be attributed to the substitution of active site Cys-469 of GH20-1 to the counterpart Tyr-903 of GH20-2. A complex model of NAGβ(1,2)Man at the active site of GH20-1 combined with activity assays of the corresponding site-directed mutants characterized two key residues Trp-443 and Tyr-482 at subsite +1 of GH20-1 (Trp-876 and Tyr-914 of GH20-2) that might determine the β(1,2) substrate specificity. Taken together, these findings shed light on the mechanism of catalytic specificity toward the β(1,2)-linked β-N-acetylglucosides.

Intrinsic disorder of the extracellular matrix

The extracellular matrix is very well organized at the supramolecular and tissue levels and little is known on the potential role of intrinsic disorder in promoting its organization. We predicted the amount of disorder and identified disordered regions in the human extracellular proteome with established computational tools. The extracellular proteome is significantly enriched in proteins comprising more than 50% of disorder compared to the complete human proteome. The enrichment is mostly due to long disordered regions containing at least 100 consecutive disordered residues. The amount of intrinsic disorder is heterogeneous in the extracellular protein families, with the most disordered being collagens and the small integrin-binding ligand N-linked glycoproteins. Although most domains found in extracellular proteins are structured, the fibronectin III domains contain a variable amount of disordered residues (up to 92%). Binding sites for heparin and integrins are found in disordered sequences of extracellular proteins. Intrinsic disorder is evenly distributed in hubs and ends in the interaction network of extracellular proteins with their extracellular partners. In contrast, extracellular hubs are significantly enriched in disorder in the network of extracellular proteins with their extracellular, membrane and intracellular partners. Disorder could thus provide the structural plasticity required for the hubs to interact with membrane and intracellular proteins. Organization and assembly of the extracellular matrix, development of mineralized tissues and cell-matrix adhesion are the biological processes overrepresented in the most disordered extracellular proteins. Extracellular disorder is associated with binding to growth factors, glycosaminoglycans and integrins at the molecular level.

Endophytic life strategies decoded by genome and transcriptome analyses of the mutualistic root symbiont Piriformospora indica

Recent sequencing projects have provided deep insight into fungal lifestyle-associated genomic adaptations. Here we report on the 25 Mb genome of the mutualistic root symbiont Piriformospora indica (Sebacinales, Basidiomycota) and provide a global characterization of fungal transcriptional responses associated with the colonization of living and dead barley roots. Extensive comparative analysis of the P. indica genome with other Basidiomycota and Ascomycota fungi that have diverse lifestyle strategies identified features typically associated with both, biotrophism and saprotrophism. The tightly controlled expression of the lifestyle-associated gene sets during the onset of the symbiosis, revealed by microarray analysis, argues for a biphasic root colonization strategy of P. indica. This is supported by a cytological study that shows an early biotrophic growth followed by a cell death-associated phase. About 10% of the fungal genes induced during the biotrophic colonization encoded putative small secreted proteins (SSP), including several lectin-like proteins and members of a P. indica-specific gene family (DELD) with a conserved novel seven-amino acids motif at the C-terminus. Similar to effectors found in other filamentous organisms, the occurrence of the DELDs correlated with the presence of transposable elements in gene-poor repeat-rich regions of the genome. This is the first in depth genomic study describing a mutualistic symbiont with a biphasic lifestyle. Our findings provide a significant advance in understanding development of biotrophic plant symbionts and suggest a series of incremental shifts along the continuum from saprotrophy towards biotrophy in the evolution of mycorrhizal association from decomposer fungi.

Guidelines for visualizing and annotating rule-based models

Rule-based modeling provides a means to represent cell signaling systems in a way that captures site-specific details of molecular interactions. For rule-based models to be more widely understood and (re)used, conventions for model visualization and annotation are needed. We have developed the concepts of an extended contact map and a model guide for illustrating and annotating rule-based models. An extended contact map represents the scope of a model by providing an illustration of each molecule, molecular component, direct physical interaction, post-translational modification, and enzyme-substrate relationship considered in a model. A map can also illustrate allosteric effects, structural relationships among molecular components, and compartmental locations of molecules. A model guide associates elements of a contact map with annotation and elements of an underlying model, which may be fully or partially specified. A guide can also serve to document the biological knowledge upon which a model is based. We provide examples of a map and guide for a published rule-based model that characterizes early events in IgE receptor (FcεRI) signaling. We also provide examples of how to visualize a variety of processes that are common in cell signaling systems but not considered in the example model, such as ubiquitination. An extended contact map and an associated guide can document knowledge of a cell signaling system in a form that is visual as well as executable. As a tool for model annotation, a map and guide can communicate the content of a model clearly and with precision, even for large models.

Visualization and biochemical analyses of the emerging mammalian 14-3-3-phosphoproteome

Hundreds of candidate 14-3-3-binding (phospho)proteins have been reported in publications that describe one interaction at a time, as well as high-throughput 14-3-3-affinity and mass spectrometry-based studies. Here, we transcribed these data into a common format, deposited the collated data from low-throughput studies in MINT (http://mint.bio.uniroma2.it/mint), and compared the low- and high-throughput data in VisANT graphs that are easy to analyze and extend. Exploring the graphs prompted questions about technical and biological specificity, which were addressed experimentally, resulting in identification of phosphorylated 14-3-3-binding sites in the mitochondrial import sequence of the iron-sulfur cluster assembly enzyme (ISCU), cytoplasmic domains of the mitochondrial fission factor (MFF), and endoplasmic reticulum-tethered receptor expression-enhancing protein 4 (REEP4), RNA regulator SMAUG2, and cytoskeletal regulatory proteins, namely debrin-like protein (DBNL) and kinesin light chain (KLC) isoforms. Therefore, 14-3-3s undergo physiological interactions with proteins that are destined for diverse subcellular locations. Graphing and validating interactions underpins efforts to use 14-3-3-phosphoproteomics to identify mechanisms and biomarkers for signaling pathways in health and disease.

Functional Diversity of the Schistosoma mansoni Tyrosine Kinases

Schistosoma mansoni, one of the causative agents of schistosomiasis, has a complex life cycle infecting over 200 million people worldwide. Such a successful and prolific parasite life cycle has been shown to be dependent on the adaptive interaction between the parasite and hosts. Tyrosine kinases (TKs) play a key role in signaling pathways as demonstrated by a large body of experimental work in eukaryotes. Furthermore, comparative genomics have allowed the identification of TK homologs and provided insights into the functional role of TKs in several biological systems. Finally, TK structural biology has provided a rational basis for obtaining selective inhibitors directed to the treatment of human diseases. This paper covers the important aspects of the phospho-tyrosine signaling network in S. mansoni, Caenorhabditis elegans, and humans, the main process of functional diversification of TKs, that is, protein-domain shuffling, and also discusses TKs as targets for the development of new anti-schistosome drugs.

Molecular evolution of the transferrin family and associated receptors

BACKGROUND

In vertebrates, serum transferrins are essential iron transporters that have bind and release Fe(III) in response to receptor binding and changes in pH. Some family members such as lactoferrin and melanotransferrin can also bind iron while others have lost this ability and have gained other functions, e.g., inhibitor of carbonic anhydrase (mammals), saxiphilin (frogs) and otolith matrix protein 1 (fish).

SCOPE OF REVIEW

This article provides an overview of the known transferrin family members and their associated receptors and interacting partners.

MAJOR CONCLUSIONS

The number of transferrin genes has proliferated as a result of multiple duplication events, and the resulting paralogs have developed a wide array of new functions. Some homologs in the most primitive metazoan groups resemble both serum and melanotransferrins, but the major yolk proteins show considerable divergence from the rest of the family. Among the transferrin receptors, the lack of TFR2 in birds and reptiles, and the lack of any TFR homologs among the insects draw attention to the differences in iron transport and regulation in those groups.

GENERAL SIGNIFICANCE

The transferrin family members are important because of their clinical significance, interesting biochemical properties, and evolutionary history. More work is needed to better understand the functions and evolution of the non-vertebrate family members. This article is part of a Special Issue entitled Molecular Mechanisms of Iron Transport and Disorders.

Eukaryotic protein kinases (ePKs) of the helminth parasite Schistosoma mansoni

BACKGROUND

Schistosomiasis remains an important parasitic disease and a major economic problem in many countries. The Schistosoma mansoni genome and predicted proteome sequences were recently published providing the opportunity to identify new drug candidates. Eukaryotic protein kinases (ePKs) play a central role in mediating signal transduction through complex networks and are considered druggable targets from the medical and chemical viewpoints. Our work aimed at analyzing the S. mansoni predicted proteome in order to identify and classify all ePKs of this parasite through combined computational approaches. Functional annotation was performed mainly to yield insights into the parasite signaling processes relevant to its complex lifestyle and to select some ePKs as potential drug targets.

RESULTS

We have identified 252 ePKs, which corresponds to 1.9% of the S. mansoni predicted proteome, through sequence similarity searches using HMMs (Hidden Markov Models). Amino acid sequences corresponding to the conserved catalytic domain of ePKs were aligned by MAFFT and further used in distance-based phylogenetic analysis as implemented in PHYLIP. Our analysis also included the ePK homologs from six other eukaryotes. The results show that S. mansoni has proteins in all ePK groups. Most of them are clearly clustered with known ePKs in other eukaryotes according to the phylogenetic analysis. None of the ePKs are exclusively found in S. mansoni or belong to an expanded family in this parasite. Only 16 S. mansoni ePKs were experimentally studied, 12 proteins are predicted to be catalytically inactive and approximately 2% of the parasite ePKs remain unclassified. Some proteins were mentioned as good target for drug development since they have a predicted essential function for the parasite.

CONCLUSIONS

Our approach has improved the functional annotation of 40% of S. mansoni ePKs through combined similarity and phylogenetic-based approaches. As we continue this work, we will highlight the biochemical and physiological adaptations of S. mansoni in response to diverse environments during the parasite development, vector interaction, and host infection.

GPS-CCD: a novel computational program for the prediction of calpain cleavage sites

As one of the most essential post-translational modifications (PTMs) of proteins, proteolysis, especially calpain-mediated cleavage, plays an important role in many biological processes, including cell death/apoptosis, cytoskeletal remodeling, and the cell cycle. Experimental identification of calpain targets with bona fide cleavage sites is fundamental for dissecting the molecular mechanisms and biological roles of calpain cleavage. In contrast to time-consuming and labor-intensive experimental approaches, computational prediction of calpain cleavage sites might more cheaply and readily provide useful information for further experimental investigation. In this work, we constructed a novel software package of GPS-CCD (Calpain Cleavage Detector) for the prediction of calpain cleavage sites, with an accuracy of 89.98%, sensitivity of 60.87% and specificity of 90.07%. With this software, we annotated potential calpain cleavage sites for hundreds of calpain substrates, for which the exact cleavage sites had not been previously determined. In this regard, GPS-CCD 1.0 is considered to be a useful tool for experimentalists. The online service and local packages of GPS-CCD 1.0 were implemented in JAVA and are freely available at: http://ccd.biocuckoo.org/.

Toward a complete in silico, multi-layered embryonic stem cell regulatory network

Recent efforts in systematically profiling embryonic stem (ES) cells have yielded a wealth of high-throughput data. Complementarily, emerging databases and computational tools facilitate ES cell studies and further pave the way toward the in silico reconstruction of regulatory networks encompassing multiple molecular layers. Here, we briefly survey databases, algorithms, and software tools used to organize and analyze high-throughput experimental data collected to study mammalian cellular systems with a focus on ES cells. The vision of using heterogeneous data to reconstruct a complete multi-layered ES cell regulatory network is discussed. This review also provides an accompanying manually extracted dataset of different types of regulatory interactions from low-throughput experimental ES cell studies available at http://amp.pharm.mssm.edu/iscmid/literature.

RNF13: an emerging RING finger ubiquitin ligase important in cell proliferation

Protein ubiquitination mediated by ubiquitin ligases plays a very important role in a wide spectrum of biological processes including development and disease pathogenesis. RING finger protein 13 (RNF13) is a recently identified ubiquitin ligase which contains an N-terminal protease-associated domain and a C-terminal RING finger domain separated by a transmembrane region. RNF13 is an evolutionarily conserved protein. Most interestingly, RNF13 expression is developmentally regulated during myogenesis and is upregulated in various human tumors. These data suggest that RNF13, acting as an ubiquitin ligase, might have profound biological functions during development and disease. This minireview summarizes recent work on RNF13 functions related to cell proliferation, differentiation and cancer development.

Cellular roles of neuronal calcium sensor-1 and calcium/calmodulin-dependent kinases in fungi

The neuronal calcium sensor-1 (NCS-1) possesses a consensus signal for N-terminal myristoylation and four EF-hand Ca(2+)-binding sites, and mediates the effects of cytosolic Ca(2+). Minute changes in free intracellular Ca(2+) are quickly transformed into changes in the activity of several kinases including calcium/calmodulin-dependent protein kinases (Ca(2+)/CaMKs) that are involved in regulating many eukaryotic cell functions. However, our current knowledge of NCS-1 and Ca(2+)/CaMKs comes mostly from studies of the mammalian enzymes. Thus far very few fungal homologues of NCS-1 and Ca(2+)/CaMKs have been characterized and little is known about their cellular roles. In this minireview, we describe the known sequences, interactions with target proteins and cellular roles of NCS-1 and Ca(2+)/CaMKs in fungi.

GPS 2.1: enhanced prediction of kinase-specific phosphorylation sites with an algorithm of motif length selection

As the most important post-translational modification of proteins, phosphorylation plays essential roles in all aspects of biological processes. Besides experimental approaches, computational prediction of phosphorylated proteins with their kinase-specific phosphorylation sites has also emerged as a popular strategy, for its low-cost, fast-speed and convenience. In this work, we developed a kinase-specific phosphorylation sites predictor of GPS 2.1 (Group-based Prediction System), with a novel but simple approach of motif length selection (MLS). By this approach, the robustness of the prediction system was greatly improved. All algorithms in GPS old versions were also reserved and integrated in GPS 2.1. The online service and local packages of GPS 2.1 were implemented in JAVA 1.5 (J2SE 5.0) and freely available for academic researches at: http://gps.biocuckoo.org.

More than 1,001 problems with protein domain databases: transmembrane regions, signal peptides and the issue of sequence homology

Large-scale genome sequencing gained general importance for life science because functional annotation of otherwise experimentally uncharacterized sequences is made possible by the theory of biomolecular sequence homology. Historically, the paradigm of similarity of protein sequences implying common structure, function and ancestry was generalized based on studies of globular domains. Having the same fold imposes strict conditions over the packing in the hydrophobic core requiring similarity of hydrophobic patterns. The implications of sequence similarity among non-globular protein segments have not been studied to the same extent; nevertheless, homology considerations are silently extended for them. This appears especially detrimental in the case of transmembrane helices (TMs) and signal peptides (SPs) where sequence similarity is necessarily a consequence of physical requirements rather than common ancestry. Thus, matching of SPs/TMs creates the illusion of matching hydrophobic cores. Therefore, inclusion of SPs/TMs into domain models can give rise to wrong annotations. More than 1001 domains among the 10,340 models of Pfam release 23 and 18 domains of SMART version 6 (out of 809) contain SP/TM regions. As expected, fragment-mode HMM searches generate promiscuous hits limited to solely the SP/TM part among clearly unrelated proteins. More worryingly, we show explicit examples that the scores of clearly false-positive hits, even in global-mode searches, can be elevated into the significance range just by matching the hydrophobic runs. In the PIR iProClass database v3.74 using conservative criteria, we find that at least between 2.1% and 13.6% of its annotated Pfam hits appear unjustified for a set of validated domain models. Thus, false-positive domain hits enforced by SP/TM regions can lead to dramatic annotation errors where the hit has nothing in common with the problematic domain model except the SP/TM region itself. We suggest a workflow of flagging problematic hits arising from SP/TM-containing models for critical reconsideration by annotation users.

Sequence homology is a fundamental principle of biology. It implies common phylogenetic ancestry of genes and, subsequently, similarity of their protein products with regard to amino acid sequence, three-dimensional structure and molecular and cellular function. Originally an esoteric concept, homology with the proxy of sequence similarity is used to justify the transfer of functional annotation from well-studied protein examples to new sequences. Yet, functional annotation via sequence similarity seems to have hit a plateau in recent years since relentless annotation transfer led to error propagation across sequence databases; thus, leading experimental follow-up work astray. It must be emphasized that the trinity of sequence, 3D structural and functional similarity has only been proven for globular segments of proteins. For non-globular regions, similarity of sequence is not necessarily a result of divergent evolution from a common ancestor but the consequence of amino acid sequence bias. In our investigation, we found that protein domain databases contain many domain models with transmembrane regions and signal peptides, non-globular segments of proteins having hydrophobic bias. Many proteins have inherited completely wrong function assignments from these domain models. We fear that future function predictions will turn out futile if this issue is not immediately addressed.

PhosSNP for systematic analysis of genetic polymorphisms that influence protein phosphorylation

We are entering the era of personalized genomics as breakthroughs in sequencing technology have made it possible to sequence or genotype an individual person in an efficient and accurate manner. Preliminary results from HapMap and other similar projects have revealed the existence of tremendous genetic variations among world populations and among individuals. It is important to delineate the functional implication of such variations, i.e. whether they affect the stability and biochemical properties of proteins. It is also generally believed that the genetic variation is the main cause for different susceptibility to certain diseases or different response to therapeutic treatments. Understanding genetic variation in the context of human diseases thus holds the promise for "personalized medicine." In this work, we carried out a genome-wide analysis of single nucleotide polymorphisms (SNPs) that could potentially influence protein phosphorylation characteristics in human. Here, we defined a phosphorylation-related SNP (phosSNP) as a non-synonymous SNP (nsSNP) that affects the protein phosphorylation status. Using an in-house developed kinase-specific phosphorylation site predictor (GPS 2.0), we computationally detected that approximately 70% of the reported nsSNPs are potential phosSNPs. More interestingly, approximately 74.6% of these potential phosSNPs might also induce changes in protein kinase types in adjacent phosphorylation sites rather than creating or removing phosphorylation sites directly. Taken together, we proposed that a large proportion of the nsSNPs might affect protein phosphorylation characteristics and play important roles in rewiring biological pathways. Finally, all phosSNPs were integrated into the PhosSNP 1.0 database, which was implemented in JAVA 1.5 (J2SE 5.0). The PhosSNP 1.0 database is freely available for academic researchers.

New insights into the active site structure and catalytic mechanism of tyrosinase and its related proteins

Tyrosinases are widely distributed in nature. They are copper-containing oxidases belonging to the type 3 copper protein family, together with catechol oxidases and haemocyanins. Tyrosinases are essential enzymes in melanin biosynthesis and therefore responsible for pigmentation of skin and hair in mammals, where two more enzymes, the tyrosinase-related proteins (Tyrps), participate in the pathway. The structure and catalytic mechanism of mammalian tyrosinases have been extensively studied but they are not completely understood because of the lack of information on the tertiary structure. The availability of crystallographic data of one plant catechol oxidase and one bacterial tyrosinase has improved the model of the three-dimensional structure of the active site of the enzyme. Furthermore, sequence comparison of tyrosinase and the Tyrps reveals that the three orthologue proteins share many key structural features, because of their common origin from an ancestral gene, although the specific residues responsible for their different catalytic capabilities have not been identified yet. This review summarizes our current knowledge of tyrosinase and Tyrps structure and function and describes the catalytic mechanism of tyrosinase and Dct/Tyrp2, which are better characterized.

A fratricidal mechanism is responsible for eDNA release and contributes to biofilm development of Enterococcus faecalis

Extracellular DNA (eDNA), a by-product of cell lysis, was recently established as a critical structural component of the Enterococcus faecalis biofilm matrix. Here, we describe fratricide as the governing principle behind gelatinase (GelE)-mediated cell death and eDNA release. GFP reporter assays confirmed that GBAP (gelatinase biosynthesis-activating pheromone) quorum non-responders (GelE-SprE-) were a minority subpopulation of prey cells susceptible to the targeted fratricidal action of the quorum responsive predatorial majority (GelE+SprE+). The killing action is dependent on GelE, and the GelE producer population is protected from self-destruction by the co-production of SprE as an immunity protein. Targeted gene inactivation and protein interaction studies demonstrate that extracellular proteases execute their characteristic effects following downstream interactions with the primary autolysin, AtlA. Finally, we address a mechanism by which GelE and SprE may modify the cell wall affinity of proteolytically processed AtlA resulting in either a pro- or anti-lytic outcome.