Module-based analysis of robustness tradeoffs in the heat shock response system.

Biological systems have evolved complex regulatory mechanisms, even in situations where much simpler designs seem to be sufficient for generating nominal functionality. Using module-based analysis coupled with rigorous mathematical comparisons, we propose that in analogy to control engineering architectures, the complexity of cellular systems and the presence of hierarchical modular structures can be attributed to the necessity of achieving robustness. We employ the Escherichia coli heat shock response system, a strongly conserved cellular mechanism, as an example to explore the design principles of such modular architectures. In the heat shock response system, the sigma-factor sigma32 is a central regulator that integrates multiple feedforward and feedback modules. Each of these modules provides a different type of robustness with its inherent tradeoffs in terms of transient response and efficiency. We demonstrate how the overall architecture of the system balances such tradeoffs. An extensive mathematical exploration nevertheless points to the existence of an array of alternative strategies for the existing heat shock response that could exhibit similar behavior. We therefore deduce that the evolutionary constraints facing the system might have steered its architecture toward one of many robustly functional solutions.

CADLIVE dynamic simulator: direct link of biochemical networks to dynamic models.

We have developed the CADLIVE (Computer-Aided Design of LIVing systEms) Simulator that provided a rule-based automatic way to convert biochemical network maps into dynamic models, which enables simulating their dynamics without going through all of the reactions down to the details of exact kinetic parameters. The simulator supports the biochemical reaction maps that are generated by the previously developed GUI editor. Notice that the part of the GUI editor had been previously published, but, as yet, not the simulator. To directly link biochemical network maps to dynamic simulation, we have created the strategy of three layers and two stages with the efficient conversion rules in an XML representation. This strategy divides a molecular network into three layers, i.e., gene, protein, and metabolic layers, and partitions the conversion process into two stages. Once a biochemical map is provided, CADLIVE automatically builds a mathematical model, thereby facilitating one to simulate and analyze it. In order to demonstrate the feasibility of CADLIVE, we analyzed the Escherichia coli nitrogen-assimilation system (64 equations with 64 variables) that consists of multiple and complicated negative and positive feedback loops. CADLIVE predicted that the glnK gene is responsible for hysteresis or reversibility of nitrogen-related (Ntr) gene expression with respect to the ammonia concentration, supporting the experimental observation of the runaway expression of the Ntr genes.