Search

Growth is a fundamental process of life. Growth requirements are well-characterized experimentally for many microbes; however, we lack a unified model for cellular growth. Such a model must be predictive of events at the molecular scale and capable of explaining the high-level behavior of the cell as a whole. Here, we construct an ME-Model for Escherichia coli--a genome-scale model that seamlessly integrates metabolic and gene product expression pathways. The model computes ~80% of the functional proteome (by mass), which is used by the cell to support growth under a given condition. Metabolism and gene expression are interdependent processes that affect and constrain each other. We formalize these constraints and apply the principle of growth optimization to enable the accurate prediction of multi-scale phenotypes, ranging from coarse-grained (growth rate, nutrient uptake, by-product secretion) to fine-grained (metabolic fluxes, gene expression levels). Our results unify many existing principles developed to describe bacterial growth.

Background: Network Component Analysis (NCA) has been used to deduce the activities of transcription factors (TFs) from gene expression data and the TF-gene binding relationship. However, the TF-gene interaction varies in different environmental conditions and tissues, but such information is rarely available and cannot be predicted simply by motif analysis. Thus, it is beneficial to identify key TF-gene interactions under the experimental condition based on transcriptome data. Such information would be useful in identifying key regulatory pathways and gene markers of TFs in further studies., Results: We developed an algorithm to trim network connectivity such that the important regulatory interactions between the TFs and the genes were retained and the regulatory signals were deduced. Theoretical studies demonstrated that the regulatory signals were accurately reconstructed even in the case where only three independent transcriptome datasets were available. At least 80% of the main target genes were correctly predicted in the extreme condition of high noise level and small number of datasets. Our algorithm was tested with transcriptome data taken from mice under rapamycin treatment. The initial network topology from the literature contains 70 TFs, 778 genes, and 1423 edges between the TFs and genes. Our method retained 1074 edges (i.e. 75% of the original edge number) and identified 17 TFs as being significantly perturbed under the experimental condition. Twelve of these TFs are involved in MAPK signaling or myeloid leukemia pathways defined in the KEGG database, or are known to physically interact with each other. Additionally, four of these TFs, which are Hif1a, Cebpb, Nfkb1, and Atf1, are known targets of rapamycin. Furthermore, the trimmed network was able to predict Eno1 as an important target of Hif1a; this key interaction could not be detected without trimming the regulatory network., Conclusions: The advantage of our new algorithm, relative to the original NCA, is that our algorithm can identify the important TF-gene interactions. Identifying the important TF-gene interactions is crucial for understanding the roles of pleiotropic global regulators, such as p53. Also, our algorithm has been developed to overcome NCA's inability to analyze large networks where multiple TFs regulate a single gene. Thus, our algorithm extends the applicability of NCA to the realm of mammalian regulatory network analysis.

Biological signalling networks allow living organisms to issue an integrated response to current conditions and make limited predictions about future environmental changes. Small-scale dynamic models of signalling cascades, including mitogen-activated protein kinase cascades, have been developed to generate hypotheses about signal transduction. Owing to technical limitations, these models and the hypotheses they generate have focused on a limited subset of signalling molecules. Now that we can simultaneously measure a substantial portion of the molecular components of a cell, we can begin to develop and test systems-level models of cellular signalling and regulatory processes, therefore gaining insights into the 'thought' processes of a cell.

Nitric oxide (NO) is used by mammalian immune systems to counter microbial invasions and is produced by bacteria during denitrification. As a defense, microorganisms possess a complex network to cope with NO. Here we report a combined transcriptomic, chemical, and phenotypic approach to identify direct NO targets and construct the biochemical response network. In particular, network component analysis was used to identify transcription factors that are perturbed by NO. Such information was screened with potential NO reaction mechanisms and phenotypic data from genetic knockouts to identify active chemistry and direct NO targets in Escherichia coli. This approach identified the comprehensive E. coli NO response network and evinced that NO halts bacterial growth via inhibition of the branched-chain amino acid biosynthesis enzyme dihydroxyacid dehydratase. Because mammals do not synthesize branched-chain amino acids, inhibition of dihydroxyacid dehydratase may have served to foster the role of NO in the immune arsenal.