Computational differentiation of N-terminal signal peptides and transmembrane helices

Signal peptides and transmembrane helices both contain a stretch of hydrophobic amino acids. This common feature makes it difficult for signal peptide and transmembrane helix predictors to correctly assign identity to stretches of hydrophobic residues near the N-terminal methionine of a protein sequence. The inability to reliably distinguish between N-terminal transmembrane helix and signal peptide is an error with serious consequences for the prediction of protein secretory status or transmembrane topology. In this study, we report a new method for differentiating protein N-terminal signal peptides and transmembrane helices. Based on the sequence features extracted from hydrophobic regions (amino acid frequency, hydrophobicity, and the start position), we set up discriminant functions and examined them on non-redundant datasets with jackknife tests. This method can incorporate other signal peptide prediction methods and achieve higher prediction accuracy. For Gram-negative bacterial proteins, 95.7% of N-terminal signal peptides and transmembrane helices can be correctly predicted (coefficient 0.90). Given a sensitivity of 90%, transmembrane helices can be identified from signal peptides with a precision of 99% (coefficient 0.92). For eukaryotic proteins, 94.2% of N-terminal signal peptides and transmembrane helices can be correctly predicted with coefficient 0.83. Given a sensitivity of 90%, transmembrane helices can be identified from signal peptides with a precision of 87% (coefficient 0.85). The method can be used to complement current transmembrane protein prediction and signal peptide prediction methods to improve their prediction accuracies. (C) 2003 Elsevier Inc. All rights reserved.

Can power laws help us understand gene and proteome information?

Proteins are biochemical entities consisting of one or more blocks typically folded in a 3D pattern. Each block (a polypeptide) is a single linear sequence of amino acids that are biochemically bonded together. The amino acid sequence in a protein is defined by the sequence of a gene or several genes encoded in the DNA-based genetic code. This genetic code typically uses twenty amino acids, but in certain organisms the genetic code can also include two other amino acids. After linking the amino acids during protein synthesis, each amino acid becomes a residue in a protein, which is then chemically modified, ultimately changing and defining the protein function. In this study, the authors analyze the amino acid sequence using alignment-free methods, aiming to identify structural patterns in sets of proteins and in the proteome, without any other previous assumptions. The paper starts by analyzing amino acid sequence data by means of histograms using fixed length amino acid words (tuples). After creating the initial relative frequency histograms, they are transformed and processed in order to generate quantitative results for information extraction and graphical visualization. Selected samples from two reference datasets are used, and results reveal that the proposed method is able to generate relevant outputs in accordance with current scientific knowledge in domains like protein sequence/proteome analysis.

Cultivo de Eucalyptus urograndis em atmosfera enriquecida com CO2: mudanças no proteoma cloroplastidial

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Multiple adenosine 3',5'-cyclic [corrected] monophosphate response element DNA-binding proteins generated by gene diversification and alternative exon splicing.

The protein sequence deduced from the open reading frame of a human placental cDNA encoding a cAMP-responsive enhancer (CRE)-binding protein (CREB-327) has structural features characteristic of several other transcriptional transactivator proteins including jun, fos, C/EBP, myc, and CRE-BP1. Results of Southwestern analysis of nuclear extracts from several different cell lines show that there are multiple CRE-binding proteins, which vary in size in cell lines derived from different tissues and animal species. To examine the molecular diversity of CREB-327 and related proteins at the nucleic acid level, we used labeled cDNAs from human placenta that encode two different CRE-binding proteins (CREB-327 and CRE-BP1) to probe Northern and Southern blots. Both probes hybridized to multiple fragments on Southern blots of genomic DNA from various species. Alternatively, when a human placental c-jun probe was hybridized to the same blot, a single fragment was detected in most cases, consistent with the intronless nature of the human c-jun gene. The CREB-327 probe hybridized to multiple mRNAs, derived from human placenta, ranging in size from 2-9 kilobases. In contrast, the CRE-BP1 probe identified a single 4-kilobase mRNA. Sequence analyses of several overlapping human genomic cosmid clones containing CREB-327 sequences in conjunction with polymerase chain reaction indicates that the CREB-327/341 cDNAs are composed of at least eight or nine exons, and analyses of human placental cDNAs provide direct evidence for at least one alternatively spliced exon. Analyses of mouse/hamster-human hybridoma DNAs by Southern blotting and polymerase chain reaction localizes the CREB-327/341 gene to human chromosome 2. The results indicate that there is a dichotomy of CREB-like proteins, those that are related by overall structure and DNA-binding specificity as well as those that are related by close similarities of primary sequences.

Molecular and antigenic characterisation of ribosomal phosphoprotein P0 from Babesia bovis

Babesia bovis is a tick-borne pathogen that remains an important constraint for the development of cattle industries in tropical and subtropical regions of the world. Effective control can be achieved by vaccination with live attenuated phenotypes of the parasite. However, these phenotypes have a number of drawbacks, which justifies the search for new, more efficient immunogens based mainly on recombinant protein technology. In the present paper, ribosomal phosphoprotein P0 from a Brazilian isolate of B. bovis was produced and evaluated with regard to conservation and antigenicity. The protein sequence displayed high conservation between different Brazilian isolates of B. bovis and several Apicomplexa parasites such as Theileria, Neospora and Toxoplasma. IgG from cattle experimentally and naturally infected with B. bovisas well as IgG1 and IgG2 from naturally infected cattle reacted with the recombinant protein. IgG from cattle experimentally infected with Babesia bigemina cross-reacted with B. bovis recombinant P0. These characteristics suggest that P0 is a potential antigen for recombinant vaccine preparations against bovine babesiosis.

Molecular analysis of non-melanoma skin cancer

SUMMARY:Cylindroma, trichoepithelioma and spiradenoma are benign tumors of hair follicle. They are caused by mutations and loss of heterozygosity in the CYLD gene. CYLD is a ubiquitously expressed, but the tumors are restricted to skin, suggesting that the tumorigenesis is influenced by skin-specific regulators and probably by mutations in other genes. The objectives of the thesis were to analyze the molecular mechanisms leading to the aforementioned tumors. In the first project, we have identified five new mutations in CYLD gene in tive families affected with different combinations of these skin appendage tumors. F our of these mutations caused the introduction of a premature stop codon in CYLD protein sequence, but one was a missense mutation changing aspartic acid 681 into glycine (D68lG), in patients exhibiting multiple trichoepitheliomas. CYLD is a deubiquitinase which can downregulate NF-κB and INK pathways through the deubiquitination of TRAF2, for example. We showed that the CYLD-D681G mutant was unable to remove polyubiquitin chains from TRAF2. We also proved that CYLD-D68lG could not inhibit TRAP 2- or TNFα- mediated NF-κB or INK activations in 293T cells. These results underlined the importance of the D68l residue for the enzymatic activity of CYLD. TRAP-interacting protein (TRIP), which is a E3-Ubiquitin ligase, is a partner of CYLD. In the second project of the thesis, we studied the function of TRIP in the epidermis. We found that TRIP was a nucleolar protein in cultured human primary keratinocytes (HEK) and HeLa cells, and was detected in the midbody of HeLa cells. Moreover, TRIP expression was shown to be downregulated through a PKC-dependent mechanism before induction of keratinocyte differentiation. We also proved that TRIP was upregulated in basal cell carcinomas. Furthermore, TRIP was found to be important for keratinocyte survival and proliferation through the regulation of the Gl/S transition. Our results suggest that TRIP may be involved in keratinocyte tumorigenesis.RÉSUMÉ :Les cylindromes, trichoépithéliomes et spiradénomes sont des tumeurs bénignes du follicule pileux causées par des mutations et une perte d'hétérozygotie du gène CYLD. CYLD est ubiquitaire mais les tumeurs sont limitées à la peau, suggérant que la tumorigénèse est influencée par des protéines spécifiques de la peau et par des mutations dans d'autres gènes. Les objectifs de la thèse étaient d'2malyser les mécanismes moléculaires aboutissant à la formation de ces tumeurs. Dans le premier projet, cinq nouvelles mutations du gène CYLD ont été identifiées chez cinq familles présentant différentes combinaisons des tumeurs citées ci- dessus. Quatre de ces mutations causaient I' introduction d'un codon stop prématuré dans la séquence protéique, mais une était une mutation «misser1se» changeant l'aspartate 681 en résidu glycine (D68lG) chez des patients présentant des trichoépithéliomes multiples. CYLD est une déubiquitinase qui inhibe les voies de signalisation de NF-κB et JNK, en déubiquitinant notamment TRAF2. Nous avons montré que la protéine mutante CYLD- D68lG ne pouvait pas cliver la chaîne de poly-ubiquitines liée à TRAF2. CYLD-D68lG était aussi incapable d'inhiber l'activation de NF-κB ou de JNK induite par TRAF2 ou TNF-o dans les cellules 293T. Ces résultats ont donc souligné l'impo1tance du résidu D68l pour l'activité de CYLD. «TRAF-interacting protein (TRIP)», qui est une «E3-ubiquitin-ligase», est un partenaire de CYLD. Dans le second proj et de la thèse, nous avons étudié la fonction de TRIP dans l'épidenne. Nous avons montrépque TRIP était nucléolaire dans les cellules HeLa et les kératinocytes primaires humains en culture et était détectée dans le «midbody» des cellules HeLa. Nous avons prouvé que l'ARNm de TRIP était diminué avant l'induction de la différentiation des kératinocytes, par un mécanisme dépendent de la protéine kinase C, tandis qu'il était augmenté dans les carcinomes baso-cellulaires. Nous avons aussi montré que TRIP influençait la prolifération et la survie des kératinocytes en régulant la transition G1/S, Nos résultats suggèrent que TRIP est peut-être impliquée dans la tumorigénèse des kératinocytes. 7

SelenoDB 1.0 : a database of selenoprotein genes, proteins and SECIS elements

Selenoproteins are a diverse group of proteinsusually misidentified and misannotated in sequencedatabases. The presence of an in-frame UGA (stop)codon in the coding sequence of selenoproteingenes precludes their identification and correctannotation. The in-frame UGA codons are recodedto cotranslationally incorporate selenocysteine,a rare selenium-containing amino acid. The developmentof ad hoc experimental and, more recently,computational approaches have allowed the efficientidentification and characterization of theselenoproteomes of a growing number of species.Today, dozens of selenoprotein families have beendescribed and more are being discovered in recentlysequenced species, but the correct genomic annotationis not available for the majority of thesegenes. SelenoDB is a long-term project that aims toprovide, through the collaborative effort of experimentaland computational researchers, automaticand manually curated annotations of selenoproteingenes, proteins and SECIS elements. Version 1.0 ofthe database includes an initial set of eukaryoticgenomic annotations, with special emphasis on thehuman selenoproteome, for immediate inspectionby selenium researchers or incorporation into moregeneral databases. SelenoDB is freely available athttp://www.selenodb.org.

Expression and one-step purification of Plasmodium proteins in dictyostelium.

Nearly full-length Circumsporozoite protein (CSP) from Plasmodium falciparum, the C-terminal fragments from both P. falciparm and P. yoelii CSP and a fragment comprising 351 amino acids of P.vivax MSPI were expressed in the slime mold Dictyostelium discoideum. Discoidin-tag expression vectors allowed both high yields of these proteins and their purification by a nearly single-step procedure. We exploited the galactose binding activity of Discoidin Ia to separate the fusion proteins by affinity chromatography on Sepharose-4B columns. Inclusion of a thrombin recognition site allowed cleavage of the Discoidin-tag from the fusion protein. Partial secretion of the protein was obtained via an ER independent pathway, whereas routing the recombinant proteins to the ER resulted in glycosylation and retention. Yields of proteins ranged from 0.08 to 3 mg l(-1) depending on the protein sequence and the purification conditions. The recognition of purified MSPI by sera from P. vivax malaria patients was used to confirm the native conformation of the protein expressed in Dictyostelium. The simple purification procedure described here, based on Sepharose-4B, should facilitate the expression and the large-scale purification of various Plasmodium polypeptides.

Type I signal peptidase from Leishmania is a target of the immune response in human cutaneous and visceral leishmaniasis.

The gene encoding type I signal peptidase (Lmjsp) has been cloned from Leishmania major. Lmjsp encodes a protein of 180 amino residues with a predicted molecular mass of 20.5 kDa. Comparison of the protein sequence with those of known type I signal peptidases indicates homology in five conserved domains A-E which are known to be important, or essential, for catalytic activity. Southern blot hybridisation analysis indicates that there is a single copy of the Lmjsp gene. A recombinant SPase protein and a synthetic peptide of the L. major signal peptidase were used to examine the presence of specific antibodies in sera from either recovered or active individuals of both cutaneous and visceral leishmaniasis. This evaluation demonstrated that sera from cutaneous and visceral forms of leishmaniasis are highly reactive to both the recombinant and synthetic signal peptidase antigens. Therefore, the Leishmania signal peptidase, albeit localised intracellularly, is a significant target of the Leishmania specific immune response and highlights its potential use for serodiagnosis of cutaneous and visceral leishmaniasis.

Mutations leading to X-linked hypohidrotic ectodermal dysplasia affect three major functional domains in the tumor necrosis factor family member ectodysplasin-A.

Mutations in the epithelial morphogen ectodysplasin-A (EDA), a member of the tumor necrosis factor (TNF) family, are responsible for the human disorder X-linked hypohidrotic ectodermal dysplasia (XLHED) characterized by impaired development of hair, eccrine sweat glands, and teeth. EDA-A1 and EDA-A2 are two splice variants of EDA, which bind distinct EDA-A1 and X-linked EDA-A2 receptors. We identified a series of novel EDA mutations in families with XLHED, allowing the identification of the following three functionally important regions in EDA: a C-terminal TNF homology domain, a collagen domain, and a furin protease recognition sequence. Mutations in the TNF homology domain impair binding of both splice variants to their receptors. Mutations in the collagen domain can inhibit multimerization of the TNF homology region, whereas those in the consensus furin recognition sequence prevent proteolytic cleavage of EDA. Finally, a mutation affecting an intron splice donor site is predicted to eliminate specifically the EDA-A1 but not the EDA-A2 splice variant. Thus a proteolytically processed, oligomeric form of EDA-A1 is required in vivo for proper morphogenesis.

Analyses of six homologous proteins of Protochlamydia amoebophila UWE25 encoded by large GC-rich genes (lgr): a model of evolution and concatenation of leucine-rich repeats.

BACKGROUND: Along the chromosome of the obligate intracellular bacteria Protochlamydia amoebophila UWE25, we recently described a genomic island Pam100G. It contains a tra unit likely involved in conjugative DNA transfer and lgrE, a 5.6-kb gene similar to five others of P. amoebophila: lgrA to lgrD, lgrF. We describe here the structure, regulation and evolution of these proteins termed LGRs since encoded by "Large G+C-Rich" genes. RESULTS: No homologs to the whole protein sequence of LGRs were found in other organisms. Phylogenetic analyses suggest that serial duplications producing the six LGRs occurred relatively recently and nucleotide usage analyses show that lgrB, lgrE and lgrF were relocated on the chromosome. The C-terminal part of LGRs is homologous to Leucine-Rich Repeats domains (LRRs). Defined by a cumulative alignment score, the 5 to 18 concatenated octacosapeptidic (28-meric) LRRs of LGRs present all a predicted alpha-helix conformation. Their closest homologs are the 28-residue RI-like LRRs of mammalian NODs and the 24-meres of some Ralstonia and Legionella proteins. Interestingly, lgrE, which is present on Pam100G like the tra operon, exhibits Pfam domains related to DNA metabolism. CONCLUSION: Comparison of the LRRs, enable us to propose a parsimonious evolutionary scenario of these domains driven by adjacent concatenations of LRRs. Our model established on bacterial LRRs can be challenged in eucaryotic proteins carrying less conserved LRRs, such as NOD proteins and Toll-like receptors.

DAI/ZBP1 recruits RIP1 and RIP3 through RIP homotypic interaction motifs to activate NF-kappaB.

Detection of viral nucleic acids is central to antiviral immunity. Recently, DAI/ZBP1 (DNA-dependent activator of IRFs/Z-DNA binding protein 1) was identified as a cytoplasmic DNA sensor and shown to activate the interferon regulatory factor (IRF) and nuclear factor-kappa B (NF-kappaB) transcription factors, leading to type-I interferon production. DAI-induced IRF activation depends on TANK-binding kinase 1 (TBK1), whereas signalling pathways and molecular components involved in NF-kappaB activation remain elusive. Here, we report the identification of two receptor-interacting protein (RIP) homotypic interaction motifs (RHIMs) in the DAI protein sequence, and show that these domains relay DAI-induced NF-kappaB signals through the recruitment of the RHIM-containing kinases RIP1 and RIP3. We show that knockdown of not only RIP1, but also RIP3 affects DAI-induced NF-kappaB activation. Importantly, RIP recruitment to DAI is inhibited by the RHIM-containing murine cytomegalovirus (MCMV) protein M45. These findings delineate the DAI signalling pathway to NF-kappaB and suggest a possible new immune modulation strategy of the MCMV.

Physicochemical and biological characterization of single-walled and double-walled carbon nanotubes in biological media

To study the toxicity of nanoparticles under relevant conditions, it is important to reproducibly disperse nanoparticles in biological media in in vitro and in vivo studies. Here, single-walled nanotubes (SWNTs) and double-walled nanotubes (DWNTs) were physicochemically and biologically characterized when dispersed in phosphate-buffered saline (PBS) and bovine serum albumin (BSA). BSA-SWNT/DWNT interaction resulted in a reduction of aggregation and an increase in particle stabilization. Based on the protein sequence coverage and protein binding results, DWNTs exhibited higher protein binding than SWNTs. SWNT and DWNT suspensions in the presence of BSA increased interleukin-6 (IL-6) levels and reduced tumor necrosis factor-alpha (TNF-α) levels in A549 cells as compared to corresponding samples in the absence of BSA. We next determined the effects of SWNTs and DWNTs on pulmonary protein modification using bronchoalveolar lavage fluid (BALF) as a surrogate collected form BALB/c mice. The BALF proteins bound to SWNTs (13 proteins) and DWNTs (11 proteins), suggesting that these proteins were associated with blood coagulation pathways. Lastly, we demonstrated the importance of physicochemical and biological alterations of SWNTs and DWNTs when dispersed in biological media, since protein binding may result in the misinterpretation of in vitro results and the activation of protein-regulated biological responses.

Characterization of alternatively spliced products and tissue-specific isoforms of USP28 and USP25

Background: The ubiquitin-dependent protein degradation pathway is essential for the proteolysis of intracellular proteins and peptides. Deubiquitinating enzymes constitute a complex protein family involved in a multitude of cellular processes. The ubiquitin-specific proteases (UBP) are a group of enzymes whose predicted function is to reverse the ubiquitinating reaction by removing ubiquitin from a large variety of substrates. We have lately reported the characterization of human USP25, a specific-ubiquitin protease gene at 21q11.2, with a specific pattern of expression in murine fetal brains and adult testis. Results: Database homology searches at the DNA and protein levels and cDNA library screenings led to the identification of a new UBP member in the human genome, named USP28, at 11q23. This novel gene showed preferential expression in heart and muscle. Moreover, cDNA, expressed sequence tag and RT-PCR analyses provided evidence for alternatively spliced products and tissue-specific isoforms. Concerning function, USP25 overexpression in Down syndrome fetal brains was shown by real-time PCR. Conclusions: On the basis of the genomic and protein sequence as well as the functional data, USP28 and USP25 establish a new subfamily of deubiquitinating enzymes. Both genes have alternatively spliced exons that could generate protein isoforms with distinct tissue-specific activity. The overexpression of USP25 in Down syndrome fetal brains supports the gene-dosage effects suggested for other UBP members related to aneuploidy syndromes.

ORMDL proteins are a conserved new family of endoplasmic reticulum membrane proteins

Background: Annotations of completely sequenced genomes reveal that nearly half of the genes identified are of unknown function, and that some belong to uncharacterized gene families. To help resolve such issues, information can be obtained from the comparative analysis of homologous genes in model organisms. Results: While characterizing genes from the retinitis pigmentosa locus RP26 at 2q31-q33, we have identified a new gene, ORMDL1, that belongs to a novel gene family comprising three genes in humans (ORMDL1, ORMDL2 and ORMDL3), and homologs in yeast, microsporidia, plants, Drosophila, urochordates and vertebrates. The human genes are expressed ubiquitously in adult and fetal tissues. The Drosophila ORMDL homolog is also expressed throughout embryonic and larval stages, particularly in ectodermally derived tissues. The ORMDL genes encode transmembrane proteins anchored in the endoplasmic reticulum (ER). Double knockout of the two Saccharomyces cerevisiae homologs leads to decreased growth rate and greater sensitivity to tunicamycin and dithiothreitol. Yeast mutants can be rescued by human ORMDL homologs. Conclusions: From protein sequence comparisons we have defined a novel gene family, not previously recognized because of the absence of a characterized functional signature. The sequence conservation of this family from yeast to vertebrates, the maintenance of duplicate copies in different lineages, the ubiquitous pattern of expression in human and Drosophila, the partial functional redundancy of the yeast homologs and phenotypic rescue by the human homologs, strongly support functional conservation. Subcellular localization and the response of yeast mutants to specific agents point to the involvement of ORMDL in protein folding in the ER.