Modeling biogeochemical dynamics in porous media: Practical considerations of pore scale variability, reaction networks, and microbial population dynamics in a sandy aquifer.

Prediction of the fate and environmental impacts of groundwater contaminants requires the identification of relevant biogeochemical processes and necessitates the macroscopic representation of microbial activity occurring at the microscale. Using a well-studied sandy aquifer environment, we evaluate the importance of pore distribution on organic matter respiration in a porous medium environment by performing spatially explicit simulations of microbial metabolism at the sub-millimeter scale. Model results using an idealized porous medium under non-biofilm forming conditions indicate that while some heterogeneity is observed for flow rates, distributions of microbes and dissolved organic substrates remain relatively homogenous at the grain scale. At the macroscale in the same environment, we assess the impact of a comprehensive reaction network description for a phenolic contaminant plume, and compare the findings to a setting describing organic matter breakdown in a coastal marine sediment. This comparison reveals the importance of reactions recycling reduced metabolites at redox interfaces, leading to a competition for oxidants. When the spatio-temporal dynamics of microbial groups are accounted for, our simulations show the importance of reaction energetics and nutrient limitations such as microbial nitrogen demands.

In silico Geobacter sulfurreducens metabolism and its representation in reactive transport models.

Microbial activity governs elemental cycling and the transformation of many anthropogenic substances in aqueous environments. Through the development of a dynamic cell model of the well-characterized, versatile, and abundant Geobacter sulfurreducens, we showed that a kinetic representation of key components of cell metabolism matched microbial growth dynamics observed in chemostat experiments under various environmental conditions and led to results similar to those from a comprehensive flux balance model. Coupling the kinetic cell model to its environment by expressing substrate uptake rates depending on intra- and extracellular substrate concentrations, two-dimensional reactive transport simulations of an aquifer were performed. They illustrated that a proper representation of growth efficiency as a function of substrate availability is a determining factor for the spatial distribution of microbial populations in a porous medium. It was shown that simplified model representations of microbial dynamics in the subsurface that only depended on extracellular conditions could be derived by properly parameterizing emerging properties of the kinetic cell model.

Chromosome segregation in Escherichia coli division: a free energy-driven string model.

Although the mechanisms of eukaryotic chromosome segregation and cell division have been elucidated to a certain extent, those for bacteria remain largely unknown. Here we present a computational string model for simulating the dynamics of Escherichia coli chromosome segregation. A novel thermal-average force field accounting for stretching, bending, volume exclusion, friction and random fluctuation is introduced. A Langevin equation is used to simulate the chromosome structural changes. The mechanism of chromosome segregation is thereby postulated as a result of free energy-driven structural optimization with replication introduced chromosomal mass increase. Predictions of the model agree well with observations of fluorescence labeled chromosome loci movement in living cells. The results demonstrate the possibility of a mechanism of chromosome segregation that does not involve cytoskeletal guidance or advanced apparatus in an E. coli cell. The model also shows that DNA condensation of locally compacted domains is a requirement for successful chromosome segregation. Simulations also imply that the shape-determining protein MreB may play a role in the segregation via modification of the membrane pressure.

Cancer onset and progression: a genome-wide, nonlinear dynamical systems perspective on onconetworks.

It is hypothesized that the many human cell types corresponding to multiple states is supported by an underlying nonlinear dynamical system (NDS) of transcriptional regulatory network (TRN) processes. This hypothesis is validated for epithelial cells whose TRN is found to support an extremely complex array of states that we term a "bifurcation nexus", for which we introduce a quantitative measure of complexity. The TRN used is constructed and analyzed by integrating a database of TRN information, cDNA microarray data analyzers, bioinformatics modules, a transcription/translation/post-translation kinetic model, and NDS analysis software. Results of this genome-wide approach suggest that a cell can be induced to persist in one state or to transition between distinct states; apparently irreversible transitions can be reversed when the high dimensional space of extracellular and intracellular parameters is understood. As conditions change, certain cellular states (cell lines) are no longer supported, new ones emerge, and transitions (cell differentiation or death) occur. The accumulation of simulated point mutations (minor changes which individually are insignificant) lead to occasional dramatic transitions. The genome-wide scope of many of these transitions is shown to arise from the cross-linked TRN structure. These notions imply that studying individual oncogenes may not be sufficient to understand cancer; rather, "onconetworks" (subsets of strongly coupled genes supporting multiple cell states) should be considered. Our approach reveals several epithelial onconetworks, each involving oncogenes and anti-tumor and supporting genes.

The Karyote physico-chemical genomic, proteomic, metabolic cell modeling system.

Modeling approaches to the dynamics of a living cell are presented that are strongly based on its underlying physical and chemical processes and its hierarchical spatio-temporal organization. Through the inclusion of a broad spectrum of processes and a rigorous analysis of the multiple scale nature of cellular dynamics, we are attempting to advance cell modeling and its applications. The presentation focuses on our cell modeling system, which integrates data archiving and quantitative physico-chemical modeling and information theory to provide a seamless approach to the modeling/data analysis endeavor. Thereby the rapidly growing mess of genomic, proteomic, metabolic, and cell physiological data can be automatically used to develop and calibrate a predictive cell model. The discussion focuses on the Karyote cell modeling system and an introduction to the CellX and VirusX models. The Karyote software system integrates three elements: (1) a model-building and data archiving module that allows one to define a cell type to be modeled through its reaction network, structure, and transport processes as well as to choose the surrounding medium and other parameters of the phenomenon to be modeled; (2) a genomic, proteomic, metabolic cell simulator that solves the equations of metabolic reaction, transcription/translation polymerization and the exchange of molecules between parts of the cell and with the surrounding medium; and (3) an information theory module (ITM) that automates model calibration and development, and integrates a variety of data types with the cell dynamic computations. In Karyote, reactions may be fast (equilibrated) or slow (finite rate), and the special effects of enzymes and other minority species yielding steady-state cycles of arbitrary complexities are accounted for. These features of the dynamics are handled via rigorous multiple scale analysis. A user interface allows for an automated generation and solution of the equations of multiple timescale, compartmented dynamics. Karyote is based on a fixed intracellular structure. However, cell response to changes in the host medium, damage, development or transformation to abnormality can involve dramatic changes in intracellular structure. As this changes the nature of the cellular dynamics, a new model, CellX, is being developed based on the spatial distribution of concentration and other variables. This allows CellX to capture the self-organizing character of cellular behavior. The self-assembly of organelles, viruses, and other subcellular bodies is being addressed in a second new model, VirusX, that integrates molecular mechanics and continuum theory. VirusX is designed to study the influence of a host medium on viral self-assembly, structural stability, infection of a single cell, and transmission of disease.

Toward Automated Cell Model Development through Information Theory†.

The objective of this paper is to present a methodology for developing and calibrating models of complex reaction/transport systems. In particular, the complex network of biochemical reaction/transport processes and their spatial organization make the development of a predictive model of a living cell a grand challenge for the 21st century. However, advances in reaction/transport modeling and the exponentially growing databases of genomic, proteomic, metabolic, and bioelectric data make cell modeling feasible, if these two elements can be automatically integrated in an unbiased fashion. In this paper, we present a procedure to integrate data with a new cell model, Karyote, that accounts for many of the physical processes needed to attain the goal of predictive modeling. Our integration methodology is based on the use of information theory. The model is integrated with a variety of types and qualities of experimental data using an objective error assessment approach. Data that can be used in this approach include NMR, spectroscopy, microscopy, and electric potentiometry. The approach is demonstrated on the well-studied Trypanosoma brucei system. A major obstacle for the development of a predictive cell model is that the complexity of these systems makes it unlikely that any model presently available will soon be complete in terms of the set of processes accounted for. Thus, one is faced with the challenge of calibrating and running an incomplete model. We present a probability functional method that allows the integration of experimental data and soft information such as choice of error measure, a priori information, and physically motivated regularization to address the incompleteness challenge.