The tricarboxylic acid cycle in Dictyostelium discoideum. III. Analysis of steady state and dynamic behavior

A Step-by-Step Guide to Using BioNetFit

BioNetFit is a software tool designed for solving parameter identification problems that arise in the development of rule-based models. It solves these problems through curve fitting (i.e., nonlinear regression). BioNetFit is compatible with deterministic and stochastic simulators that accept BioNetGen language (BNGL)-formatted files as inputs, such as those available within the BioNetGen framework. BioNetFit can be used on a laptop or stand-alone multicore workstation as well as on many Linux clusters, such as those that use the Slurm Workload Manager to schedule jobs. BioNetFit implements a metaheuristic population-based global optimization procedure, an evolutionary algorithm (EA), to minimize a user-defined objective function, such as a residual sum of squares (RSS) function. BioNetFit also implements a bootstrapping procedure for determining confidence intervals for parameter estimates. Here, we provide step-by-step instructions for using BioNetFit to estimate the values of parameters of a BNGL-encoded model and to define bootstrap confidence intervals. The process entails the use of several plain-text files, which are processed by BioNetFit and BioNetGen. In general, these files include (1) one or more EXP files, which each contains (experimental) data to be used in parameter identification/bootstrapping; (2) a BNGL file containing a model section, which defines a (rule-based) model, and an actions section, which defines simulation protocols that generate GDAT and/or SCAN files with model predictions corresponding to the data in the EXP file(s); and (3) a CONF file that configures the fitting/bootstrapping job and that defines algorithmic parameter settings.

Evaluation of an S-system root-finding method for estimating parameters in a metabolic reaction model

In a mathematical model, estimation of parameters from time-series data of metabolic concentrations in cells is a challenging task. However, it seems that a promising approach for such estimation has not yet been established. Biochemical Systems Theory (BST) is a powerful methodology to construct a power-law type model for a given metabolic reaction system and to then characterize it efficiently. In this paper, we discuss the use of an S-system root-finding method (S-system method) to estimate parameters from time-series data of metabolite concentrations. We demonstrate that the S-system method is superior to the Newton-Raphson method in terms of the convergence region and iteration number. We also investigate the usefulness of a translocation technique and a complex-step differentiation method toward the practical application of the S-system method. The results indicate that the S-system method is useful to construct mathematical models for a variety of metabolic reaction networks.

Mathematical Modeling and Dynamic Simulation of Metabolic Reaction Systems Using Metabolome Time Series Data

The high-throughput acquisition of metabolome data is greatly anticipated for the complete understanding of cellular metabolism in living organisms. A variety of analytical technologies have been developed to acquire large-scale metabolic profiles under different biological or environmental conditions. Time series data are useful for predicting the most likely metabolic pathways because they provide important information regarding the accumulation of metabolites, which implies causal relationships in the metabolic reaction network. Considerable effort has been undertaken to utilize these data for constructing a mathematical model merging system properties and quantitatively characterizing a whole metabolic system in toto. However, there are technical difficulties between benchmarking the provision and utilization of data. Although, hundreds of metabolites can be measured, which provide information on the metabolic reaction system, simultaneous measurement of thousands of metabolites is still challenging. In addition, it is nontrivial to logically predict the dynamic behaviors of unmeasurable metabolite concentrations without sufficient information on the metabolic reaction network. Yet, consolidating the advantages of advancements in both metabolomics and mathematical modeling remain to be accomplished. This review outlines the conceptual basis of and recent advances in technologies in both the research fields. It also highlights the potential for constructing a large-scale mathematical model by estimating model parameters from time series metabolome data in order to comprehensively understand metabolism at the systems level.

An Efficient and Very Accurate Method for Calculating Steady-State Sensitivities in Metabolic Reaction Systems

Stability and sensitivity analyses of biological systems require the ad hocwriting of computer code, which is highly dependent on the particular model and burdensome for large systems. We propose a very accurate strategy to overcome this challenge. Its core concept is the conversion of the model into the format of biochemical systems theory (BST), which greatly facilitates the computation of sensitivities. First, the steady state of interest is determined by integrating the model equations toward the steady state and then using a Newton-Raphson method to fine-tune the result. The second step of conversion into the BST format requires several instances of numerical differentiation. The accuracy of this task is ensured by the use of a complex-variable Taylor scheme for all differentiation steps. The proposed strategy is implemented in a new software program, COSMOS, which automates the stability and sensitivity analysis of essentially arbitrary ODE models in a quick, yet highly accurate manner. The methods underlying the process are theoretically analyzed and illustrated with four representative examples: a simple metabolic reaction model; a model of aspartate-derived amino acid biosynthesis; a TCA-cycle model; and a modified TCA-cycle model. COSMOS has been deposited to https://github.com/BioprocessdesignLab/COSMOS.

PENDISC: a simple method for constructing a mathematical model from time-series data of metabolite concentrations

The availability of large-scale datasets has led to more effort being made to understand characteristics of metabolic reaction networks. However, because the large-scale data are semi-quantitative, and may contain biological variations and/or analytical errors, it remains a challenge to construct a mathematical model with precise parameters using only these data. The present work proposes a simple method, referred to as PENDISC (arameter stimation in a on-mensionalized -system with onstraints), to assist the complex process of parameter estimation in the construction of a mathematical model for a given metabolic reaction system. The PENDISC method was evaluated using two simple mathematical models: a linear metabolic pathway model with inhibition and a branched metabolic pathway model with inhibition and activation. The results indicate that a smaller number of data points and rate constant parameters enhances the agreement between calculated values and time-series data of metabolite concentrations, and leads to faster convergence when the same initial estimates are used for the fitting. This method is also shown to be applicable to noisy time-series data and to unmeasurable metabolite concentrations in a network, and to have a potential to handle metabolome data of a relatively large-scale metabolic reaction system. Furthermore, it was applied to aspartate-derived amino acid biosynthesis in Arabidopsis thaliana plant. The result provides confirmation that the mathematical model constructed satisfactorily agrees with the time-series datasets of seven metabolite concentrations.

Biochemical Systems Theory: A Review

Biochemical systems theory (BST) is the foundation for a set of analytical andmodeling tools that facilitate the analysis of dynamic biological systems. This paper depicts major developments in BST up to the current state of the art in 2012. It discusses its rationale, describes the typical strategies and methods of designing, diagnosing, analyzing, and utilizing BST models, and reviews areas of application. The paper is intended as a guide for investigators entering the fascinating field of biological systems analysis and as a resource for practitioners and experts.

Steady-state global optimization of metabolic non-linear dynamic models through recasting into power-law canonical models

Background

Design of newly engineered microbial strains for biotechnological purposes would greatly benefit from the development of realistic mathematical models for the processes to be optimized. Such models can then be analyzed and, with the development and application of appropriate optimization techniques, one could identify the modifications that need to be made to the organism in order to achieve the desired biotechnological goal. As appropriate models to perform such an analysis are necessarily non-linear and typically non-convex, finding their global optimum is a challenging task. Canonical modeling techniques, such as Generalized Mass Action (GMA) models based on the power-law formalism, offer a possible solution to this problem because they have a mathematical structure that enables the development of specific algorithms for global optimization.

Results

Based on the GMA canonical representation, we have developed in previous works a highly efficient optimization algorithm and a set of related strategies for understanding the evolution of adaptive responses in cellular metabolism. Here, we explore the possibility of recasting kinetic non-linear models into an equivalent GMA model, so that global optimization on the recast GMA model can be performed. With this technique, optimization is greatly facilitated and the results are transposable to the original non-linear problem. This procedure is straightforward for a particular class of non-linear models known as Saturable and Cooperative (SC) models that extend the power-law formalism to deal with saturation and cooperativity.

Conclusions

Our results show that recasting non-linear kinetic models into GMA models is indeed an appropriate strategy that helps overcoming some of the numerical difficulties that arise during the global optimization task.

Highly accurate computation of dynamic sensitivities in metabolic reaction systems by a Taylor series method

We have previously developed the software for calculation of dynamic sensitivities, SoftCADS, in which one can calculate dynamic sensitivities with high accuracy by just setting the differential equations for metabolite concentrations. However, SoftCADS did not always provide calculated values with the machine accuracy of a computer, although a Taylor series method was employed to numerically solve the differential equations. This is because numerical derivatives calculated from an approximate formula were directly used in the derivation of the differential equations for sensitivities from those for metabolite concentrations. The present work therefore attempts to further enhance the performance of SoftCADS, including not only the accuracies of the calculated values but also the calculation time. To overcome the problem, the approximate formula is expanded into a Taylor series in time and the first-term value of the series is replaced by the exact coefficient on the second term of the flux function expanded into a Taylor series in an independent or dependent variable. The result reveals that this replacement certainly provides not only numerical derivatives but also dynamic sensitivities with superhigh accuracies comparable to the machine accuracy, regardless of the degree of stiffness of the differential equations. Moreover, a comparison indicates that the improved SoftCADS shortens the calculation time of the dynamic sensitivities without reducing their accuracies, even when the simplest approximate derivative formula is used.

A reliable Taylor series-based computational method for the calculation of dynamic sensitivities in large-scale metabolic reaction systems: algorithm and software evaluation

Dynamic sensitivity analysis has become an important tool to successfully characterize all sorts of biological systems. However, when the analysis is carried out on large scale systems, it becomes imperative to employ a highly accurate computational method in order to obtain reliable values. Furthermore, the preliminary laborious mathematical operations required by current software before the computation of dynamic sensitivities makes it inconvenient for a significant number of unacquainted users. To satisfy these needs, the present work investigates a newly developed algorithm consisting of a combination of Taylor series method that can directly execute Taylor expansions for simultaneous non-linear-differential equations and a simple but highly-accurate numerical differentiation method based on finite-difference formulas. Applications to three examples of biochemical systems indicate that the proposed method makes it possible to compute the dynamic sensitivity values with highly-reliable accuracies and also allows to readily compute them by setting up only the differential equations for metabolite concentrations in the computer program. Also, it is found that the Padé approximation introduced in the Taylor series method shortens the computation time greatly because it stabilizes the computation so that it allows us to use larger stepsizes in the numerical integration. Consequently, the calculated results suggest that the proposed computational method, in addition to being user-friendly, makes it possible to perform dynamic sensitivity analysis in large-scale metabolic reaction systems both efficiently and reliably.

A gradual update method for simulating the steady-state solution of stiff differential equations in metabolic circuits

Numerical simulation of differential equation systems plays a major role in the understanding of how metabolic network models generate particular cellular functions. On the other hand, the classical and technical problems for stiff differential equations still remain to be solved, while many elegant algorithms have been presented. To relax the stiffness problem, we propose new practical methods: the gradual update of differential-algebraic equations based on gradual application of the steady-state approximation to stiff differential equations, and the gradual update of the initial values in differential-algebraic equations. These empirical methods show a high efficiency for simulating the steady-state solutions for the stiff differential equations that existing solvers alone cannot solve. They are effective in extending the applicability of dynamic simulation to biochemical network models.

Coordination of the dynamics of yeast sphingolipid metabolism during the diauxic shift

BACKGROUND

The diauxic shift in yeast requires cells to coordinate a complicated response that involves numerous genes and metabolic processes. It is unknown whether responses of this type are mediated in vivo through changes in a few "key" genes and enzymes, which are mathematically characterized by high sensitivities, or whether they are based on many small changes in genes and enzymes that are not particularly sensitive. In contrast to global assessments of changes in gene or protein interaction networks, we study here control aspects of the diauxic shift by performing a detailed analysis of one specific pathway-sphingolipid metabolism-which is known to have signaling functions and is associated with a wide variety of stress responses.

RESULTS

The approach uses two components: publicly available sets of expression data of sphingolipid genes and a recently developed Generalized Mass Action (GMA) mathematical model of the sphingolipid pathway. In one line of exploration, we analyze the sensitivity of the model with respect to enzyme activities, and thus gene expression. Complementary to this approach, we convert the gene expression data into changes in enzyme activities and then predict metabolic consequences by means of the mathematical model. It was found that most of the sensitivities in the model are low in magnitude, but that some stand out as relatively high. This information was then deployed to test whether the cell uses a few of the very sensitive pathway steps to mount a response or whether the control is distributed throughout the pathway. Pilot experiments confirm qualitatively and in part quantitatively the predictions of a group of metabolite simulations.

CONCLUSION

The results indicate that yeast coordinates sphingolipid mediated changes during the diauxic shift through an array of small changes in many genes and enzymes, rather than relying on a strategy involving a few select genes with high sensitivity. This study also highlights a novel approach in coupling data mining with mathematical modeling in order to evaluate specific metabolic pathways.

Optimization of biotechnological systems through geometric programming

Background

In the past, tasks of model based yield optimization in metabolic engineering were either approached with stoichiometric models or with structured nonlinear models such as S-systems or linear-logarithmic representations. These models stand out among most others, because they allow the optimization task to be converted into a linear program, for which efficient solution methods are widely available. For pathway models not in one of these formats, an Indirect Optimization Method (IOM) was developed where the original model is sequentially represented as an S-system model, optimized in this format with linear programming methods, reinterpreted in the initial model form, and further optimized as necessary.

Results

A new method is proposed for this task. We show here that the model format of a Generalized Mass Action (GMA) system may be optimized very efficiently with techniques of geometric programming. We briefly review the basics of GMA systems and of geometric programming, demonstrate how the latter may be applied to the former, and illustrate the combined method with a didactic problem and two examples based on models of real systems. The first is a relatively small yet representative model of the anaerobic fermentation pathway in S. cerevisiae, while the second describes the dynamics of the tryptophan operon in E. coli. Both models have previously been used for benchmarking purposes, thus facilitating comparisons with the proposed new method. In these comparisons, the geometric programming method was found to be equal or better than the earlier methods in terms of successful identification of optima and efficiency.

Conclusion

GMA systems are of importance, because they contain stoichiometric, mass action and S-systems as special cases, along with many other models. Furthermore, it was previously shown that algebraic equivalence transformations of variables are sufficient to convert virtually any types of dynamical models into the GMA form. Thus, efficient methods for optimizing GMA systems have multifold appeal.

A simple and highly accurate numerical differentiation method for sensitivity analysis of large-scale metabolic reaction systems

Numerical differentiation is known to be one of the most difficult numerical calculation methods to obtain reliable calculated values at all times. A simple numerical differentiation method using a combination of finite-difference formulas, derived by approximation of Taylor-series equations, is investigated in order to efficiently perform the sensitivity analysis of large-scale metabolic reaction systems. A result of the application to four basic mathematical functions reveals that the use of the eight-point differentiation formula with a non-dimensionalized stepsize close to 0.01 mostly provides more than 14 digits of accuracy in double precision for the numerical derivatives. Moreover, a result of the application to the modified TCA cycle model indicates that the numerical differentiation method gives the calculated values of steady-state metabolite concentrations within a range of round-off error and also makes it possible to transform the Michaelis-Menten equations into the S-system equations having the kinetic orders whose accuracies are mostly more than 14 significant digits. Because of the simple structure of the numerical differentiation formula and its promising high accuracy, it is evident that the present numerical differentiation method is useful for the analysis of large-scale metabolic reaction systems according to the systematic procedure of BST.

Mechanistic model of fuel selection in the muscle

Fuel selection in human muscle is key to explaining insulin resistance. In obesity and type 2 diabetes mellitus, there is an increased content of lipid within and around muscle fibers. Changes in muscle fuel partitioning of lipid, between oxidation and storage of fat, contribute to the accumulation of intramuscular triglycerides and to the pathogenesis of both obesity and type 2 diabetes mellitus. A mathematical model of the aggregated metabolism in skeletal muscle was developed and the effects of fuel selection for lean and obese individuals under fasting conditions, insulin-stimulated conditions, and oscillating insulin conditions were examined. Model results were complementary to prior observations that elevated lipid oxidation during insulin-stimulated conditions is correlated with insulin resistance. The model also adequately simulated metabolic inflexibility between fat and glucose oxidation in the obese individual. A novel sensitivity analysis indicated the strong interaction effects of parameters of glucose and lipid oxidation pathways on the variables of each pathway.

Mathematical model predicts a critical role for osteoclast autocrine regulation in the control of bone remodeling

Bone remodeling occurs asynchronously at multiple sites in the adult skeleton and involves resorption by osteoclasts, followed by formation of new bone by osteoblasts. Disruptions in bone remodeling contribute to the pathogenesis of disorders such as osteoporosis, osteoarthritis, and Paget's disease. Interactions among cells of osteoblast and osteoclast lineages are critical in the regulation of bone remodeling. We constructed a mathematical model of autocrine and paracrine interactions among osteoblasts and osteoclasts that allowed us to calculate cell population dynamics and changes in bone mass at a discrete site of bone remodeling. The model predicted different modes of dynamic behavior: a single remodeling cycle in response to an external stimulus, a series of internally regulated cycles of bone remodeling, or unstable behavior similar to pathological bone remodeling in Paget's disease. Parametric analysis demonstrated that the mode of dynamic behavior in the system depends strongly on the regulation of osteoclasts by autocrine factors, such as transforming growth factor beta. Moreover, simulations demonstrated that nonlinear dynamics of the system may explain the differing effects of immunosuppressants on bone remodeling in vitro and in vivo. In conclusion, the mathematical model revealed that interactions among osteoblasts and osteoclasts result in complex, nonlinear system behavior, which cannot be deduced from studies of each cell type alone. The model will be useful in future studies assessing the impact of cytokines, growth factors, and potential therapies on the overall process of remodeling in normal bone and in pathological conditions such as osteoporosis and Paget's disease.

Issues of exploration: human health and wellbeing during a mission to Mars

Today, the tools are in our hands to enable us to travel away from our home planet and become citizens of the solar system. Even now, we are seriously beginning to develop the robust infrastructure that will make the 21st century the Century of Space Travel. But this bold step must be taken with due concern for the health, safety and wellbeing of future space explorers. Our long experience with space biomedical research convinces us that, if we are to deal effectively with the medical and biomedical issues of exploration, then dramatic and bold steps are also necessary in this field. We can no longer treat the human body as if it were composed of muscles, bones, heart and brain acting independently. Instead, we must lead the effort to develop a fully integrated view of the body, with all parts connected and fully interacting in a realistic way. This paper will present the status of current (2000) plans by the National Space Biomedical Research Institute to initiate research in this area of integrative physiology and medicine. Specifically, three example projects are discussed as potential stepping stones towards the ultimate goal of producing a digital human. These projects relate to developing a functional model of the human musculoskeletal system and the heart.

Modeling, optimization and experimental assessment of continuous L-(-)-carnitine production by Escherichia coli cultures

In a previous paper Cánovas et al. (Biotechnol Bioeng 2002;77:764-775) presented a model for L-(-)-carnitine production using Escherichia coli O44 K74, in a cell-recycle bioreactor for the biotransformation of crotonobetaine into L-carnitine. In this work we optimize this biotechnological setup and experimentally verify the predicted optimal parameter profiles. Provided with a reliable and robust S-system description of the cell-bioreactor combined system, we applied the Indirect Optimization Method described by Torres et al. (Biotechnol Bioeng 1997;55(5):758-772; Food Technol Biotechnol 1998;36(3):177-184). This optimization approach provides different parameter value profiles, all of which are compatible with the cell physiology and the bioreactor operating conditions, that yield increased rates of L-(-)-carnitine production. Three parameters were seen to be of critical importance for maximizing L-(-)-carnitine production: the dilution rate, the initial crotonobetaine concentration, and the carnitine dehydratase activity. When the first two were changed in the experimental setup, there was a 74% increase in the L-(-)-carnitine production rate, performance that was in close agreement with the predictions of the model. In accordance with the optimized solution, a further improvement (90% increase in the L-(-)-carnitine production rate) could be attained by over-expressing up to 5 times the carnitine dehydratase basal activity. Thus the optimization approach shown herein provides experimental evidence of a new strategy which demonstrates the possible variables that can be subjected to modifications compatible with the cell physiology and bioreactor operating conditions, and which are able to yield increased rates of L-(-)-carnitine production.

Modeling and simulation: tools for metabolic engineering

Mathematical modeling is one of the key methodologies of metabolic engineering. Based on a given metabolic model different computational tools for the simulation, data evaluation, systems analysis, prediction, design and optimization of metabolic systems have been developed. The currently used metabolic modeling approaches can be subdivided into structural models, stoichiometric models, carbon flux models, stationary and nonstationary mechanistic models and models with gene regulation. However, the power of a model strongly depends on its basic modeling assumptions, the simplifications made and the data sources used. Model validation turns out to be particularly difficult for metabolic systems. The different modeling approaches are critically reviewed with respect to their potential and benefits for the metabolic engineering cycle. Several tools that have emerged from the different modeling approaches including structural pathway synthesis, stoichiometric pathway analysis, metabolic flux analysis, metabolic control analysis, optimization of regulatory architectures and the evaluation of rapid sampling experiments are discussed.

Cardiac energy metabolism: models of cellular respiration

The heart requires a large amount of energy to sustain both ionic homeostasis and contraction. Under normal conditions, adenosine triphosphate (ATP) production meets this demand. Hence, there is a complex regulatory system that adjusts energy production to meet this demand. However, the mechanisms for this control are a topic of active debate. Energy metabolism can be divided into three main stages: substrate delivery to the tricarboxylic acid (TCA) cycle, the TCA cycle, and oxidative phosphorylation. Each of these processes has multiple control points and exerts control over the other stages. This review discusses the basic stages of energy metabolism, mechanisms of control, and the mathematical and computational models that have been used to study these mechanisms.

Power-law modeling based on least-squares criteria: consequences for system analysis and simulation

The power-law formalism was initially derived as a Taylor series approximation in logarithmic space for kinetic rate-laws. The resulting models, either as generalized mass action (GMA) or as S-systems models, allow to characterize the target system and to simulate its dynamical behavior in response to external perturbations and parameter changes. This approach has been succesfully used as a modeling tool in many applications from cell metabolism to population dynamics. Without leaving the general formalism, we recently proposed to derive the power-law representation in an alternative way that uses least-squares (LS) minimization instead of the traditional derivation based on Taylor series [B. Hernández-Bermejo, V. Fairén, A. Sorribas, Math. Biosci. 161 (1999) 83-94]. It was shown that the resulting LS power-law mimics the target rate-law in a wider range of concentration values than the classical power-law, and that the prediction of the steady-state using the LS power-law is closer to the actual steady-state of the target system. However, many implications of this alternative approach remained to be established. We explore some of them in the present work. Firstly, we extend the definition of the LS power-law within a given operating interval in such a way that no preferred operating point is selected. Besides providing an alternative to the classical Taylor power-law, that can be considered a particular case when the operating interval is reduced to a single point, the LS power-law so defined is consistent with the results that can be obtained by fitting experimental data points. Secondly, we show that the LS approach leads to a system description, either as an S-system or a GMA model, in which the systemic properties (such as the steady-state prediction or the log-gains) appear averaged over the corresponding interval when compared with the properties that can be computed from Taylor-derived models in different operating points within the considered operating range. Finally, we also show that the LS description leads to a global, accurate description of the system when it is submitted to external forcing.

Power-law modeling based on least-squares minimization criteria

The power-law formalism has been successfully used as a modeling tool in many applications. The resulting models, either as Generalized Mass Action or as S-systems models, allow one to characterize the target system and to simulate its dynamical behavior in response to external perturbations and parameter changes. The power-law formalism was first derived as a Taylor series approximation in logarithmic space for kinetic rate-laws. The especial characteristics of this approximation produce an extremely useful systemic representation that allows a complete system characterization. Furthermore, their parameters have a precise interpretation as local sensitivities of each of the individual processes and as rate-constants. This facilitates a qualitative discussion and a quantitative estimation of their possible values in relation to the kinetic properties. Following this interpretation, parameter estimation is also possible by relating the systemic behavior to the underlying processes. Without leaving the general formalism, in this paper we suggest deriving the power-law representation in an alternative way that uses least-squares minimization. The resulting power-law mimics the target rate-law in a wider range of concentration values than the classical power-law. Although the implications of this alternative approach remain to be established, our results show that the predicted steady-state using the least-squares power-law is closest to the actual steady-state of the target system.

Development of fractal kinetic theory for enzyme-catalysed reactions and implications for the design of biochemical pathways

Recent evidence has shown that elementary bimolecular reactions under dimensionally-restricted conditions, such as those that might occur within cells when reactions are confined to two-dimensional membranes and one-dimensional channels, do not follow traditional mass-action kinetics, but fractal kinetics. The power-law formalism, which provides the context for examining the kinetics under these conditions, is used here to examine the implications of fractal kinetics in a simple pathway of reversible reactions. Starting with elementary chemical kinetics, we proceed to characterise the equilibrium behaviour of a simple bimolecular reaction, derive a generalised set of conditions for microscopic reversibility, and develop the fractal kinetic rate law for a reversible Michaelis-Menten mechanism. Having established this fractal kinetic framework, we go on to analyse the steady-state behaviour and temporal response of a pathway characterised by both the fundamental and quasi-steady-state equations. These results are contrasted with those for the fundamental and quasi-steady-state equations based on traditional mass-action kinetics. Finally, we compare the accuracy of three local representations based on both fractal and mass-action kinetics. The results with fractal kinetics show that the equilibrium ratio is a function of the amount of material in a closed system, and that the principle of microscopic reversibility has a more general manifestation that imposes new constraints on the set of fractal kinetic orders. Fractal kinetics in a biochemical pathway allow an increase in flux to occur with less accumulation of pathway intermediates and a faster temporal response than is the case with traditional kinetics. These conclusions are obtained regardless of the level of representation considered. Thus, fractal kinetics provide a novel means to achieve important features of pathway design.

Redesigning metabolic networks using mathematical programming

The design of new generation bioprocessing plants is increasingly dependent on the design of process-compatible microorganisms. The latter, whether through genetic or physiological manipulations, can be greatly assisted by metabolic engineering. An emerging powerful tool in metabolic engineering research is computer-assisted cell design using mathematical programming. In this work, the problem of optimizing cellular metabolic networks has been formulated as a Mixed Integer Nonlinear Programming (MINLP) model. The model can assist genetic engineers to identify which cellular enzymes should be modified, and the new levels of activity required to produce an optimal network. Results are presented from the tricarboxylic acid cycle in Dictyostelium discoideum. Copyright 1998 John Wiley & Sons, Inc.

Application of biochemical systems theory to metabolism in human red blood cells. Signal propagation and accuracy of representation

Human erythrocytes are among the simplest of cells. Many of their enzymes have been characterized kinetically using steady-state methods in vitro, and several investigators have assembled this kinetic information into mathematical models of the integrated system. However, despite its relative simplicity, the integrated behavior of erythrocyte metabolism is still complex and not well understood. Errors will inevitably be encountered in any such model because of this complexity; thus, the construction of an integrative model must be considered an iterative process of assessment and refinement. In a previous study, we selected a recent model of erythrocyte metabolism as our starting point and took it through three stages of model assessment and refinement using systematic strategies provided by biochemical systems theory. At each stage deficiencies were diagnosed, putative remedies were identified, and modifications consistent with existing experimental evidence were incorporated into the working model. In this paper we address two issues: the propagation of biochemical signals within the metabolic network, and the accuracy of kinetic representation. The analysis of signal propagation reveals the importance of glutathione peroxidase, transaldolase, and the concentration of total glutathione in determining systemic behavior. It also reveals a highly amplified diversion of flux between the pathways of pentose phosphate and nucleotide metabolism. In determining the range of accurate representation based on alternative kinetic formalisms we discovered large discrepancies. These were identified with the behavior of the model represented within the Michaelis-Menten formalism. This model fails to exhibit a nominal steady state when the activity of glutathione peroxidase is decreased by as little as 9%. Our current understanding, as embodied in this working model, is in need of further refinement, and the results presented in this paper suggest areas of the model where such effort might profitably be concentrated.

Comparative characterization of the fermentation pathway of Saccharomyces cerevisiae using biochemical systems theory and metabolic control analysis: model definition and nomenclature

Mathematical tools that involve the determination of systemic responses to small changes in metabolites or enzymes have demonstrated their utility for analyzing metabolic pathways. The different methodologies based on these ideas allow for modeling and analyzing biochemical pathways focusing on the coordinate behavior of the whole system. However, one must become familiar with the difference in nomenclature and methodology to relate the models and results obtained by applying these techniques and to appreciate their potential for answering fundamental questions about biochemical systems. In the following three papers we show how this can be facilitated by comparing the nomenclature, methodology, and results of the two leading techniques in this area, metabolic control analysis and biochemical systems theory, using a model of the fermentation pathway in Saccharomyces cerevisiae as a reference system. In the present paper we review the nomenclature, technical concepts, and related experimental measurements while creating a practical dictionary for the reference system that makes the relatedness of the two approaches more apparent. In the second paper, subtitled Steady-State Analysis, we show that both approaches give the same picture for many systemic responses of the reference system. In the third paper of this series, subtitled Model Validation and Dynamic Behavior, we show that the quality of the model can be assessed by studying the sensitivity to changes in the system parameters. We hope to illustrate the usefulness of these tools in providing an interpretation of the experimental measurements in a specific metabolic pathway.

Comparative characterization of the fermentation pathway of Saccharomyces cerevisiae using biochemical systems theory and metabolic control analysis: model validation and dynamic behavior

In the first two papers of this series (immediately preceding, this issue), we characterized the steady-state properties of a model of a fermentation pathway in Saccharomyces cerevisiae in four experimental conditions. In each of these conditions, the pictures obtained by metabolic control analysis and biochemical systems theory were coincident, which illustrates the relatedness of the two approaches. In this paper we analyze the quality of this description by means of the tools available within biochemical systems theory, and we show that in some of the experimental conditions studied the system is poorly characterized. The most critical condition corresponds to the immobilization of the cells at pH 5.5, in which the kinetic characterization appears to be inaccurate. Furthermore, sensitivity analysis and the study of the local steady-state stability identify the most critical parameters. The results of these analyses are confirmed by the predictions of the dynamic response of the model using its S-system representation. This illustrates the utility of these tools and warns against using the steady-state characterization without testing its validity.