Constraint-based analysis concentration constraints

Despite its widely recognized limitations, flux balance analysis has resulted in a large number of successful applications [35, 67, 72 74], including several extensions and refinements. See Ref. [247] for a recent review. Of particular interest are recent efforts to augment the stoichiometric balance equations with thermodynamic constraints providing a link between concentration and flux in the constraint-based analysis of metabolic networks [74, 149, 150]. For a more comprehensive review, we refer to the very readable monograph of Palsson [50]. 

Analysis of biochemical systems, with their behaviors constrained by the known system stoichiometry, falls under the broad heading constraint-based analysis, a methodology that allows us to explore computationally metabolic fluxes and concentrations constrained by the physical chemical laws of mass conservation and thermodynamics. This chapter introduces the mathematical formulation of the constraints on reaction fluxes and reactant concentrations that arise from the stoichiometry of an integrated network and are the basis of constraint-based analysis. 

Introducing the chemical potential (or free energy) and the thermodynamic constraint provides a solid physical chemistry foundation for the constraint-based analysis approach to metabolic systems analysis. Treatment of the network thermodynamics not only improves the accuracy of the predictions on the steady state fluxes, but can also be used to make predictions on the steady state concentrations of metabolites. To see this, we substitute the relation between reaction Gibbs free energy (ArG ) of the th reaction and the concentrations of biochemical reactants... 

Mathematical optimization deals with determining values for a set of unknown variables x, X2, , x , which best satisfy (optimize) some mathematical objective quantified by a scalar function of the unknown variables, F(xi, X2, , xn). The function F is termed the objective function bounds on the variables, along with mathematical dependencies between them, are termed constraints. Constraint-based analysis of metabolic systems requires definition of the constraints acting on biochemical variables (fluxes, concentrations, enzyme activities) and determining appropriate objective functions useful in determining the behavior of metabolic systems. 

Applying mass-balance and thermodynamic constraints typically leaves one without a precisely defined (unique) solution for reaction fluxes and reactant concentration, but instead with a mathematically constrained feasible space for these variables. Exploration of this feasible space is the purview of constraint-based analysis. It has so far been left unstated that any application in this area starts with the determination of the reactions in a system, from which the stoichiometric matrix arises. This first step, network reconstruction, integrates genomic and proteomic data to determine carefully the enzymes present in an organism, cell, or subcellular compartment. The network reconstruction process is described elsewhere [107]. 

As we shall see, linear algebraic constraints arising from steady state mass balance form the basis of metabolic flux analysis (MFA) and flux balance analysis (FBA). Thermodynamic laws, while introducing inherent non-linearities into the mathematical description of the feasible flux space, allow determination of feasible reaction directions and facilitate the introduction of reactant concentrations to the constraint-based framework. 

The foregoing analysis was mainly based on time taken as the primary constraints and concentration taken as the secondary constraints. It has been shown that this primary constraints can be bypassed by merely using storage facilities, as long as the secondary constraints is met. However, it is also possible to take concentration as the primary constraints and time as the secondary constraints. The procedure and the associated outcome of such an analysis form the subject of the following section. 

One can view biochemical systems as represented at the most basic level as networks of given stoichiometry. Whether the steady state or the kinetic behavior is explored, the stoichiometry constrains the feasible behavior according to mass balance and the laws of thermodynamics. As we have seen in this chapter, some analysis is possible based solely on the stoichiometric structure of a given system. Mass balance provides linear constraints on reaction fluxes non-linear thermodynamic constraints provide information about feasible flux directions and reactant concentrations. 

Geochemical constraints on acid deposition come from the foram record of a spike in Sr/ Sr across the K/T boundaiy and from base cation leaching at terrestrial K/T sites in eastern Montana. A statistical analysis of the foram data of Martin and Macdougall [1] yields an enhancement in Sr/ Sr of 25 x 10 9x10 We considered two sources of this enhancement Sr leaching from impact ejecta and increased continental weathering due to acid deposition. Impact ejecta with the Sr concentration and ratio of Chiexulub melt rock, 300... 

This paper reports cryoscopic measurements for NiS04 aqueous solutions (0.00 > t /°C > - 0.21), and demonstrates the dependence of the calculated value of K on the assumptions in the complexation model. Based on the arguments of Brown and Prue, it is clear that use of the SIT as is done in the TDB project imposes constraints on the data analysis, and that these constraints should lead to a specific value for K. However, recalculations show that the association constant is also strongly dependent on the concentration at which the data set is truncated. For the purposes of the present review, recalculations were done using the SIT, but only results for molalities < 0.03 were used in the selection of the value of log . When data for higher concentrations were included, the calculated value for the association constant increased, and log,Q = 300 was obtained if the entire data set was used. 

Since a blend containing high concentration of cavitated rubber particles becomes cellular solid (porous) rather than continuous material, Eq. 11.15 does not apply to it any longer and any analysis of the plastic zone size must be based on yield criteria appropriate for the porous solid. Free from the constraints of continuum mechanics, the cavitated plastic zones formed in polymer blends are able to increase substantially in radius even under plane-strain conditions (Bucknall and Paul 2009). 

In the last twenty years, various non-deterministic methods have been developed to deal with optimum design under environmental uncertainties. These methods can be classified into two main branches, namely reliability-based methods and robust-based methods. The reliability methods, based on the known probabiUty distribution of the random parameters, estimate the probability distribution of the system s response, and are predominantly used for risk analysis by computing the probability of system failure. However, variation is not minimized in reliability approaches (Siddall, 1984) because they concentrate on rare events at the tail of the probability distribution (Doltsinis and Kang, 2004). The robust design methods are commonly based on multiobjective minimization problems. The are commonly indicated as Multiple Objective Robust Optimization (MORO) and find a set of optimal solutions that optimise a performance index in terms of mean value and, at the same time, minimize its resulting dispersion due to input parameters uncertainty. The final solution is less sensitive to the parameters variation but eventually maintains feasibility with regards probabilistic constraints. This is achieved by the optimization of the design vector in order to make the performance minimally sensitive to the various causes of variation. 

---