Details of a Researcher - Tomita, Masaru

Papers - Tomita, Masaru

Division display 381 - 400 of about 520 ／ All the affair displays >>

Validating the significance of genomic properties of Chi sites from the distribution of all octamers in Escherichia coli

Arakawa K., Uno R., Nakayama Y., Tomita M.

Gene 392 ( 1-2 ) 239 - 246 2007.05

ISSN 03781119

Chi sites (5′-GCTGGTGG-3′) are homologous recombinational hotspot octamer sequences, which attenuate the exonuclease activity of RecBCD in Escherichia coli. They are overrepresented in the genome (1008 occurrences), preferentially located within coding regions (98%), oriented in the direction of replication (75%), and occur most commonly on the mRNA-synonymous sense strand of the double helix (79%). Previous statistical studies of the genome sequence suggested that these genomic properties of Chi sites appear to be related to their role in recombinational repair and therefore to replication and transcription. In this study, we employ three mathematical models to predict the properties of Chi sites from single nucleotide and multi-nucleotide compositions, and validate them statistically using the distribution of all octamer sequences in the entire genome, or exclusively within ORFs. The model based on the overall distribution of all octamers provided better predictions than the single nucleotide composition model, and the ORF and sense strand preference of Chi sites were shown to be within the standard deviation of all octamers. In contrast, the orientation bias of the Chi sites in the direction of replication was significant, although the bias was not as pronounced as with the single nucleotide composition model, suggesting a selective pressure related to the role of RecBCD in replication. © 2007 Elsevier B.V. All rights reserved.

Access to Document (DOI)

Roles of hemoglobin allostery in hypoxia-induced metabolic alterations in erythrocytes: Simulation and its verification by metabolome analysis

Kinoshita A., Tsukada K., Soga T., Hishiki T., Ueno Y., Nakayama Y., Tomita M., Suematsu M.

Journal of Biological Chemistry 282 ( 14 ) 10731 - 10741 2007.04

ISSN 00219258

When erythrocytes are exposed to hypoxia, hemoglobin (Hb) stabilizes in theT-state by capturing 2,3-bisphosphoglycerate. This process could reduce the intracellular pool of glycolytic substrates, jeopardizing cellular energetics. Recent observations suggest that hypoxia-induced activation of glycolytic enzymes is correlated with their release from Band III (BIII) on the cell membrane. Based on these data, we developed a mathematical model of erythrocyte metabolism and compared hypoxia-induced differences in predicted activities of the enzymes, their products, and cellular energetics between models with and without the interaction of Hb with BIII. The models predicted that the allostery-dependent Hb interaction with BIII accelerates consumption of upstream glycolytic substrates such as glucose 6-phosphate and increases downstream products such as phosphoenolpyruvate. This prediction was consistent with metabolomic data from capillary electrophoresis mass spectrometry. The hypoxia-induced alterations in the metabolites resulted from acceleration of glycolysis, as judged by increased conversion of [13C]glucose to [13C]lactate. The allostery-dependent interaction of Hb with BIII appeared to contribute not only to maintenance of energy charge but also to further synthesis of 2,3-bisphosphoglycerate, which could help sustain stabilization of T-state Hb during hypoxia. Furthermore, such an activation of glycolysis was not observed when Hb was stabilized in R-state by treating the cells with CO. These results suggest that Hb allostery in erythrocytes serves as an O2-sensing trigger that drives glycolytic acceleration to stabilize intracellular energetics and promote the ability to release O 2 from the cells. © 2007 by The American Society for Biochemistry and Molecular Biology, Inc.

Access to Document (DOI)

Simulation study of methemoglobin reduction in erythrocytes: Differential contributions of two pathways to tolerance to oxidative stress

Kinoshita A., Nakayama Y., Kitayama T., Tomita M.

FEBS Journal 274 ( 6 ) 1449 - 1458 2007.03

ISSN 1742464X

Methemoglobin (metHb), an oxidized form of hemoglobin, is unable to bind and carry oxygen. Erythrocytes are continuously subjected to oxidative stress and nitrite exposure, which results in the spontaneous formation of metHb. To avoid the accumulation of metHb, reductive pathways mediated by cytochrome b5 or flavin, coupled with NADH-dependent or NADPH-dependent metHb reductases, respectively, keep the level of metHb in erythrocytes at less than 1% of the total hemoglobin under normal conditions. In this work, a mathematical model has been developed to quantitatively assess the relative contributions of the two major metHb-reducing pathways, taking into consideration the supply of NADH and NADPH from central energy metabolism. The results of the simulation experiments suggest that these pathways have different roles in the reduction of metHb; one has a high response rate to hemoglobin oxidation with a limited reducing flux, and the other has a low response rate with a high capacity flux. On the basis of the results of our model, under normal oxidative conditions, the NADPH-dependent system, the physiological role of which to date has been unclear, is predicted to be responsible for most of the reduction of metHb. In contrast, the cytochrome b5-NADH pathway becomes dominant under conditions of excess metHb accumulation, only after the capacity of the flavin-NADPH pathway has reached its limit. We discuss the potential implications of a system designed with two metHb-reducing pathways in human erythrocytes. © 2007 The Authors.

Access to Document (DOI)

A mathematical model for the Kai-protein-based chemical oscillator and clock gene expression rhythms in cyanobacteria

Miyoshi, Fumihiko ; Nakayama, Yoichi ; Kaizu, Kazunari ; Iwasaki, Hideo ; Tomita, Masaru

JOURNAL OF BIOLOGICAL RHYTHMS 22 ( 1 ) 69-80 - 80 2007.02

Research paper (scientific journal), Joint Work, ISSN 0748-7304

In the cyanobacterium, Synechococcus elongatus, most promoters are regulated by a circadian clock under continuous light (LL) conditions. Nevertheless, the basic circadian oscillation is primarily generated by alternating KaiC phosphorylation/dephosphorylation reactions at the posttranslational level. Indeed, the KaiC phosphorylation cycle was recently reconstituted in vitro by incubating KaiA, KaiB, and KaiC proteins with ATP. However, the molecular dynamics of this chemical oscillation and the mechanism that drives the circadian transcription/translation rhythms remain unknown. In this report, the KaiC phosphorylation cycle and the gene regulatory network in the cyanobacterial circadian system have been modeled. The model reproduces the robust KaiC phosphorylation cycle in the absence of de novo gene expression as is observed in vitro, as well as its coupling to transcriptional/translational feedback in LL conditions in vivo. Moreover, the model is consistent with most previous experiments, including various combinations of genetic knockout or overexpression of kai genes. It also predicts that multiple KaiC phosphorylation states and dynamic Kai protein interactions may be required for the cyanobacterial circadian system. © 2007 Sage Publications.

Access to Document (DOI)

Noise-reduction filtering for accurate detection of replication termini in bacterial genomes

Arakawa, Kazuharu ; Saito, Rintaro ; Tomita, Masaru

FEBS LETTERS 581 ( 2 ) 253-258 - 258 2007.01

Research paper (scientific journal), Joint Work, ISSN 0014-5793

Bacterial chromosomes are highly polarized in their nucleotide composition through mutational selection related to replication. Using compositional skews such as the GC skew, replication origin and terminus can be predicted in silico by observing the shift points. However, the genome sequence is affected by myriad functional requirements and selection on numerous subgenomic features, and elimination of this "noise" should lead to better predictions. Here, we present a noise-reduction approach that uses low-pass filtering through Fast Fourier transform coupled with cumulative skew graphs. It increases the prediction accuracy of the replication termini compared with previously documented methods based on genomic base composition. © 2006 Federation of European Biochemical Societies.

Access to Document (DOI)

eXpanda: an Integrated Platform for Network Analysis and Visualization

Negishi, Y., Nakamura, H., Yachie, N., Saito, R. and Tomita, M.

In silico Biology 7 ( 2 ) 141 - 144 2007

ISSN 13866338

Analysis and visualization of biological networks, such as protein-protein and protein-DNA interactions, are crucially important toward obtaining a thorough understanding of living systems. Here, we present an integrative software platform, eXpanda, which enables an analysis of a very broad range of biological networks, with a special focus on the extraction of characteristic topologies which potentially function as units in the networks. eXpanda is provided as a Perl library which gives full-automatic connections to various biological databases via a Perl programmable interface and can perform topological analysis based on graph theory. The results of these analyses are visualizable by vector graphics. eXpanda is under GNU General Public License. Software package, detailed documentations, source codes, and some sample scripts are downloadable at http://medcd.iab.keio.ac.jp/expanda/. © 2007 IOS Press. All rights reserved.

Direct measurement of isotopomer of intracellular metabolites using capillary electrophoresis time-of-flight mass spectrometry for efficient metabolic flux analysis

Toya, Y., Ishii, N., Hirasawa, T., Tomita, M. and Soga, T.

Journal of Chromatography A 1159 ( 1-2 ) 134–141 - 141 2007

Research paper (scientific journal), Joint Work, ISSN 00219673

We have developed a metabolic flux analysis method that is based on 13C-labeling patterns of the intracellular metabolites directly measured by capillary electrophoresis time-of-flight mass spectrometry (CE-TOFMS). The flux distribution of the central carbon metabolism in Escherichia coli was determined by this new approach and the results were compared with findings obtained by conventional GC-MS analysis based on isotopomer of the proteinogenic amino acids. There were some differences in estimation results between new approach using CE-TOFMS and conventional approach using GC-MS. These were thought to be attributable to variations in measured mass distributions between amino acids and the corresponding precursors and to differences in the sensitivity of the exchange coefficients to mass distributions. However, our CE-TOFMS method facilitates high-throughput flux analysis without requiring complicated sample preparation such as hydrolysis of proteins and derivatization of amino acids. © 2007 Elsevier B.V. All rights reserved.

Access to Document (DOI)

Analysis of nucleotides by pressure-assisted capillary electrophoresis　mass spectrometry using silanol mask technique

Soga, T., Ishikawa, T., Igarashi, S., Sugawara, K., Kakazu, Y. and Tomita, M.

Journal of Chromatography A 1159 ( 1-2 ) 125–133 - 133 2007

Research paper (scientific journal), Joint Work, ISSN 00219673

A method for the determination of nucleotides based on pressure-assisted capillary electrophoresis-electrospray ionization mass spectrometry (PACE-MS) is described. To prevent multi-phosphorylated species from adsorbing onto the fused-silica capillary, silanol groups were masked with phosphate ions by preconditioning the capillary with the background electrolyte containing phosphate. During preconditioning, nebulizer gas was turned off to avoid contamination of MS detector with phosphate ions. To detect nucleotides using the CE positive mode at a pH 7.5, it was necessary to apply air pressure to the inlet capillary during electrophoresis to supplement the electroosmotic flow (EOF) toward the cathode. Moreover, we exchanged the running electrolyte every analysis using the buffer replenishment system to obtain the required reproducibility. Under the optimized conditions, 14 phosphorylated species such as nucleotides, nicotinamide-adenine dinucleotides and coenzyme A (CoA) compounds were well determined in less than 20 min. The relative standard deviations (n = 6) of the method were better than 0.9% for migration times and between 1.7% and 8.1% for peak areas. The detection limits for these species were between 0.5 and 1.7 μmol/L with pressure injection of 50 mbar for 30 s (30 nL) at a signal-to-noise ratio of 3. This approach is robust and quantitative compared to the previous method, and its utility is demonstrated by the analysis of intracellular nucleotides and CoA compounds extracted from Escherichia coli wild type, pfkA and pfkB knockout mutants. The methodology was used to suggest that pfkA is the main functional enzyme. © 2007 Elsevier B.V. All rights reserved.

Access to Document (DOI)

Conjugational transfer system to shuttle giant DNA cloned by Bacillus subtilis genome (BGM) vector

Kuroki, A., Ohtani, N., Tsuge, K., Tomita, M. and Itaya, M.

Gene 399 ( 1 ) 72-80 - 80 2007

Research paper (scientific journal), Joint Work, ISSN 03781119

The Bacillus subtilis GenoMe (BGM) vector was designed as a versatile vector for the cloning of giant DNA segments. Cloned DNA in the BGM can be retrieved to a plasmid using our Bacillus recombinational transfer (BReT) method that takes advantage of competent cell transformation. However, delivery of the plasmid to a different B. subtilis strain by the normal transformation method is hampered by DNA size-related inefficiency. Therefore, we designed a novel method, conjugational plasmid-mediated DNA retrieval and transfer (CReT) from the BGM vector, and investigated conjugational transmission to traverse DNA between cells to circumvent the transformation-induced size limitation. pLS20, a 65-kb plasmid capable of conjugational transfer between B. subtilis strains, was modified to retrieve DNA cloned in the BGM vector by homologous recombination during normal culture. As the plasmid copy number was estimated to be 3, the retrieval plasmid was selected using increased numbers of marker genes derived from the retrieved DNA. We applied this method to retrieve Synechocystis genome segments up to 90 kb in length. We observed retrieved plasmid transfers between B. subtilis strains by conjugation in the absence of structural alterations in the DNA fragment. Our observations extend DNA transfer protocols over previously exploited size ranges. © 2007 Elsevier B.V. All rights reserved.

Access to Document (DOI)

Proteome-Wide Prediction of Novel DNA/RNA-Binding Proteins Using Amino Acid Composition and Periodicity in the Hyperthermophilic Archaeon Pyrococcus furiosus

Fujishima, K., Komasa, M., Kitamura, S., Suzuki, H., Tomita, M. and Kanai, A.

DNA Research 14 ( 3 ) 91-102 - 102 2007

Research paper (scientific journal), Joint Work, ISSN 13402838

Proteins play a critical role in complex biological systems, yet about half of the proteins in publicly available databases are annotated as functionally unknown. Proteome-wide functional classification using bioinformatics approaches thus is becoming an important method for revealing unknown protein functions. Using the hyperthermophilic archaeon Pyrococcus furiosus as a model species, we used the support vector machine (SVM) method to discriminate DNA/RNA-binding proteins from proteins with other functions, using amino acid composition and periodicities as feature vectors. We defined this value as the composition score (CO) and periodicity score (PD). The P. furiosus proteins were classified into three classes (I-III) on the basis of the two-dimensional correlation analysis of CO score and PD score. As a result, approximately 87 of the functionally known proteins categorized as class I proteins (CO score + PD score > 0.6) were found to be DNA/RNA-binding proteins. Applying the two-dimensional correlation analysis to the 994 hypothetical proteins in P. furiosus, a total of 151 proteins were predicted to be novel DNA/RNA-binding protein candidates. DNA/RNA-binding activities of randomly chosen hypothetical proteins were experimentally verified. Six out of seven candidate proteins in class I possessed DNA/RNA-binding activities, supporting the efficacy of our method. © The Author 2007. Kazusa DNA Research Institute.

Access to Document (DOI)

Evolution of ribonuclease H genes in prokaryotes to avoid inheritance of redundant genes

Kochiwa, H., Tomita, M. and Kanai, A.

BMC Evolutionary Biology 7 128 2007

Research paper (scientific journal), Joint Work

Background. A theoretical model of genetic redundancy has proposed that the fates of redundant genes depend on the degree of functional redundancy, and that functionally redundant genes will not be inherited together. However, no example of actual gene evolution has been reported that can be used to test this model. Here, we analyzed the molecular evolution of the ribonuclease H (RNase H) family in prokaryotes and used the results to examine the implications of functional redundancy for gene evolution. Results. In prokaryotes, RNase H has been classified into RNase HI, HII, and HIII on the basis of amino acid sequences. Using 353 prokaryotic genomes, we identified the genes encoding the RNase H group and examined combinations of these genes in individual genomes. We found that the RNase H group may have evolved in such a way that the RNase HI and HIII genes will not coexist within a single genome - in other words, these genes are inherited in a mutually exclusive manner. Avoiding the simultaneous inheritance of the RNase HI and HIII genes is remarkable when RNase HI contains an additional non-RNase H domain, double-stranded RNA, and an RNA-DNA hybrid-binding domain, which is often observed in eukaryotic RNase H1. This evolutionary process may have resulted from functional redundancy of these genes, because the substrate preferences of RNase HI and RNase HIII are similar. Conclusion. We provide two possible evolutionary models for RNase H genes in which functional redundancy contributes to the exclusion of redundant genes from the genome of a species. This is the first empirical study to show the effect of functional redundancy on changes in gene constitution during the course of evolution. © 2007 Kochiwa et al; licensee BioMed Central Ltd.

Access to Document (DOI)

Indeterminacy of Reverse Engineering of Gene Regulatory Networks: The Curse of Gene Elasticity

Krishnan, A., Giuliani, A. and Tomita, M

PLoS ONE 2 ( 6 ) e562 2007

Research paper (scientific journal), Joint Work

Background. Gene Regulatory Networks (GRNs) have become a major focus in interest in recent years. A number of reverse engineering approaches have been developed to help uncover the regulatory networks giving rise to the observed gene expression profiles. However, this is an overspecified problem due to the fact that more than one genotype (network wiring) can give rise to the same phenotype. We refer to this phenomenon as "gene elasticity." In this work, we study the effect of this particular problem on the pure, data-driven inference of gene regulatory networks. Methodology. We simulated a four-gene network in order to produce "data" (protein levels) that we use in lieu of real experimental data. We then optimized the network connections between the four genes with a view to obtain the original network that gave rise to the data. We did this for two different cases: one in which only the network connections were optimized and the other in which both the network connections as well as the kinetic parameters (given as reaction probabilities in our case) were estimated. We observed that multiple genotypes gave rise to very similar protein levels. Statistical experimentation indicates that it is impossible to differentiate between the different networks on the basis of both equilibrium as well as dynamic data. Conclusions. We show explicitly that reverse engineering of GRNs from pure expression data is an indeterminate problem. Our results suggest the unsuitability of an inferential, purely data-driven approach for the reverse engineering transcriptional networks in the case of gene regulatory networks displaying a certain level of complexity. © 2007 Krishnan et al.

Access to Document (DOI)

Restauro-G: A Rapid Genome Re-Annotation System for Comparative Genomics

Tamaki, S., Arakawa, K., Kono, N. and Tomita, M.

Genomics Proteomics Bioinformatics 5; 53-58(2007) 5 ( 1 ) 53-58 - 58 2007

Research paper (scientific journal), Joint Work, ISSN 16720229

Annotations of complete genome sequences submitted directly from sequencing projects are diverse in terms of annotation strategies and update frequencies. These inconsistencies make comparative studies difficult. To allow rapid data preparation of a large number of complete genomes, automation and speed are important for genome re-annotation. Here we introduce an open-source rapid genome re-annotation software system, Restauro-G, specialized for bacterial genomes. Restauro-G re-annotates a genome by similarity searches utilizing the BLAST-Like Alignment Tool, referring to protein databases such as UniProt KB, NCBI nr, NCBI COGs, Pfam, and PSORTb. Re-annotation by Restauro-G achieved over 98% accuracy for most bacterial chromosomes in comparison with the original manually curated annotation of EMBL releases. Restauro-G was developed in the generic bioinformatics workbench G-language Genome Analysis Environment and is distributed at http://restauro-g.iab.keio.ac.jp/under the GNU General Public License. © 2007 Beijing Genomics Institute.

Access to Document (DOI)

Computational analysis of microRNA-mediated antiviral defense in humans

Watanabe, Y., Kishi, A., Yachie, N., Kanai, A. and Tomita, M.

FEBS Letters 581 ( 24 ) 4603-4610 - 4610 2007

Research paper (scientific journal), Joint Work, ISSN 00145793

Recent studies have proposed the interesting perspective that viral gene expression is downregulated by host microRNAs (miRNAs), small non-coding RNAs well known as post-transcriptional gene regulators. We computationally predicted human miRNA target sites within 228 human-infecting and 348 invertebrate-infecting virus genomes, and we observed that human-infecting viruses were more likely than invertebrate-infecting ones to be targeted by human miRNAs. We listed 62 possible human miRNA-targeted viruses from 6 families, most of which consisted of single-stranded RNA viruses. These results suggest that miRNAs extensively mediate antiviral defenses in humans. © 2007 Federation of European Biochemical Societies.

Access to Document (DOI)

Polymer Entrapment in Polymerized Silicate for Preparing Highly Stable Capillary Coatings for CE and CE-MS

Monton, M., Tomita, M., Soga, T. and Ishihama, Y.

Analytical Chemistry 79 ( 20 ) 7838-7844 - 7844 2007

Research paper (scientific journal), Joint Work, ISSN 00032700

An easy-to-implement capillary coating strategy based on polymer entrapment in the network of polymerized silicate is described. In this manner, cationic polymers are tightiy fixed onto the inner wall of the capillary for electroosmotic flow control without necessitating complex surface modification chemistries. The resulting coated capillary exhibited good stability over a wide range of pH, good reproducibility, strong endurance in more than 300 electrophoretic runs, and tolerance of commonly employed organic solvent additives in CE. Applications in CE-MS analysis of biologically important anions as well as sample enrichment are shown. Additionally, it was used as a durable base for attachment of multiple layers of charged polymers on the wall, via electrostatic interaction with the preceding layer. Thus, two novel types of highly stable coated capillaries, one with anodic EOF and the other cathodic, were developed. © 2007 American Chemical Society.

Access to Document (DOI)

Phosphopeptide Enrichment by Aliphatic Hydroxy Acid-modified Metal Oxide Chromatography for Nano-LC-MS/MS in Proteomics Applications

Sugiyama, N., Masuda, T., Shinoda, K., Nakamura, A., Tomita, M. and Ishihama, Y.

Mol Cell Proteomics 6 ( 6 ) 1103-1109 - 1109 2007

Research paper (scientific journal), Joint Work, ISSN 15359476

We developed novel methods for phosphopeptide enrichment using aliphatic hydroxy acid-modified metal oxide chromatography (MOC). Titania and zirconia were successfully applied to enrich phosphopeptides with the aid of aliphatic hydroxy acids, such as lactic acid and β-hydroxypropanoic acid, to reduce the interaction between acidic non-phosphopeptides and the metal oxides. These methods removed the vast majority of non-phosphopeptides from phosphoprotein standard digests, and large numbers of phosphopeptides could be readily identified. The methods were coupled with nano-LC-MS/MS systems without difficulty. Recovery of phosphopeptides in MOC varied greatly from peptide to peptide, ranging from a few percent to 100%, and the average was almost 50%. Repeatability and linearity were satisfactory. In an examination of the cytoplasmic fraction of HeLa cells, more than 1000 phosphopeptides were identified using lactic acid-modified titania MOC and β-hydroxypropanoic acid-modified zirconia MOC, respectively. The overlap between phosphopeptides enriched by these two methods was 40%, and the combined results provided 1646 unique phosphopeptides. To our knowledge, this is the first successful application of a single MOC-based approach to phosphopeptide enrichment from complex biological samples such as cell lysates. © 2007 by The American Society for Biochemistry and Molecular Biology, Inc.

Access to Document (DOI)

Statistical Analysis of Gene Expression in Innate Immune Responses: Dynamic Interactions between MicroRNA and Signaling Molecules

Piras, V., Selvarajoo, K., Fujikawa, N., Choi, S., Tomita, M., Giuliani, A. and Tsuchiya, M.

Genomics and Informatics 5 107-112 2007

Research paper (scientific journal), Joint Work

Multiple high throughput analyses monitor the response of E. coli to perturbations

Ishii, N., Nakahigashi, K., Baba, T., Robert, M., Soga, T., Kanai, A., Hirasawa, T., Naba, M., Hirai, K., Hoque, A., Yee, Ho, Pei., Kakazu, Y., Sugawara, K., Igarashi, S., Harada, S., Masuda, T., Sugiyama, N., Togashi, T., Hasegawa, M., Takai, T., Yugi, K., Arakawa, K., Iwata, N., Toya, Y., Nakayama, Y., Nishioka, T., Shimizu, K., Mori, H. and Tomita, M.

Science 316 ( 5824 ) 593-597 - 597 2007

Research paper (scientific journal), Joint Work, ISSN 00368075

Analysis of cellular components at multiple levels of biological information can provide valuable functional insights. We performed multiple high-throughput measurements to study the response of Escherichia coli cells to genetic and environmental perturbations. Analysis of metabolic enzyme gene disruptants revealed unexpectedly small changes in messenger RNA and proteins for most disruptants. Overall, metabolite levels were also stable, reflecting the rerouting of fluxes in the metabolic network. In contrast, E. coli actively regulated enzyme levels to maintain a stable metabolic state in response to changes in growth rate. E. coli thus seems to use complementary strategies that result in a metabolic network robust against perturbations.

Access to Document (DOI)

Network scaling invariants help to elucidate basic topological principles of proteins

Krishnan, A., Giuliani, A., Zbilut, J.P. and Tomita, M.

Journal of Proteome Research 6 ( 10 ) 3924-3934 - 3934 2007

Research paper (scientific journal), Joint Work, ISSN 15353893

The structural architecture of proteins continues to be an area of active research. Despite the difference in models dealing with the way proteins fold into their tertiary structures, it is recognized that small regions of proteins tend to fold independently and are then stabilized by interactions between these distinct subunits. However, there are a number of different definitions of what comprises an independent subunit. In the belief that an unequivocal definition of a domain must be based on the most fundamental property of protein 3D structure, namely, the adjacency matrix of inter-residues contact, we adopt a network representation of the protein. In this work, we used a well-established, global method for identifying modules in networks, without any specific reference to the kind of network being analyzed. The algorithm converges toward the maximization of the modularity of the given protein network and, in doing so, allows the representation of the residues of the protein in terms of their intramodule degree, z, and participation coefficient, P. We demonstrate that the labeling of residues in terms of these invariants allows for information-rich representations of the studied proteins as well as to sketch a new way to link sequence, structure, and the dynamical properties of proteins. We discovered a strong invariant character of protein molecules in terms of P/z characterization, pointing to a common topological design of all protein structures. This invariant representation, applied to different protein systems, enabled us to identify the possible functional role of high P/z residues during the folding process. Additionally, we observe a hierarchical behavior of protein structural organization that provides a sequence-secondary-tertiary structure link. The discovery of similar and repeatable scaling laws at different level of definitions going from hydrophobicity patterning along the sequence up to the size of an autonomous folding unit (AFU) and general contact distribution of the entire molecule suggest a hierarchical-like behavior of protein architecture. This implies the possibility to select different privileged scales of observation for deriving useful information on protein systems. © 2007 American Chemical Society.

Access to Document (DOI)

Inferring rules of Escherichia coli translational efficiency using an artificial neural network

Mori, K., Saito, R., Kikuchi, S. and Tomita, M.

BioSystems 90 ( 2 ) 414-420 - 420 2007

Research paper (scientific journal), Joint Work, ISSN 03032647

Although the machinery for translation initiation in Escherichia coli is very complicated, the translational efficiency has been reported to be predictable from upstream oligonucleotide sequences. Conventional models have difficulties in their generalization ability and prediction nonlinearity and in their ability to deal with a variety of input attributions. To address these issues, we employed structural learning by artificial neural networks to infer general rules for translational efficiency. The correlation between translational activities measured by biological experiments and those predicted by our method in the test data was significant (r = 0.78), and our method uncovered underlying rules of translational activities and sequence patterns from the obtained skeleton structure. The significant rules for predicting translational efficiency were (1) G- and A-rich oligonucleotide sequences, resembling the Shine-Dalgarno sequence, at positions -10 to -7; (2) first base A in the initiation codon; (3) transport/binding or amino acid metabolism gene function; (4) high binding energy between mRNA and 16S rRNA at positions -15 to -5. An additional inferred novel rule was that C at position -1 increases translational efficiency. When our model was applied to the entire genomic sequence of E. coli, translational activities of genes for metabolism and translational were significantly high. © 2006 Elsevier Ireland Ltd. All rights reserved.

Access to Document (DOI)

　 Previous page - Next page 　