Search

Summary Bacterial and viral pathogens affect the homeostasis of their hosts in complex and oftentimes unintuitive ways. Mechanistic insights into these host–pathogen interactions may one day allow us to discover more effective antimicrobials and antiviral therapies. The mechanisms that affect overall infection outcomes in the tissue and organ scales are often dictated by behaviors at the cellular and subcellular scales, such as cell metabolism, reproduction, aggregation, and survival. Computational modeling allows us to determine these large-scale outcomes by modeling smaller-scale phenomena, and this field of “multiscale computational modeling” has recently seen tremendous growth and innovation. In this article, we review recent studies of multiscale computational modeling of infection and their impact on the mechanistic understanding of host–pathogen interactions toward developing better therapies and multiscale modeling frameworks.

Building structural models of entire cells has been a long-standing cross-discipline challenge for the research community, as it requires an unprecedented level of integration between multiple sources of biological data and enhanced methods for computational modeling and visualization. Here, we present the first 3D structural models of an entire Mycoplasma genitalium (MG) cell, built using the CellPACK suite of computational modeling tools. Our model recapitulates the data described in recent whole-cell system biology simulations and provides a structural representation for all MG proteins, DNA and RNA molecules, obtained by combining experimental and homology-modeled structures and lattice-based models of the genome. We establish a framework for gathering, curating and evaluating these structures, exposing current weaknesses of modeling methods and the boundaries of MG structural knowledge, and visualization methods to explore functional characteristics of the genome and proteome. We compare two approaches for data gathering, a manually-curated workflow and an automated workflow that uses homologous structures, both of which are appropriate for the analysis of mesoscale properties such as crowding and volume occupancy. Analysis of model quality provides estimates of the regularization that will be required when these models are used as starting points for atomic molecular dynamics simulations.

Our ability to build computational models that account for all known gene functions in a cell has increased dramatically. But why build whole-cell models, and how can they best be used? In this review we enumerate several areas in which whole-cell modeling can significantly impact research and technology.

SummaryIncreasing evidence has shown that population dynamics are qualitatively different from single-cell behaviors. Reporters to probe dynamic, single-cell behaviors are desirable yet relatively scarce. Here, we describe an easy-to-implement and generalizable technology to generate reporters of kinase activity for individual cells. Our technology converts phosphorylation into a nucleocytoplasmic shuttling event that can be measured by epifluorescence microscopy. Our reporters reproduce kinase activity for multiple types of kinases and allow for calculation of active kinase concentrations via a mathematical model. Using this technology, we made several experimental observations that had previously been technicallyunfeasible, including stimulus-dependent patterns of c-Jun N-terminal kinase (JNK) and nuclear factor kappa B (NF-κB) activation. We also measured JNK, p38, and ERK activities simultaneously, finding that p38 regulates the peak number, but not the intensity, of ERK fluctuations. Our approach opens the possibility of analyzing a wide range of kinase-mediated processes in individual cells.

We present two modifications of the flux balance analysis (FBA) metabolic modeling framework which relax implicit assumptions of the biomass reaction. Our flexible flux balance analysis (flexFBA) objective removes the fixed proportion between reactants, and can therefore produce a subset of biomass reactants. Our time-linked flux balance analysis (tFBA) simulation removes the fixed proportion between reactants and byproducts, and can therefore describe transitions between metabolic steady states. Used together, flexFBA and tFBA model a time scale shorter than the regulatory and growth steady state encoded by the biomass reaction. This combined short-time FBA method is intended for integrated modeling applications to enable detailed and dynamic depictions of microbial physiology such as whole-cell modeling. For example, when modeling Escherichia coli, it avoids artifacts caused by low-copy-number enzymes in single-cell models with kinetic bounds. Even outside integrated modeling contexts, the detailed predictions of flexFBA and tFBA complement existing FBA techniques. We show detailed metabolite production of in silico knockouts used to identify when correct essentiality predictions are made for the wrong reason.

Constraints-based models have been effectively used to analyse, interpret, and predict the function of reconstructed genome-scale metabolic models. The first generation of these models used "hard" non-adjustable constraints associated with network connectivity, irreversibility of metabolic reactions, and maximal flux capacities. These constraints restrict the allowable behaviors of a network to a convex mathematical solution space whose edges are extreme pathways that can be used to characterize the optimal performance of a network under a stated performance criterion. The development of a second generation of constraints-based models by incorporating constraints associated with regulation of gene expression was described in a companion paper published in this journal, using flux-balance analysis to generate time courses of growth and by-product secretion using a skeleton representation of core metabolism. The imposition of these additional restrictions prevents the use of a subset of the extreme pathways that are derived from the "hard" constraints, thus reducing the solution space and restricting allowable network functions. Here, we examine the reduction of the solution space due to regulatory constraints using extreme pathway analysis. The imposition of environmental conditions and regulatory mechanisms sharply reduces the number of active extreme pathways. This approach is demonstrated for the skeleton system mentioned above, which has 80 extreme pathways. As regulatory constraints are applied to the system, the number of feasible extreme pathways is reduced to between 26 and 2 extreme pathways, a reduction of between 67.5 and 97.5%. The method developed here provides a way to interpret how regulatory mechanisms are used to constrain network functions and produce a small range of physiologically meaningful behaviors from all allowable network functions.

Full genome sequences enable the construction of genome-scale in silico models of complex cellular functions. Genome-scale constraints-based models of Escherichia coli metabolism have been constructed and used to successfully interpret and predict cellular behavior under a range of conditions. These previous models do not account for regulation of gene transcription and thus cannot accurately predict some organism functions. Here we present an in silico model of the central E. coli metabolism that accounts for regulation of gene expression. This model accounts for 149 genes, the products of which include 16 regulatory proteins and 73 enzymes. These enzymes catalyze 113 reactions, 45 of which are controlled by transcriptional regulation. The combined metabolic/regulatory model can predict the ability of mutant E. coli strains to grow on defined media as well as time courses of cell growth, substrate uptake, metabolic by-product secretion, and qualitative gene expression under various conditions, as indicated by comparison with experimental data under a variety of environmental conditions. The in silico model may also be used to interpret dynamic behaviors observed in cell cultures. This combined metabolic/regulatory model is thus an important step toward the goal of synthesizing genome-scale models that accurately represent E. coli behavior.

Genome-scale metabolic networks can now be reconstructed based on annotated genomic data augmented with biochemical and physiological information about the organism. Mathematical analysis can be performed to assess the capabilities of these reconstructed networks. The constraints-based framework, with #ux balance analysis (FBA), has been used successfully to predict time course of growth and by-product secretion, e!ects of mutation and knock-outs, and gene expression pro"les. However, FBA leads to incorrect predictions in situations where regulatory e!ects are a dominant in#uence on the behavior of the organism. Thus, there is a need to include regulatory events within FBA to broaden its scope and predictive capabilities. Here we represent transcriptional regulatory events as time-dependent constraints on the capabilities of a reconstructed metabolic network to further constrain the space of possible network functions. Using a simpli"ed metabolic/regulatory network, growth is simulated under various conditions to illustrate systemic e!ects such as catabolite repression, the aerobic/anaerobic diauxic shift and amino acid biosynthesis pathway repression. The incorporation of transcriptional regulatory events in FBA enables us to interpret, analyse and predict the e!ects of transcriptional regulation on cellular metabolism at the systemic level. ( 2001 Academic Press

The large volume of genome-scale data that is being produced and made available in databases on the World Wide Web is demanding the development of integrated mathematical models of cellular processes. The analysis of reconstructed metabolic networks as systems leads to the development of an in silico or computer representation of collections of cellular metabolic constituents, their interactions and their integrated function as a whole. The use of quantitative analysis methods to generate testable hypotheses and drive experimentation at a whole-genome level signals the advent of a systemic modeling approach to cellular and molecular biology.

A central challenge in biology is to understand how cellular life emerges from individual biochemical interactions. To address this challenge we have developed a novel computational framework which facilitates the integration of multiple disparate biochemical networks and data into a single unified model. Using this framework we have developed a detailed computational model of the complete life cycle of the smallest known freely-living organism, Mycoplasma genitalium. The model describes the life cycle of a single cell including DNA, RNA, and protein synthesis, metabolism, and cell division. The model accounts for the specific function of every annotated gene product, and simulates the dynamics of every molecular species. We have validated the model using several publically available experimental datasets. Currently we are using the model to gain insight into the control and regulation of cellular growth by exploring the effects of genomic and environmental perturbations on cellular behavior. In addition, we are developing an open-source, open-access web-based platform to facilitate broad application of our model.