Dataset Information

A mathematical framework for yield (vs. rate) optimization in constraint-based modeling and applications in metabolic engineering.

ABSTRACT:

Background

The optimization of metabolic rates (as linear objective functions) represents the methodical core of flux-balance analysis techniques which have become a standard tool for the study of genome-scale metabolic models. Besides (growth and synthesis) rates, metabolic yields are key parameters for the characterization of biochemical transformation processes, especially in the context of biotechnological applications. However, yields are ratios of rates, and hence the optimization of yields (as nonlinear objective functions) under arbitrary linear constraints is not possible with current flux-balance analysis techniques. Despite the fundamental importance of yields in constraint-based modeling, a comprehensive mathematical framework for yield optimization is still missing.

Results

We present a mathematical theory that allows one to systematically compute and analyze yield-optimal solutions of metabolic models under arbitrary linear constraints. In particular, we formulate yield optimization as a linear-fractional program. For practical computations, we transform the linear-fractional yield optimization problem to a (higher-dimensional) linear problem. Its solutions determine the solutions of the original problem and can be used to predict yield-optimal flux distributions in genome-scale metabolic models. For the theoretical analysis, we consider the linear-fractional problem directly. Most importantly, we show that the yield-optimal solution set (like the rate-optimal solution set) is determined by (yield-optimal) elementary flux vectors of the underlying metabolic model. However, yield- and rate-optimal solutions may differ from each other, and hence optimal (biomass or product) yields are not necessarily obtained at solutions with optimal (growth or synthesis) rates. Moreover, we discuss phase planes/production envelopes and yield spaces, in particular, we prove that yield spaces are convex and provide algorithms for their computation. We illustrate our findings by a small example and demonstrate their relevance for metabolic engineering with realistic models of E. coli.

Conclusions

We develop a comprehensive mathematical framework for yield optimization in metabolic models. Our theory is particularly useful for the study and rational modification of cell factories designed under given yield and/or rate requirements.

SUBMITTER: Klamt S

PROVIDER: S-EPMC5992331 | biostudies-literature | 2018 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

A mathematical framework for yield (vs. rate) optimization in constraint-based modeling and applications in metabolic engineering.

Klamt Steffen S Müller Stefan S Regensburger Georg G Zanghellini Jürgen J

Metabolic engineering 20180207

<h4>Background</h4>The optimization of metabolic rates (as linear objective functions) represents the methodical core of flux-balance analysis techniques which have become a standard tool for the study of genome-scale metabolic models. Besides (growth and synthesis) rates, metabolic yields are key parameters for the characterization of biochemical transformation processes, especially in the context of biotechnological applications. However, yields are ratios of rates, and hence the optimization ...[more]

PMID: 29427605

Similar Datasets

Project description:The complement system is part of the innate immune system and plays an important role in the host defense against infectious pathogens. One of the main effects is the opsonization of foreign invaders and subsequent uptake by phagocytosis. Due to the continuous default basal level of active complement molecules, a tight regulation is required to protect the body's own cells (self cells) from opsonization and from complement damage. A major complement regulator is Factor H, which is recruited from the fluid phase and attaches to cell surfaces where it effectively controls complement activation. Besides self cells, pathogens also have the ability to bind Factor H; they can thus escape opsonization and phagocytosis causing severe infections. In order to advance our understanding of the opsonization process at a quantitative level, we developed a mathematical model for the dynamics of the complement system-termed DynaCoSys model-that is based on ordinary differential equations for cell surface-bound molecules and on partial differential equations for concentration profiles of the fluid phase molecules in the environment of cells. This hybrid differential equation approach allows to model the complement cascade focusing on the role of active C3b in the fluid phase and on the cell surface as well as on its inactivation on the cell surface. The DynaCoSys model enables us to quantitatively predict the conditions under which Factor H mediated complement evasion occurs. Furthermore, investigating the quantitative impact of model parameters by a sensitivity analysis, we identify the driving processes of complement activation and regulation in both the self and non-self regime. The two regimes are defined by a critical Factor H concentration on the cell surface and we use the model to investigate the differential impact of complement model parameters on this threshold value. The dynamic modeling on the surface of pathogens are further relevant to understand pathophysiological situations where Factor H mutants and defective Factor H binding to target surfaces results in pathophysiology such as renal and retinal disease. In the future, this DynaCoSys model will be extended to also enable evaluating treatment strategies of complement-related diseases.

Project description:BACKGROUND: An important step in strain optimization is to identify reactions whose activities should be modified to achieve the desired cellular objective. Preferably, these reactions are identified systematically, as the number of possible combinations of reaction modifications could be very large. Over the last several years, a number of computational methods have been described for identifying combinations of reaction modifications. However, none of these methods explicitly address uncertainties in implementing the reaction activity modifications. In this work, we model the uncertainties as probability distributions in the flux carrying capacities of reactions. Based on this model, we develop an optimization method that identifies reactions for flux capacity modifications to predict outcomes with high statistical likelihood. RESULTS: We compare three optimization methods that select an intervention set comprising up- or down-regulation of reaction flux capacity: CCOpt (Chance constrained optimization), DetOpt (Deterministic optimization), and MCOpt (Monte Carlo-based optimization). We evaluate the methods using a Monte Carlo simulation-based method, MCEval (Monte Carlo Evaluations). We present two case studies analyzing a CHO cell and an adipocyte model. The flux capacity distributions required for our methods were estimated from maximal reaction velocities or elementary mode analysis. The intervention set selected by CCOpt consistently outperforms the intervention set selected by DetOpt in terms of tolerance to flux capacity variations. MCEval shows that the optimal flux predicted based on the CCOpt intervention set is more likely to be obtained, in a probabilistic sense, than the flux predicted by DetOpt. The intervention sets identified by CCOpt and MCOpt were similar; however, the exhaustive sampling required by MCOpt incurred significantly greater computational cost. CONCLUSIONS: Maximizing tolerance to variable engineering outcomes (in modifying enzyme activities) can identify intervention sets that statistically improve the desired cellular objective.

Project description:Improved prediction of optimal N fertilizer rates for corn (Zea mays L.) can reduce N losses and increase profits. We tested the ability of the Agricultural Production Systems sIMulator (APSIM) to simulate corn and soybean (Glycine max L.) yields, the economic optimum N rate (EONR) using a 16-year field-experiment dataset from central Iowa, USA that included two crop sequences (continuous corn and soybean-corn) and five N fertilizer rates (0, 67, 134, 201, and 268 kg N ha-1) applied to corn. Our objectives were to: (a) quantify model prediction accuracy before and after calibration, and report calibration steps; (b) compare crop model-based techniques in estimating optimal N rate for corn; and (c) utilize the calibrated model to explain factors causing year to year variability in yield and optimal N. Results indicated that the model simulated well long-term crop yields response to N (relative root mean square error, RRMSE of 19.6% before and 12.3% after calibration), which provided strong evidence that important soil and crop processes were accounted for in the model. The prediction of EONR was more complex and had greater uncertainty than the prediction of crop yield (RRMSE of 44.5% before and 36.6% after calibration). For long-term site mean EONR predictions, both calibrated and uncalibrated versions can be used as the 16-year mean differences in EONR's were within the historical N rate error range (40-50 kg N ha-1). However, for accurate year-by-year simulation of EONR the calibrated version should be used. Model analysis revealed that higher EONR values in years with above normal spring precipitation were caused by an exponential increase in N loss (denitrification and leaching) with precipitation. We concluded that long-term experimental data were valuable in testing and refining APSIM predictions. The model can be used as a tool to assist N management guidelines in the US Midwest and we identified five avenues on how the model can add value toward agronomic, economic, and environmental sustainability.

Dataset Information

A mathematical framework for yield (vs. rate) optimization in constraint-based modeling and applications in metabolic engineering.

Background

Results

Conclusions

Publications

A mathematical framework for yield (vs. rate) optimization in constraint-based modeling and applications in metabolic engineering.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets