Published papers

A unified mathematical description that expresses the characteristics of whole systems is necessary for an understanding of signal transduction cascades. In this study we explore an algebraic method, named extreme signaling flow, enhanced from the concept of extreme pathway, to analyze signal transduction systems. This method enables us to represent the long-term potentiation (LTP) and the long-term depression (LTD) of hippocampal neuronal plasticity in an integrated simulation model. The model is validated by comparing the results of redundancy, reaction participation, and in silico knockout analysis with biological knowledge available from the literature. The following properties are assumed in these computational analyses: (1) LTP is fault-tolerant under network modification, (2) protein kinase C and MAPK have numerous routes to LTP induction, (3) calcium-calmodulin kinase II has a few routes to LTP induction, and (4) calcineurin has many routes to LTD induction. These results demonstrate that our approach produces an integrated framework for analyzing properties of large-scale systems with complicated signal transduction. © 2008 Massachusetts Institute of Technology.

The circadian clock of Drosophila is a model pathway for research in biological dock mechanisms, both with traditional experimental approaches and with emerging systems biology approaches utilizing mathematical modeling and in silico computer simulation. Dynamic diurnal oscillations are achieved by the complex interaction of components as a system, and mathematical reconstruction his proven to be an invaluable means of understanding such systematic behavior. In this study, we implemented eight published models of the Drosqpbila circadian clock in Systems Biology Markup Language (SBML) for comparative systems biology studies using E-Cell Simulation Environment version 3, to examine the system-level requirements for the clock mechanism to be robust, by calculating the period and amplitude sensitivity coefficients with simulation experiments. While all models were generally robust as determined by the network topology of the oscillatory feedback loop structure, existing models place relatively strong emphasis on transcription regulation, although this is a limitation on robustness. We suggest that more comprehensive modeling including protein phosphorylation, polymerization, and nuclear transport with regard to amplitude sensitivity will be necessary for understanding the light entrainment and temperature compensation of circadian clocks. © 2008 Massachusetts Institute of Technology.

An enormous body of information has been obtained by molecular and cellular biology in the last half century, However, even these powerful approaches are not adequate when it comes to higher-level biological structures, such as tissues, organs, and individual organisms, because of the complexities involved. Thus, accumulation of data at the higher levels supports and broadens the context for that obtained on the molecular and cellular levels. Under such auspices, an attempt to clucidate mesoscopic and macroscopic subjects based on plentiful nanoscopic and microscopic data is ot great potential value. On the other hand, fully realistic simulation is impracticable because of the extensive cost entailed and enormous amount of data required. Abstraction and modeling that balance the dual requirements of prediction accuracy and manageable calculation cost are of great importance for systems biology. We have constructed an ammonia metabolism model of the hepatic lobule, a histological component of the liver, based on a single-hepatocyte model that consists of the biochemical kinetics of enzymes and transporters. To bring the calculation cost within reason, the porto-central axis, which is an elemental structure of the lobule, is defined as the systems biological unit of the liver, and is accordingly modeled. A model including both histological structure and position-specific gene expression of major enzymes largely represents the physiological dynamics of the hepatic lobule in nature. In addition, heterogeneous gene expression is suggested to have evolved to optimize the energy efficiency of ammonia detoxificattion at the macroscopic level, implying that approaches like this may elucidate how properties at the molecular and cellular levels, such as regulated gene expression, modify higher-level phenomena of multicellular tissue, organs, and organisms. © 2008 Massachusetts Institute of Technology.

Protein phosphorylation regulates a wide range of cellular processes. Here, we report the proteome-wide mapping of in vivo phosphorylation sites in Arabidopsis by using complementary phosphopeptide enrichment techniques coupled with high-accuracy mass spectrometry. Using unfractionated whole cell lysates of Arabidopsis, we identified 2597 phosphopeptides with 2172 high-confidence, unique phosphorylation sites from 1346 proteins. The distribution of phosphoserine, phosphothreonine, and phosphotyrosine sites was 85.0, 10.7, and 4.3%. Although typical tyrosine-specific protein kinases are absent in Arabidopsis, the proportion of phosphotyrosines among the phospho-residues in Arabidopsis is similar to that in humans, where over 90 tyrosine-specific protein kinases have been identified. In addition, the tyrosine phosphoproteome shows features distinct from those of the serine and threonine phosphoproteomes. Taken together, we highlight the extent and contribution of tyrosine phosphorylation in plants. © 2008 EMBO and Nature Publishing Group.

Proteins are a major regulatory component in complex biological systems. Among them, DNA/RNA-binding proteins, the key components of the central dogma of molecular biology, and membrane proteins, which are necessary for both signal transduction and metabolite transport, are suggested to be the most important protein families that arose in the early stage of life. In this study, we computationally analyzed the whole proteome data of six model species to overview the protein diversity in the three domains of life (Bacteria, Archaea and Eukaryota), especially focusing on the above two protein families. To compare the protein distribution among the six model species, we calculated various protein profiles: hydropathy, molecular weight, amino acid composition and periodicity for each protein. We found a domain-specific distribution of the proteome based on 2D correlation analysis of hydropathy and molecular weight. Further, the merged protein distribution of Archaea and other domains revealed many membrane proteins localized in Bacteria-specific regions with a high ratio of hydropathy and many DNA/RNA-binding proteins localized in Eukaryota-specific regions with a low ratio of hydropathy. Since about half of the proteins encoded in the genome are still functionally unknown, we further conducted Support Vector Machine (SVM)-based functional prediction using amino acid composition (CO score) and periodicity (PD score) as feature vectors to predict the overall number of DNA/RNA-binding proteins and membrane proteins in the proteome. Our estimation indicated that two functional categories occupy approximately 60% to 80% of the proteome, and further, the proportion of the two categories varied among the three domains of life, suggesting that the proteome has gone through different selective pressure during evolution.

Sulfite plays an important role in beer flavor stability. Although breeding of bottom-fermenting Saccharomyces strains that produce high levels of SO 2 is desirable, it is complicated by the fact that undesirable H 2S is produced as an intermediate in the same pathway. Here, we report the development of a high-level SO2-producing bottom-fermenting yeast strain by integrated metabolome and transcriptome analysis. This analysis revealed that O-acetylhomoserine (OAH) is the rate-limiting factor for the production of SO2 and H2S. Appropriate genetic modifications were then introduced into a prototype strain to increase metabolic fluxes from aspartate to OAH and from sulfate to SO 2, resulting in high SO2 and low H2S production. Spontaneous mutants of an industrial strain that were resistant to both methionine and threonine analogs were then analyzed for similar metabolic fluxes. One promising mutant produced much higher levels of SO2 than the parent but produced parental levels of H2S. Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Background. It has long been known that small regions of proteins tend to fold independently and are then stabilized by interactions between these distinct subunits or modules. Such units, also known as autonomous folding units (AFUs) or"foldons" play a key role in protein folding. A knowledge of such early folding units has diverse applications in protein engineering as well as in developing an understanding of the protein folding process. Such AFUs can also be used as model systems in order to study the structural organization of proteins. Results. In an earlier work, we had utilized a global network partitioning algorithm to identify modules in proteins. We had shown that these modules correlate well with AFUs. In this work, we have developed a webserver, GANDivAWeb, to identify early folding units or "foldons" in networks using the algorithm described earlier. The website has three functionalities: (a) It is able to display information on the modularity of a database of 1420 proteins used in the original work, (b) It can take as input an uploaded PDB file, identify the modules using the GANDivA algorithm and email the results back to the user and (c) It can take as input an uploaded PDB file and a results file (obtained from functionality (b)) and display the results using the embedded viewer. The results include the module decomposition of the protein, plots of cartoon representations of the protein colored by module identity and connectivity as well as contour plots of the hydrophobicity and relative accessible surface area (RASA) distributions. Conclusion. We believe that the GANDivAWeb server, will be a useful tool for scientists interested in the phenomena of protein folding as well as in protein engineering. Our tool not only provides a knowledge of the AFUs through a natural graph partitioning approach but is also able to identify residues that are critical during folding. It is our intention to use this tool to study the topological determinants of protein folding by analyzing the topological changes in proteins over the unfolding/folding pathways. © 2008 Laborde et al; licensee BioMed Central Ltd.

A novel class of small, noncoding RNAs called microRNAs (miRNAs) was recently identified, and up to 20%-30% of human genes are said to be regulated by the members of this class of RNA. miRNAs regulate their target genes by imperfect binding to the 3′ untranslated region (UTR) of the gene transcript, and computational predictions of the binding between miRNAs and target gene transcripts have been performed in order to determine which genes are regulated by these RNAs. Phylogenetic analysis has also been used as a powerful tool to predict miRNA target genes. Much emphasis has therefore been placed on studying the phylogenetic conservation and evolution of this novel type of gene regulator, although there still is much that is not known. Here, we propose a hypothesis of how human miRNAs optimized their gene regulation by adjusting their transcript levels, and how they evolved through specific selection of their target genes. We analyzed the correlation between the conservation of miRNAs among species and three features: the number of transcripts, the formation of duplicates, and the number of target genes. The number of miRNA transcripts and the formation of duplicates increased as the conservation rate increased. In contrast, the number of target genes decreased as the conservation rate increased. Therefore, we propose that miRNAs gradually gain an ability to regulate specific target genes when such regulation has a positive effect on the organism. As its pool of target genes is refined, the ability of an miRNA to regulate the genes may be stabilized by an increase in the miRNA transcript number and the formation of duplicates.

We have developed a new protocol for digesting hydrophobic proteins using trypsin with the aid of phase-transfer surfactants (PTS), such as sodium deoxycholate (SDC). SDC increases the solubility of hydrophobic proteins, enhances the activity of trypsin, and improves the accessibility to trypsin of proteins denatured during the extraction process. After digestion, SDC was successfully removed from the acidified solution containing tryptic peptides by adding a water-immiscible organic solvent, into which SDC was predominantly transferred, while the digested peptides remained in the aqueous phase. Compared with a protocol using an acid-labile surfactant, this PTS protocol increased the number of identified proteins and the recovery of hydrophobic peptides in the analysis of 400 ng of a membrane-enriched fraction of Escherichia coli. Application of the PTS protocol to 9.0 μg of a membrane-enriched pellet from human cervical cancer HeLa cells resulted in identification of a total of 1450 proteins, of which 764 (53%) were membrane proteins, by two-dimensional strong cation exchange (SCX)-C18 LC-MSMS with 5 SCX fractions. The distribution of the number of transmembrane domains in proteins identified in this study was in agreement with that in the IPI human database, suggesting that the PTS protocol can provide unbiased digestion of the membrane proteome. © 2008 American Chemical Society.

An exacerbated sensitivity to apparently minor stimuli and a general resilience of the entire system stay together side-by-side in biological systems. This apparent paradox can be explained by the consideration of biological systems as very strongly interconnected network systems. Some nodes of these networks, thanks to their peculiar location in the network architecture, are responsible for the sensitivity aspects, while the large degree of interconnection is at the basis of the resilience properties of the system. One relevant feature of the high degree of connectivity of gene regulation networks is the emergence of collective ordered phenomena influencing the entire genome and not only a specific portion of transcripts. The great majority of existing gene regulation models give the impression of purely local 'hard-wired' mechanisms disregarding the emergence of global ordered behavior encompassing thousands of genes while the general, genome wide, aspects are less known. Here we address, on a data analysis perspective, the discrimination between local and global scale regulations, this goal was achieved by means of the examination of two biological systems: innate immune response in macrophages and oscillating growth dynamics in yeast. Our aim was to reconcile the 'hard-wired' local view of gene regulation with a global continuous and scalable one borrowed from statistical physics. This reconciliation is based on the network paradigm in which the local 'hard-wired' activities correspond to the activation of specific crucial nodes in the regulation network, while the scalable continuous responses can be equated to the collective oscillations of the network after a perturbation. © 2008 Elsevier B.V. All rights reserved.

Volvulina Playfair is a colonial volvocalean genus with usually 8 or 16 cells forming hollow ellipsoidal colonies. In this study, a strain of Volvulina was isolated from Lake Tsukui, Kanagawa Prefecture, Japan. This strain was morphologically identified as Volvulina compacta Nozaki, which was previously known only from strains isolated from Nepal. Chloroplast multigene analyses using rbcL, atpB, psaA, psaB and psbC genes showed that the Japanese strain of V. compacta is closely related to the Nepalese strain of V. compacta and to an American strain of Volvulina pringsheimii Starr, but only distantly to Volvulina steinii Playfair previously reported from Japan. Although the Japanese strain was more closely related to V. pringsheimii than to Nepalese V. compacta, it was clearly distinguished from V. pringsheimii by its square cell shape in the top view. These results indicate the existence of cryptic species within V. compacta.

In metabolomics, the rapid identification of quantitative differences between multiple biological samples remains a major challenge. While capillary electrophoresis-mass spectrometry (CE-MS) is a powerful tool to simultaneously quantify charged metabolites, reliable and easy-to-use software that is well suited to analyze CE-MS metabolic profiles is still lacking. Optimized software tools for CE-MS are needed because of the sometimes large variation in migration time between runs and the wider variety of peak shapes in CE-MS data compared with LC-MS or GC-MS. Therefore, we implemented a stand-alone application named JDAMP (Java application for Differential Analysis of Metabolite Profiles), which allows users to identify the metabolites that vary between two groups. The main features include fast calculation modules and a file converter using an original compact file format, baseline subtraction, dataset normalization and alignment, visualization on 2D plots (m/z and time axis) with matching metabolite standards, and the detection of significant differences between metabolite profiles. Moreover, it features an easy-to-use graphical user interface that requires only a few mouse-actions to complete the analysis. The interface also enables the analyst to evaluate the semiautomatic processes and interactively tune options and parameters depending on the input datasets. The confirmation of findings is available as a list of overlaid electropherograms, which is ranked using a novel difference-evaluation function that accounts for peak size and distortion as well as statistical criteria for accurate difference-detection. Overall, the JDAMP software complements other metabolomics data processing tools and permits easy and rapid detection of significant differences between multiple complex CE-MS profiles. © Springer Science+Business Media, LLC 2009.

We isolated strains of a new species of Chlorogonium (Cg.), Cg. complexum Nakada sp. nov. The genus Chlorogonium was recently emended to have several, irregularly distributed contractile vacuoles, two-to-many pyrenoids associated with fragmented starch plates and tubular thylakoid lamellae penetrating their matrices, and stigma composed of a single layer of globules. Chlorogonium complexum has these diagnostic characters of the emended genus. It is distinguished from other species of Chlorogonium, Cg. euchlorum, Cg. elongatum and Cg. capillatum, on the combination of diagnostic characters, such as the number and stability of pyrenoids and the mode of encystment. Based on a dataset combining 18S rRNA, RuBisCO large subunit (rbcL) and P700 chl a-apoprotein A2 (psaB) gene sequences, the phylogenetic position of Chlorogonium was determined. Chlorogonium complexum represented an independent lineage, sister to Cg. euchlorum and Cg. capillatum. The evolution of diagnostic characters among Chlorogonium species is discussed based on the molecular phylogenetic analyses, and a key to species of Chlorogonium and fusiform chlamydomonads with irregularly distributed contractile vacuoles is presented. © 2010 British Phycological Society.

Cephalomonas (Cp.) granulata N. L. Higinb. (Volvocales, Chlorophyceae), a rare volvocalean phytoflagellate, has recently been isolated from a paddy field in Japan and examined by light and fluorescence microscopy. The vegetative cells of Cp. granulata have granulate loricae and a characteristic mushroom-like shape with a hemispherical to spherical anterior part and narrower posterior part. Cp. granulata has been classified among the Phacotaceae along with other loricated genera. However, its phylogenetic position has not been investigated using molecular phylogeny. To evaluate the phylogenetic position of Cephalomonas, the 18S rRNA gene sequence of the Japanese strain was determined. The phylogenetic analyses revealed that Cp. granulata was not closely related to other loricated taxa with known phylogenetic affinities, such as Phacotaceae sensu stricto (Phacotus, Pteromonas, and Wislouchiella) and Dysmorphococcus globosus. This indicates at least three independent origins of loricae within the Volvocales. © 2010 Japanese Society of Phycology.

Since most bio-production processes are conducted in a batch or fed-batch manner, the evaluation of metabolism with respect to time is highly desirable. Toward this aim, we applied 13C-metabolic flux analysis to nonstationary conditions by measuring the mass isotopomer distribution of intracellular metabolites. We performed our analysis on batch cultures of wild-type Escherichia coli, as well as on Pyk and Pgi mutants, obtained the fluxes and metabolite concentrations as a function of time. Our results for the wild-type indicated that the TCA cycle flux tended to increase during growth on glucose. Following glucose exhaustion, cells controlled the branch ratio between the glyoxylate pathway and the TCA cycle, depending on the availability of acetate. In the Pyk mutant, the concentrations of glycolytic intermediates changed drastically over time due to the dumping and feedback inhibition caused by PEP accumulation. Nevertheless, the flux distribution and free amino acid concentrations changed little. The growth rate and the fluxes remained constant in the Pgi mutant and the glucose-6-phosphate dehydrogenase reaction was the rate-limiting step. The measured fluxes were compared with those predicted by flux balance analysis using maximization of biomass yield or ATP production. Our findings indicate that the objective function of biosynthesis became less important as time proceeds on glucose in the wild-type, while it remained highly important in the Pyk mutant. Furthermore, ATP production was the primary objective function in the Pgi mutant. This study demonstrates how cells adjust their metabolism in response to environmental changes and/or genetic perturbations in the batch cultivation. © 2010 American Institute of Chemical Engineers.

Chlamydomonas Ehrenb. is a unicellular volvocalean genus consisting of 400-600 species, most of which are known solely based on microscopy. In this study, a newly isolated strain of Chlamydomonas neoplanoconvexa (M. O. P. Iyengar) Nakada nom. nov. (≡Chlamydomonas planoconvexa M. O. P. Iyengar non J. W. G. Lund) was examined using light microscopy and combined 18S rRNA, rbcL and psaB gene phylogeny. The C. neoplanoconvexa strain was quite similar to Chlamydomonas perpusilla (Korshikov) Gerloff in terms of its small fusiform vegetative cells (∼10μm in length), parietal chloroplast containing a single pyrenoid, and nucleus posterior to the pyrenoid. However, C. neoplanoconvexa was distinguished from C. perpusilla based on the morphology of the papillae and pyrenoids. The two species belong to the clade Caudivolvoxa in the order Volvocales, but are distantly related to each other. Chlamydomonas perpusilla was shown to be sister to Chlorogonium Ehrenb. (clade Chlorogonia, within Caudivolvoxa), while C. neoplanoconvexa represented a second basal lineage within Caudivolvoxa, next to the clade Characiosiphonia. Although the morphology of C. neoplanoconvexa was not particularly outstanding, its unique phylogenetic position will encourage further investigation of this species and its uncultivated relatives. © 2011 Japanese Society of Phycology.

Sulfite (SO2) plays an important role in flavour stability in alcoholic beverages, whereas hydrogen sulfide (H2S) has an undesirable aroma. To discover the cellular processes that control SO2 and H2S production, we screened a library of Saccharomyces cerevisiae deletion mutants. Deletion of 12 genes led to increased H2S productivity. Ten of these genes are known to be involved in sulfur-containing amino acid metabolism, whereas UBI4 functions in the ubiquitin-proteasome system and SKP2 encodes an F-box-containing protein whose function is unknown. We found that the skp2 mutant accumulated H2S and SO2, because the adenosylphophosulfate kinase Met14p is a substrate of SCFSkp2 and more stable in the skp2 mutant than in the wild-type strain. Furthermore, the skp2 mutant grew more slowly than the wild-type strain under nutrient-limited conditions. Metabolome analysis showed that the concentration of intracellular cysteine is lower in the skp2 mutant than in the wild-type strain. The slow growth of the skp2 mutant was due to a lower concentration of intracellular cysteine, because the addition of cysteine suppressed the slow growth. In the skp2 mutant, the cysteine biosynthesis proteins Str2p, Str3p and Str4p are more stable than in the wild-type strain. Moreover, supplementation with methionine, S-adenosylmethionine, S-adenosylhomocysteine and homocysteine also suppressed the slow growth. Overexpression of STR1 or STR4 caused a more severe defect in the skp2 mutant. These results suggest that the balance of methionine and cysteine biosynthesis is important for yeast cell growth. Thus, Skp2p is one of the key components regulating this balance and H2S/SO2 production. © 2010 John Wiley & Sons, Ltd.

Salivary glands are intimately connected with the entire body via the nervous system and blood vessels, yet few studies have documented the effects of the salivary glands on the health of the whole body, although saliva has often been reported to possibly contain promising biomarkers of systemic disease states, including cancer. We therefore propose that the investigation of the effects of the salivary gland on the whole body and vice versa could introduce a new field of study, i.e., “salivary gland and health medicine.” This review examines data from novel studies on “salivary gland and health medicine” that were carried out by 3 different research groups. © 2011, Japanese Association for Oral Biology. All rights reserved.