Eivind Almaas

Global organization of metabolic fluxes in the bacterium Escherichia coli

Cellular metabolism, the integrated interconversion of thousands of metabolic substrates through enzyme-catalysed biochemical reactions, is the most investigated complex intracellular web of molecular interactions. Although the topological organization of individual reactions into metabolic networks is well understood, the principles that govern their global functional use under different growth conditions raise many unanswered questions. By implementing a flux balance analysis of the metabolism of Escherichia coli strain MG1655, here we show that network use is highly uneven. Whereas most metabolic reactions have low fluxes, the overall activity of the metabolism is dominated by several reactions with very high fluxes. E. coli responds to changes in growth conditions by reorganizing the rates of selected fluxes predominantly within this high-flux backbone. This behaviour probably represents a universal feature of metabolic activity in all cells, with potential implications for metabolic engineering.

Comparative Genome-Scale Metabolic Reconstruction and Flux Balance Analysis of Multiple Staphylococcus aureus Genomes Identify Novel Antimicrobial Drug Targets

Mortality due to multidrug-resistant Staphylococcus aureus infection is predicted to surpass that of human immunodeficiency virus/AIDS in the United States. Despite the various treatment options for S. aureus infections, it remains a major hospital- and community-acquired opportunistic pathogen. With the emergence of multidrug-resistant S. aureus strains, there is an urgent need for the discovery of new antimicrobial drug targets in the organism. To this end, we reconstructed the metabolic networks of multidrug-resistant S. aureus strains using genome annotation, functional-pathway analysis, and comparative genomic approaches, followed by flux balance analysis-based in silico single and double gene deletion experiments. We identified 70 single enzymes and 54 pairs of enzymes whose corresponding metabolic reactions are predicted to be unconditionally essential for growth. Of these, 44 single enzymes and 10 enzyme pairs proved to be common to all 13 S. aureus strains, including many that had not been previously identified as being essential for growth by gene deletion experiments in S. aureus. We thus conclude that metabolic reconstruction and in silico analyses of multiple strains of the same bacterial species provide a novel approach for potential antibiotic target identification.

The Activity Reaction Core and Plasticity of Metabolic Networks

Understanding the system-level adaptive changes taking place in an organism in response to variations in the environment is a key issue of contemporary biology. Current modeling approaches, such as constraint-based flux-balance analysis, have proved highly successful in analyzing the capabilities of cellular metabolism, including its capacity to predict deletion phenotypes, the ability to calculate the relative flux values of metabolic reactions, and the capability to identify properties of optimal growth states. Here, we use flux-balance analysis to thoroughly assess the activity of Escherichia coli, Helicobacter pylori, and Saccharomyces cerevisiae metabolism in 30,000 diverse simulated environments. We identify a set of metabolic reactions forming a connected metabolic core that carry non-zero fluxes under all growth conditions, and whose flux variations are highly correlated. Furthermore, we find that the enzymes catalyzing the core reactions display a considerably higher fraction of phenotypic essentiality and evolutionary conservation than those catalyzing noncore reactions. Cellular metabolism is characterized by a large number of species-specific conditionally active reactions organized around an evolutionary conserved, but always active, metabolic core. Finally, we find that most current antibiotics interfering with bacterial metabolism target the core enzymes, indicating that our findings may have important implications for antimicrobial drug-target discovery.

Predicting synthetic rescues in metabolic networks

An important goal of medical research is to develop methods to recover the loss of cellular function due to mutations and other defects. Many approaches based on gene therapy aim to repair the defective gene or to insert genes with compensatory function. Here, we propose an alternative, network‐based strategy that aims to restore biological function by forcing the cell to either bypass the functions affected by the defective gene, or to compensate for the lost function. Focusing on the metabolism of single‐cell organisms, we computationally study mutants that lack an essential enzyme, and thus are unable to grow or have a significantly reduced growth rate. We show that several of these mutants can be turned into viable organisms through additional gene deletions that restore their growth rate. In a rather counterintuitive fashion, this is achieved via additional damage to the metabolic network. Using flux balance‐based approaches, we identify a number of synthetically viable gene pairs, in which the removal of one enzyme‐encoding gene results in a non‐viable phenotype, while the deletion of a second enzyme‐encoding gene rescues the organism. The systematic network‐based identification of compensatory rescue effects may open new avenues for genetic interventions.

Differences in human and chimpanzee gene expression patterns define an evolving network of transcription factors in brain

Humans differ from other primates by marked differences in cognitive abilities and a significantly larger brain. These differences correlate with metabolic changes, as evidenced by the relative up-regulation of energy-related genes and metabolites in human brain. While the mechanisms underlying these evolutionary changes have not been elucidated, altered activities of key transcription factors (TFs) could play a pivotal role. To assess this possibility, we analyzed microarray data from five tissues from humans and chimpanzees. We identified 90 TF genes with significantly different expression levels in human and chimpanzee brain among which the rapidly evolving KRAB-zinc finger genes are markedly over-represented. The differentially expressed TFs cluster within a robust regulatory network consisting of two distinct but interlinked modules, one strongly associated with energy metabolism functions, and the other with transcription, vesicular transport, and ubiquitination. Our results suggest that concerted changes in a relatively small number of interacting TFs may coordinate major gene expression differences in human and chimpanzee brain.

Evolutionary cores of domain co-occurrence networks

BackgroundThe modeling of complex systems, as disparate as the World Wide Web and the cellular metabolism, as networks has recently uncovered a set of generic organizing principles: Most of these systems are scale-free while at the same time modular, resulting in a hierarchical architecture. The structure of the protein domain network, where individual domains correspond to nodes and their co-occurrences in a protein are interpreted as links, also falls into this category, suggesting that domains involved in the maintenance of increasingly developed, multicellular organisms accumulate links. Here, we take the next step by studying link based properties of the protein domain co-occurrence networks of the eukaryotes S. cerevisiae, C. elegans, D. melanogaster, M. musculus and H. sapiens.ResultsWe construct the protein domain co-occurrence networks from the PFAM database and analyze them by applying a k-core decomposition method that isolates the globally central (highly connected domains in the central cores) from the locally central (highly connected domains in the peripheral cores) protein domains through an iterative peeling process. Furthermore, we compare the subnetworks thus obtained to the physical domain interaction network of S. cerevisiae. We find that the innermost cores of the domain co-occurrence networks gradually grow with increasing degree of evolutionary development in going from single cellular to multicellular eukaryotes. The comparison of the cores across all the organisms under consideration uncovers patterns of domain combinations that are predominately involved in protein functions such as cell-cell contacts and signal transduction. Analyzing a weighted interaction network of PFAM domains of Yeast, we find that domains having only a few partners frequently interact with these, while the converse is true for domains with a multitude of partners. Combining domain co-occurrence and interaction information, we observe that the co-occurrence of domains in the innermost cores (globally central domains) strongly coincides with physical interaction. The comparison of the multicellular eukaryotic domain co-occurrence networks with the single celled of S. cerevisiae (the overlap network) uncovers small, connected network patterns.ConclusionWe hypothesize that these patterns, consisting of the domains and links preserved through evolution, may constitute nucleation kernels for the evolutionary increase in proteome complexity. Combining co-occurrence and physical interaction data we argue that the driving force behind domain fusions is a collective effect caused by the number of interactions and not the individual interaction frequency.

Radiation forces on a micrometer-sized sphere in an evanescent field

Electromagnetic wave theory is used to predict the radiation forces exerted upon a micrometer-sized spherical particle illuminated by evanescent waves penetrating across a dielectric interface. These forces are quantified for two incident beam polarizations (p and s polarization) and for different refractive-index media. The electromagnetic formalism that we use is based on theoretical studies of Barton et al. [ J. Appl. Phys.64, 1632 ( 1988);J. Appl. Phys.66, 4594 ( 1989)]. The novelty of the present research is to apply this general formalism to the calculation of forces when the evanescent field is identified with the incident field. Our theoretical results for the horizontal and vertical force components are shown graphically in nondimensional form as functions of the size parameter of the sphere. Moreover, our results are found to be in reasonable agreement with recent experimental findings of Kawata and Sugiura [ Opt. Lett.17, 772 ( 1992)].

Power laws, scale-free networks and genome biology

Power Laws in Biological Networks.- Graphical Analysis of Biocomplex Networks and Transport Phenomena.- Large-Scale Topological Properties of Molecular Networks.- The Connectivity of Large Genetic Networks.- The Drosophila Protein Interaction Network May Be neither Power-Law nor Scale-Free.- Birth and Death Models of Genome Evolution.- Scale-Free Evolution.- Gene Regulatory Networks.- Power Law Correlations in DNA Sequences.- Analytical Evolutionary Model for Protein Fold Occurrence in Genomes, Accounting for the Effects of Gene Duplication, Deletion, Acquisition and Selective Pressure.- The Protein Universes.- The Role of Computation in Complex Regulatory Networks.- Neutrality and Selection in the Evolution of Gene Families.- Scaling Laws in the Functional Content of Genomes.

Genome-level transcription data of Yersinia pestis analyzed with a New metabolic constraint-based approach

BackgroundConstraint-based computational approaches, such as flux balance analysis (FBA), have proven successful in modeling genome-level metabolic behavior for conditions where a set of simple cellular objectives can be clearly articulated. Recently, the necessity to expand the current range of constraint-based methods to incorporate high-throughput experimental data has been acknowledged by the proposal of several methods. However, these methods have rarely been used to address cellular metabolic responses to some relevant perturbations such as antimicrobial or temperature-induced stress. Here, we present a new method for combining gene-expression data with FBA (GX-FBA) that allows modeling of genome-level metabolic response to a broad range of environmental perturbations within a constraint-based framework. The method uses mRNA expression data to guide hierarchical regulation of cellular metabolism subject to the interconnectivity of the metabolic network.ResultsWe applied GX-FBA to a genome-scale model of metabolism in the gram negative bacterium Yersinia pestis and analyzed its metabolic response to (i) variations in temperature known to induce virulence, and (ii) antibiotic stress. Without imposition of any a priori behavioral constraints, our results show strong agreement with reported phenotypes. Our analyses also lead to novel insights into how Y. pestis uses metabolic adjustments to counter different forms of stress.ConclusionsComparisons of GX-FBA predicted metabolic states with fluxomic measurements and different reported post-stress phenotypes suggest that mass conservation constraints and network connectivity can be an effective representative of metabolic flux regulation in constraint-based models. We believe that our approach will be of aid in the in silico evaluation of cellular goals under different conditions and can be used for a variety of analyses such as identification of potential drug targets and their action.

Trend Motif: A Graph Mining Approach for Analysis of Dynamic Complex Networks

Complex networks have been used successfully in scientific disciplines ranging from sociology to microbiology to describe systems of interacting units. Until recently, studies of complex networks have mainly focused on their network topology. However, in many real world applications, the edges and vertices have associated attributes that are frequently represented as vertex or edge weights. Furthermore, these weights are often not static, instead changing with time and forming a time series. Hence, to fully understand the dynamics of the complex network, we have to consider both network topology and related time series data. In this work, we propose a motif mining approach to identify trend motifs for such purposes. Simply stated, a trend motif describes a recurring subgraph where each of its vertices or edges displays similar dynamics over a user- defined period. Given this, each trend motif occurrence can help reveal significant events in a complex system; frequent trend motifs may aid in uncovering dynamic rules of change for the system, and the distribution of trend motifs may characterize the global dynamics of the system. Here, we have developed efficient mining algorithms to extract trend motifs. Our experimental validation using three disparate empirical datasets, ranging from the stock market, world trade, to a protein interaction network, has demonstrated the efficiency and effectiveness of our approach.