Optimization in integrated biochemical systems

As of yet, steady-state optimization in biochemical systems has been limited to a few studies in which networks of fluxes were optimized. These networks of fluxes are represented by linear (stoichiometric) equations that are used as constraints in a linear program, and a flux or a sum of weighted fluxes is used as the objective function. In contrast to networks of fluxes, systems of metabolic processes have not been optimized in a comparable manner. The primary reason is that these types of integrated biochemical systems are full of synergisms, antagonisms, and regulatory mechanisms that can only be captured appropriately with nonlinear models. These models are mathematically complex and difficult to analyze. In most cases it is not even possible to compute, let alone optimize, steady-state solutions analytically. Rare exceptions are S-system representations. These are nonlinear and able to represent virtually all types of dynamic behaviors, but their steady states are characterized by linear equations that can be evaluated both analytically and numerically. The steady-state equations are expressed in terms of the logarithms of the original variables, and because a function and its logarithms of the original variables, and because a function and its logarithm assume their maxima for the same argument, yields or fluxes can be optimized with linear programs expressed in terms of the logarithms of the original variables.