An integrative and practical evolutionary optimization for a complex, dynamic model of biological networks

Development of an accurate kinetic model for the central carbon metabolism of Escherichia coli

Background

A kinetic model provides insights into the dynamic response of biological systems and predicts how their complex metabolic and gene regulatory networks generate particular functions. Of many biological systems, Escherichia coli metabolic pathways have been modeled extensively at the enzymatic and genetic levels, but existing models cannot accurately reproduce experimental behaviors in a batch culture, due to the inadequate estimation of a specific cell growth rate and a large number of unmeasured parameters.

Results

In this study, we developed a detailed kinetic model for the central carbon metabolism of E. coli in a batch culture, which includes the glycolytic pathway, tricarboxylic acid cycle, pentose phosphate pathway, Entner-Doudoroff pathway, anaplerotic pathway, glyoxylate shunt, oxidative phosphorylation, phosphotransferase system (Pts), non-Pts and metabolic gene regulations by four protein transcription factors: cAMP receptor, catabolite repressor/activator, pyruvate dehydrogenase complex repressor and isocitrate lyase regulator. The kinetic parameters were estimated by a constrained optimization method on a supercomputer. The model estimated a specific growth rate based on reaction kinetics and accurately reproduced the dynamics of wild-type E. coli and multiple genetic mutants in a batch culture.

Conclusions

This model overcame the intrinsic limitations of existing kinetic models in a batch culture, predicted the effects of multilayer regulations (allosteric effectors and gene expression) on central carbon metabolism and proposed rationally designed fast-growing cells based on understandings of molecular processes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-016-0511-x) contains supplementary material, which is available to authorized users.

BioFNet: biological functional network database for analysis and synthesis of biological systems

In synthetic biology and systems biology, a bottom-up approach can be used to construct a complex, modular, hierarchical structure of biological networks. To analyze or design such networks, it is critical to understand the relationship between network structure and function, the mechanism through which biological parts or biomolecules are assembled into building blocks or functional networks. A functional network is defined as a subnetwork of biomolecules that performs a particular function. Understanding the mechanism of building functional networks would help develop a methodology for analyzing the structure of large-scale networks and design a robust biological circuit to perform a target function. We propose a biological functional network database, named BioFNet, which can cover the whole cell at the level of molecular interactions. The BioFNet takes an advantage in implementing the simulation program for the mathematical models of the functional networks, visualizing the simulated results. It presents a sound basis for rational design of biochemical networks and for understanding how functional networks are assembled to create complex high-level functions, which would reveal design principles underlying molecular architectures.

Multi-objective optimal control of dynamic bioprocesses using ACADO Toolkit

The optimal design and operation of dynamic bioprocesses gives in practice often rise to optimisation problems with multiple and conflicting objectives. As a result typically not a single optimal solution but a set of Pareto optimal solutions exist. From this set of Pareto optimal solutions, one has to be chosen by the decision maker. Hence, efficient approaches are required for a fast and accurate generation of the Pareto set such that the decision maker can easily and systematically evaluate optimal alternatives. In the current paper the multi-objective optimisation of several dynamic bioprocess examples is performed using the freely available ACADO Multi-Objective Toolkit ( http://www.acadotoolkit.org ). This toolkit integrates efficient multiple objective scalarisation strategies (e.g., Normal Boundary Intersection and (Enhanced) Normalised Normal Constraint) with fast deterministic approaches for dynamic optimisation (e.g., single and multiple shooting). It has been found that the toolkit is able to efficiently and accurately produce the Pareto sets for all bioprocess examples. The resulting Pareto sets are added as supplementary material to this paper.