The compositional and evolutionary logic of metabolism

Metabolism displays striking and robust regularities in the forms of modularity and hierarchy, whose composition may be compactly described. This renders metabolic architecture comprehensible as a system, and suggests the order in which layers of that system emerged. Metabolism also serves as the foundation in other hierarchies, at least up to cellular integration including bioenergetics and molecular replication, and trophic ecology. The recapitulation of patterns first seen in metabolism, in these higher levels, suggests metabolism as a source of causation or constraint on many forms of organization in the biosphere. We identify as modules widely reused subsets of chemicals, reactions, or functions, each with a conserved internal structure. At the small molecule substrate level, module boundaries are generally associated with the most complex reaction mechanisms and the most conserved enzymes. Cofactors form a structurally and functionally distinctive control layer over the small-molecule substrate. Complex cofactors are often used at module boundaries of the substrate level, while simpler ones participate in widely used reactions. Cofactor functions thus act as "keys" that incorporate classes of organic reactions within biochemistry. The same modules that organize the compositional diversity of metabolism are argued to have governed long-term evolution. Early evolution of core metabolism, especially carbon-fixation, appears to have required few innovations among a small number of conserved modules, to produce adaptations to simple biogeochemical changes of environment. We demonstrate these features of metabolism at several levels of hierarchy, beginning with the small-molecule substrate and network architecture, continuing with cofactors and key conserved reactions, and culminating in the aggregation of multiple diverse physical and biochemical processes in cells.Comment: 56 pages, 28 figure

Computational analysis of innate and adaptive immune responses

Both innate and adaptive immune processes rely on the activation of differentiated haematopoietic stem cell lineages to affect an appropriate response to pathogens. This thesis employs a largely network biology focused approach to better understand the specificity of immune cell responses in two distinct cases of pathogenic challenge. In the context of adaptive immunity, I studied the transcriptional responses of T cells during Graft-versus-Host Disease (GvHD). GvHD represents one of the major complications to arise following allogeneic hematopoietic stem cell transplantation and yet why only particular organs are damaged as a result of this pathology is still unclear. To investigate whether key GvHD transcriptional signatures seen in effector CD8+ T cells compared to naïve T cells are triggered in target organs or the secondary lymphoid organs, a module-based association test was developed to combine the output of gene clustering algorithms with that of differential expression analysis. This methodology significantly aided the identification of skin specific effector T cell transcriptional programs believed to drive murine GvHD pathogenesis at this site. Turning to the innate immune response, I investigated the transcriptional profiles of resting and activated macrophages in the setting of Tuberculosis (TB), the second leading cause of death from infectious disease worldwide. Regression-based analyses and clustering of macrophage expression data provided insight into the variations in gene expression profiles of naïve macrophages compared to those infected with Mycobacterium tuberculosis (MTB) or a vaccine strain of mycobacteria (BCG). The availability of genotype data as part of the macrophage dataset facilitated an expression quantitative trait loci (eQTL) study which highlighted a novel association between the cytoskeleton gene BCAR1 and TB risk, together with a previously undescribed trans-eQTL module specific to MTB infected macrophages. Potential genetic variants impacting expression of the aforementioned GvHD specific T cell transcriptional signatures were additionally investigated using external trans-eQTL datasets

A systems biology understanding of protein constraints in the metabolism of budding yeasts

Fermentation technologies, such as bread making and production of alcoholic beverages, have been crucial for development of humanity throughout history. Saccharomyces cerevisiae provides a natural platform for this, due to its capability to transform sugars into ethanol. This, and other yeasts, are now used for production of pharmaceuticals, including insulin and artemisinic acid, flavors, fragrances, nutraceuticals, and fuel precursors. In this thesis, different systems biology methods were developed to study interactions between metabolism, enzymatic capabilities, and regulation of gene expression in budding yeasts. In paper I, a study of three different yeast species (S. cerevisiae, Yarrowia lipolytica and Kluyveromyces marxianus), exposed to multiple conditions, was carried out to understand their adaptation to environmental stress. Paper II revises the use of genome-scale metabolic models (GEMs) for the study and directed engineering of diverse yeast species. Additionally, 45 GEMs for different yeasts were collected, analyzed, and tested. In paper III, GECKO 2.0, a toolbox for integration of enzymatic constraints and proteomics data into GEMs, was developed and used for reconstruction of enzyme-constrained models (ecGEMs) for three yeast species and model organisms. Proteomics data and ecGEMs were used to further characterize the impact of environmental stress over metabolism of budding yeasts. On paper IV, gene engineering targets for increased accumulation of heme in S. cerevisiae cells were predicted with an ecGEM. Predictions were experimentally validated, yielding a 70-fold increase in intracellular heme. The prediction method was systematized and applied to the production of 102 chemicals in S. cerevisiae (Paper V). Results highlighted general principles for systems metabolic engineering and enabled understanding of the role of protein limitations in bio-based chemical production. Paper VI presents a hybrid model integrating an enzyme-constrained metabolic network, coupled to a gene regulatory model of nutrient-sensing mechanisms in S. cerevisiae. This model improves prediction of protein expression patterns while providing a rational connection between metabolism and the use of nutrients from the environment.This thesis demonstrates that integration of multiple systems biology approaches is valuable for understanding the connection of cell physiology at different levels, and provides tools for directed engineering of cells for the benefit of society

The functional network in predictive biology : predicting phenotype from genotype and predicting human disease from fungal phenotype

textThe ability to predict is one of the hallmarks of successful theories. Historically, the predictive power of biology has lagged behind disciplines like physics because the biological world is complex, challenging to quantify, and full of exceptions. However, in recent years the amount of available data has expanded exponentially and biological predictions based on this data become a possibility. The functional gene network is a quantitative way to integrate this data and a useful framework for making biological predictions. This study demonstrates that functional networks capture real biological insight and uses the network to predict both subcellular protein localization and the phenotypic outcome of gene knockouts. Furthermore, I use the functional network to evaluate genetic modules shared between diverse organisms that lead to orthologous phenotypes, many that are non-obvious. I show that the successful predictions of the functional network have broad applicability and implications that range from the design of large-scale biological experiments to the discovery of genes with potential roles in human disease.Institute for Cellular and Molecular Biolog

On the foundations of cancer modelling: selected topics, speculations, & perspectives

This paper presents a critical review of selected topics related to the modelling of cancer onset, evolution and growth, with the aim of illustrating, to a wide applied mathematical readership, some of the novel mathematical problems in the field. This review attempts to capture, from the appropriate literature, the main issues involved in the modelling of phenomena related to cancer dynamics at all scales which characterise this highly complex system: from the molecular scale up to that of tissue. The last part of the paper discusses the challenge of developing a mathematical biological theory of tumour onset and evolution

Integrating Symbolic and Neural Processing in a Self-Organizing Architechture for Pattern Recognition and Prediction

British Petroleum (89A-1204); Defense Advanced Research Projects Agency (N00014-92-J-4015); National Science Foundation (IRI-90-00530); Office of Naval Research (N00014-91-J-4100); Air Force Office of Scientific Research (F49620-92-J-0225

Finding the pathology of major depression through effects on gene interaction networks

The disease signature of major depressive disorder is distributed across multiple physical scales and investigative specialties, including genes, cells and brain regions. No single mechanism or pathway currently implicated in depression can reproduce its diverse clinical presentation, which compounds the difficulty in finding consistently disrupted molecular functions. We confront these key roadblocks to depression research - multi-scale and multi-factor pathology - by conducting parallel investigations at the levels of genes, neurons and brain regions, using transcriptome networks to identify collective patterns of dysfunction. Our findings highlight how the collusion of multi-system deficits can form a broad-based, yet variable pathology behind the depressed phenotype. For instance, in a variant of the classic lethality-centrality relationship, we show that in neuropsychiatric disorders including major depression, differentially expressed genes are pushed out to the periphery of gene networks. At the level of cellular function, we develop a molecular signature of depression based on cross-species analysis of human and mouse microarrays from depression-affected areas, and show that these genes form a tight module related to oligodendrocyte function and neuronal growth/structure. At the level of brain-region communication, we find a set of genes and hormones associated with the loss of feedback between the amygdala and anterior cingulate cortex, based on a novel assay of interregional expression synchronization termed "gene coordination". These results indicate that in the absence of a single pathology, depression may be created by dysynergistic effects among genes, cell-types and brain regions, in what we term the "floodgate" model of depression. Beyond our specific biological findings, these studies indicate that gene interaction networks are a coherent framework in which to understand the faint expression changes found in depression and complex neuropsychiatric disorders