Dynamic optimization of biological networks under parametric uncertainty

Abstract

BackgroundMicro-organisms play an important role in various industrial sectors (including biochemical, food and pharmaceutical industries). A profound insight in the biochemical reactions inside micro-organisms enables an improved biochemical process control. Biological networks are an important tool in systems biology for incorporating microscopic level knowledge. Biochemical processes are typically dynamic and the cells have often more than one objective which are typically conflicting, e.g., minimizing the energy consumption while maximizing the production of a specific metabolite. Therefore multi-objective optimization is needed to compute trade-offs between those conflicting objectives. In model-based optimization, one of the inherent problems is the presence of uncertainty. In biological processes, this uncertainty can be present due to, e.g., inherent biological variability. Not taking this uncertainty into account, possibly leads to the violation of constraints and erroneous estimates of the actual objective function(s). To account for the variance in model predictions and compute a prediction interval, this uncertainty should be taken into account during process optimization. This leads to a challenging optimization problem under uncertainty, which requires a robustified solution.ResultsThree techniques for uncertainty propagation: linearization, sigma points and polynomial chaos expansion, are compared for the dynamic optimization of biological networks under parametric uncertainty. These approaches are compared in two case studies: (i) a three-step linear pathway model in which the accumulation of intermediate metabolites has to be minimized and (ii) a glycolysis inspired network model in which a multi-objective optimization problem is considered, being the minimization of the enzymatic cost and the minimization of the end time before reaching a minimum extracellular metabolite concentration. A Monte Carlo simulation procedure has been applied for the assessment of the constraint violations. For the multi-objective case study one Pareto point has been considered for the assessment of the constraint violations. However, this analysis can be performed for any Pareto point.ConclusionsThe different uncertainty propagation strategies each offer a robustified solution under parametric uncertainty. When making the trade-off between computation time and the robustness of the obtained profiles, the sigma points and polynomial chaos expansion strategies score better in reducing the percentage of constraint violations. This has been investigated for a normal and a uniform parametric uncertainty distribution. The polynomial chaos expansion approach allows to directly take prior knowledge of the parametric uncertainty distribution into account.Electronic supplementary materialThe online version of this article (doi:10.1186/s12918-016-0328-6) contains supplementary material, which is available to authorized users.

An Analysis of Uncertainty Propagation Methods Applied to Breakage Population Balance

“…In the polynomial chaos expansion method (PCE), the model response is approximated by an infinite series of orthogonal basis functions [13]. For practical applications, the infinite series is truncated to a limited number of terms M.…”

“…These methods can be a complex set of equations that need to be derived and solved for each case. Non-intrusive methods are based on sampling by repeated model evaluations at the collocation points [13,17]. The number of collocation points should be higher or equal to the number of coefficients in the PCE.…”

A simple and flexible computational framework for inferring sources of heterogeneity from single-cell dynamics

“…Simplicity and generality therefore make the GTS framework ideal for future expansions to address potential practical limitations. For example, one could use further parallelization and more efficient numerical routines 42 for the expensive model simulations, integrate global optimization routines in the first stage, use alternative likelihood functions 43 to guard against outlier-corrupted data, or use filtering methods 44 integrating a stochastic model of gene expression to identify the processes underlying phenotypic heterogeneity in more detail. Second, we expect the GTS framework to find applications also in systems pharmacology, specifically in combining systems biology models with preclinical data for translational medicine, where increasingly complex models meet data of increasing quantity and quality 45 .…”

Using optimal control to understand complex metabolic pathways

“…The work of Klipp et al [72] spurred the application of optimal control to characterize cellular dynamics related with metabolism [80][81][82][83][84][85][86][87][88][89][90][91][92][93][94][95][96]. Ewald et al [77] have recently reviewed many of these studies, illustrating how dynamic optimization is a powerful approach that can be used to decipher the activation and regulation of metabolism.…”