Keiichi Homma

Binary classification of protein molecules into intrinsically disordered and ordered segments

BackgroundAlthough structural domains in proteins (SDs) are important, half of the regions in the human proteome are currently left with no SD assignments. These unassigned regions consist not only of novel SDs, but also of intrinsically disordered (ID) regions since proteins, especially those in eukaryotes, generally contain a significant fraction of ID regions. As ID regions can be inferred from amino acid sequences, a method that combines SD and ID region assignments can determine the fractions of SDs and ID regions in any proteome.ResultsIn contrast to other available ID prediction programs that merely identify likely ID regions, the DICHOT system we previously developed classifies the entire protein sequence into SDs and ID regions. Application of DICHOT to the human proteome revealed that residue-wise ID regions constitute 35%, SDs with similarity to PDB structures comprise 52%, while SDs with no similarity to PDB structures account for the remaining 13%. The last group consists of novel structural domains, termed cryptic domains, which serve as good targets of structural genomics. The DICHOT method applied to the proteomes of other model organisms indicated that eukaryotes generally have high ID contents, while prokaryotes do not. In human proteins, ID contents differ among subcellular localizations: nuclear proteins had the highest residue-wise ID fraction (47%), while mitochondrial proteins exhibited the lowest (13%). Phosphorylation and O-linked glycosylation sites were found to be located preferentially in ID regions. As O-linked glycans are attached to residues in the extracellular regions of proteins, the modification is likely to protect the ID regions from proteolytic cleavage in the extracellular environment. Alternative splicing events tend to occur more frequently in ID regions. We interpret this as evidence that natural selection is operating at the protein level in alternative splicing.ConclusionsWe classified entire regions of proteins into the two categories, SDs and ID regions and thereby obtained various kinds of complete genome-wide statistics. The results of the present study are important basic information for understanding protein structural architectures and have been made publicly available at http://spock.genes.nig.ac.jp/~genome/DICHOT.

Comparative Genome Analysis of Three Eukaryotic Parasites with Differing Abilities To Transform Leukocytes Reveals Key Mediators of Theileria-Induced Leukocyte Transformation

ABSTRACT We sequenced the genome of Theileria orientalis, a tick-borne apicomplexan protozoan parasite of cattle. The focus of this study was a comparative genome analysis of T. orientalis relative to other highly pathogenic Theileria species, T. parva and T. annulata. T. parva and T. annulata induce transformation of infected cells of lymphocyte or macrophage/monocyte lineages; in contrast, T. orientalis does not induce uncontrolled proliferation of infected leukocytes and multiplies predominantly within infected erythrocytes. While synteny across homologous chromosomes of the three Theileria species was found to be well conserved overall, subtelomeric structures were found to differ substantially, as T. orientalis lacks the large tandemly arrayed subtelomere-encoded variable secreted protein-encoding gene family. Moreover, expansion of particular gene families by gene duplication was found in the genomes of the two transforming Theileria species, most notably, the TashAT/TpHN and Tar/Tpr gene families. Gene families that are present only in T. parva and T. annulata and not in T. orientalis, Babesia bovis, or Plasmodium were also identified. Identification of differences between the genome sequences of Theileria species with different abilities to transform and immortalize bovine leukocytes will provide insight into proteins and mechanisms that have evolved to induce and regulate this process. The T. orientalis genome database is available at http://totdb.czc.hokudai.ac.jp/. IMPORTANCE Cancer-like growth of leukocytes infected with malignant Theileria parasites is a unique cellular event, as it involves the transformation and immortalization of one eukaryotic cell by another. In this study, we sequenced the whole genome of a nontransforming Theileria species, Theileria orientalis, and compared it to the published sequences representative of two malignant, transforming species, T. parva and T. annulata. The genome-wide comparison of these parasite species highlights significant genetic diversity that may be associated with evolution of the mechanism(s) deployed by an intracellular eukaryotic parasite to transform its host cell. Cancer-like growth of leukocytes infected with malignant Theileria parasites is a unique cellular event, as it involves the transformation and immortalization of one eukaryotic cell by another. In this study, we sequenced the whole genome of a nontransforming Theileria species, Theileria orientalis, and compared it to the published sequences representative of two malignant, transforming species, T. parva and T. annulata. The genome-wide comparison of these parasite species highlights significant genetic diversity that may be associated with evolution of the mechanism(s) deployed by an intracellular eukaryotic parasite to transform its host cell.

The Autophagy Database: an all-inclusive information resource on autophagy that provides nourishment for research

Autophagy is a process of self-digestion generally observed in eukaryotes and has been shown to play crucial roles for survival under starvation and removal of deleterious substances. Despite great advances that have been made, many problems in mechanisms of autophagy remain unsolved. As a large number of autophagy-related proteins are identified in each species, a database that collects data, identifies their homologs in other species and makes them available will contribute to research advancement. As no such resources exist, we built the Autophagy database (http://tp-apg.genes.nig.ac.jp/autophagy) to provide basics, up-to-date information on relevant literature, and a list of autophagy-related proteins and their homologs in 41 eukaryotes. From the database, the user can search for proteins by keywords or sequences to obtain a wealth of data including functional and structural information and find possible functional homologs of proteins whose functions have been demonstrated in other species. As proteins that bind the phospholipid, phosphatidyl inositol 3-phosphate (PI3P) are essential for autophagy to proceed, we carried out an original analysis to identify probable PI3P-binding proteins, and made the list available from the database. The database is expected to give impetus to further research on autophagy by providing basic and specialized data on the subject.

Computational Prediction of O-linked Glycosylation Sites that Preferentially Map on Intrinsically Disordered Regions of Extracellular Proteins

O-glycosylation of mammalian proteins is one of the important posttranslational modifications. We applied a support vector machine (SVM) to predict whether Ser or Thr is glycosylated, in order to elucidate the O-glycosylation mechanism. O-glycosylated sites were often found clustered along the sequence, whereas other sites were located sporadically. Therefore, we developed two types of SVMs for predicting clustered and isolated sites separately. We found that the amino acid composition was effective for predicting the clustered type, whereas the site-specific algorithm was effective for the isolated type. The highest prediction accuracy for the clustered type was 74%, while that for the isolated type was 79%. The existence frequency of amino acids around the O-glycosylation sites was different in the two types: namely, Pro, Val and Ala had high existence probabilities at each specific position relative to a glycosylation site, especially for the isolated type. Independent component analyses for the amino acid sequences around O-glycosylation sites showed the position-specific existences of the identified amino acids as independent components. The O-glycosylation sites were preferentially located within intrinsically disordered regions of extracellular proteins: particularly, more than 90% of the clustered O-GalNAc glycosylation sites were observed in intrinsically disordered regions. This feature could be the key for understanding the non-conservation property of O-glycosylation, and its role in functional diversity and structural stability.

Development of an accurate classification system of proteins into structured and unstructured regions that uncovers novel structural domains: its application to human transcription factors

BackgroundIn addition to structural domains, most eukaryotic proteins possess intrinsically disordered (ID) regions. Although ID regions often play important functional roles, their accurate identification is difficult. As human transcription factors (TFs) constitute a typical group of proteins with long ID regions, we regarded them as a model of all proteins and attempted to accurately classify TFs into structural domains and ID regions. Although an extremely high fraction of ID regions besides DNA binding and/or other domains was detected in human TFs in our previous investigation, 20% of the residues were left unassigned. In this report, we exploit the generally higher sequence divergence in ID regions than in structural regions to completely divide proteins into structural domains and ID regions.ResultsThe new dichotomic system first identifies domains of known structures, followed by assignment of structural domains and ID regions with a combination of pre-existing tools and a newly developed program based on sequence divergence, taking un-aligned regions into consideration. The system was found to be highly accurate: its application to a set of proteins with experimentally verified ID regions had an error rate as low as 2%. Application of this system to human TFs (401 proteins) showed that 38% of the residues were in structural domains, while 62% were in ID regions. The preponderance of ID regions makes a sharp contrast to TFs of Escherichia coli (229 proteins), in which only 5% fell in ID regions. The method also revealed that 4.0% and 11.8% of the total length in human and E. coli TFs, respectively, are comprised of structural domains whose structures have not been determined.ConclusionThe present system verifies that sequence divergence including information of unaligned regions is a good indicator of ID regions. The system for the first time estimates the complete fractioning of structured/un-structured regions in human TFs, also revealing structural domains without homology to known structures. These predicted novel structural domains are good targets of structural genomics. When applied to other proteins, the system is expected to uncover more novel structural domains.

The GTOP database in 2009: updated content and novel features to expand and deepen insights into protein structures and functions

The Genomes TO Protein Structures and Functions (GTOP) database (http://spock.genes.nig.ac.jp/~genome/gtop.html) freely provides an extensive collection of information on protein structures and functions obtained by application of various computational tools to the amino acid sequences of entirely sequenced genomes. GTOP contains annotations of 3D structures, protein families, functions, and other useful data of a protein of interest in user-friendly ways to give a deep insight into the protein structure. From the initial 1999 version, GTOP has been continually updated to reap the fruits of genome projects and augmented to supply novel information, in particular intrinsically disordered regions. As intrinsically disordered regions constitute a considerable fraction of proteins and often play crucial roles especially in eukaryotes, their assignments give important additional clues to the functionality of proteins. Additionally, we have incorporated the following features into GTOP: a platform independent structural viewer, results of HMM searches against SCOP and Pfam, secondary structure predictions, color display of exon boundaries in eukaryotic proteins, assignments of gene ontology terms, search tools, and master files.

Parvulin (Par14), a Peptidyl-Prolyl cis-trans Isomerase, Is a Novel rRNA Processing Factor That Evolved in the Metazoan Lineage

Although parvulin (Par14/eukaryotic parvulin homolog), a peptidyl-prolyl cis-trans isomerase, is found associated with the preribosomal ribonucleoprotein (pre-rRNP) complexes, its roles in ribosome biogenesis remain undetermined. In this study, we describe a comprehensive proteomics analysis of the Par14-associated pre-rRNP complexes using LC-MS/MS and a knockdown analysis of Par14. Together with our previous results, we finally identified 115 protein components of the complexes, including 39 ribosomal proteins and 54 potential trans-acting factors whose yeast homologs are found in the pre-rRNP complexes formed at various stages of ribosome biogenesis. We give evidence that, although Par14 exists in both the phosphorylated and unphosphorylated forms in the cell, only the latter form is associated with the pre-40 S and pre-60 S ribosomal complexes. We also show that Par14 co-localizes with the nucleolar protein B23 during the interphase and in the spindle apparatus during mitosis and that actinomycin D treatment results in the exclusion of Par14 from the nucleolus. Finally we demonstrate that knockdown of Par14 mRNA decelerates the processing of pre-rRNA to 18 and 28 S rRNAs. We propose that Par14 is a component of the pre-rRNA complexes and functions as an rRNA processing factor in ribosome biogenesis. As the amino acid sequence of Par14 including that in the amino-terminal pre-rRNP binding region is conserved only in metazoan homologs, we suggest that its roles in ribosome biogenesis have evolved in the metazoan lineage.

A systematic investigation identifies a significant number of probable pseudogenes in the Escherichia coli genome.

Pseudogenes are open reading frames (ORFs) encoding dysfunctional proteins with high homology to known protein-coding genes. Although pseudogenes were reported to exist in the genomes of many eukaryotes and bacteria, no systematic search for pseudogenes in the Escherichia coli genome has been carried out. Genome comparisons of E. coli strains K-12 and O157 revealed that many protein-coding sequences have prematurely terminated orthologs encoding unstable proteins. To systematically screen for pseudogenes, we selected ORFs generated by premature termination of the orthologous protein-coding genes and subsequently excluded those possibly arising from sequence errors. Lastly we eliminated those with close homologs in this and other species, as these shortened ORFs may actually have functions. The process produced 95 and 101 pseudogene candidates in K-12 and O157, respectively. The assigned three-dimensional structures suggest that most of the encoded proteins cannot fold properly and thus are dysfunctional, indicating that they are probably pseudogenes. Therefore, the existence of a significant number of probable pseudogenes in E. coli is predicted, awaiting experimental verification. Most of them were found to be genes with paralogs or horizontally transferred genes or both. We suggest that pseudogenes constitute a small fraction of the genomes of free-living bacteria in general, reflecting the faster elimination than production of pseudogenes.

Splicing Factor 2-Associated Protein p32 Participates in Ribosome Biogenesis by Regulating the Binding of Nop52 and Fibrillarin to Preribosome Particles

Ribosome biogenesis starts with transcription of the large ribosomal RNA precursor (47S pre-rRNA), which soon combines with numerous factors to form the 90S pre-ribosome in the nucleolus. Although the subsequent separation of the pre-90S particle into pre-40S and pre-60S particles is critical for the production process of mature small and large ribosomal subunits, its molecular mechanisms remain undetermined. Here, we present evidence that p32, fibrillarin (FBL), and Nop52 play key roles in this separation step. Mass-based analyses combined with immunoblotting showed that p32 associated with 155 proteins including 31 rRNA-processing factors (of which nine were components of small subunit processome, and six were those of RIX1 complex), 13 chromatin remodeling components, and six general transcription factors required for RNA polymerase III-mediated transcription. Of these, a late rRNA-processing factor Nop52 interacted directly with p32. Immunocytochemical analyses demonstrated that p32 colocalized with an early rRNA-processing factor FBL or Nop52 in the nucleolus and Cajal bodies, but was excluded from the nucleolus after actinomycin D treatment. p32 was present in the pre-ribosomal fractions prepared by cell fractionation or separated by ultracentrifugation of the nuclear extract. p32 also associated with pre-rRNAs including 47S/45S and 32S pre-rRNAs. Furthermore, knockdown of p32 with a small interfering RNA slowed the early processing from 47S/45S pre-rRNAs to 18S rRNA and 32S pre-rRNA. Finally, Nop52 was found to compete with FBL for binding to p32 probably in the nucleolus. Given the fact that FBL and Nop52 are associated with pre-ribosome particles distinctly different from each other, we suggest that p32 is a new rRNA maturation factor involved in the remodeling from pre-90S particles to pre-40S and pre-60S particles that requires the exchange of FBL for Nop52.

Mechanistic insights into ectodomain shedding: Susceptibility of CADM1 adhesion molecule is determined by alternative splicing and O-glycosylation

Ectodomain shedding (shedding) is a post-translational modification, which liberates the extracellular domain of membrane proteins through juxtamembrane processing executed mainly by the ADAM (a disintegrin and metalloprotease) family of metalloproteases. Because shedding alters characteristics of cells in a rapid and irreversible manner, it should be strictly regulated. However, the molecular mechanisms determining membrane protein susceptibility to shedding (shedding susceptibility) are largely unknown. Here we report that alternative splicing can give rise to both shedding-susceptible and shedding-resistant CADM1 (cell adhesion molecule 1) variant proteins. We further show that O-glycans adjacent to the shedding cleavage site interfere with CADM1 shedding, and the only 33-bp alternative exon confers shedding susceptibility to CADM1 by inserting five non-glycosylatable amino acids between interfering O-glycans and the shedding cleavage site. These results demonstrate that shedding susceptibility of membrane protein can be determined at two different levels of its biosynthesis pathway, alternative splicing and O-glycosylation.