Investigation of kinetic-order sensitivities in metabolic reaction networks.

Kinetic-order sensitivity (the ratio of relative change in a dependent variable to the relative change in a kinetic order in a power-law-type differential equation) has recently become an important indicator in metabolic pathway analysis using mathematical models with parameter values determined from time-series data on cellular metabolite concentrations. Here, we discuss a potential problem in calculating kinetic-order sensitivities. When the steady-state metabolite concentration is less than unity, a slight increase in the kinetic order changes the metabolite concentration in the incorrect direction, yielding a kinetic-order sensitivity value with an incorrect sign. This is caused by a property of the power-law function (y=Xn): when X is less than unity, y decreases for a larger positive n or for a smaller absolute value of negative n. We propose two solutions. The first is to directly calculate the kinetic-order sensitivities and then reverse the sign of the relevant value if a steady-state metabolite concentration less than unity is involved. The second involves calculation of the kinetic-order sensitivities after setting all metabolite concentrations to values greater than unity (e.g., by changing the units from mM to µM). The latter method changes the absolute values of the kinetic-order sensitivities according to the magnitude of a multiplication factor, because kinetic-order sensitivities do not have unique values. Nevertheless, since the normalized absolute values exhibit an almost identical distribution, it should not be difficult to identify which kinetic order has greater effect, although kinetic order rankings may change slightly under different calculation conditions.

A new parametric method to smooth time-series data of metabolites in metabolic networks.

Mathematical modeling of large-scale metabolic networks usually requires smoothing of metabolite time-series data to account for measurement or biological errors. Accordingly, the accuracy of smoothing curves strongly affects the subsequent estimation of model parameters. Here, an efficient parametric method is proposed for smoothing metabolite time-series data, and its performance is evaluated. To simplify parameter estimation, the method uses S-system-type equations with simple power law-type efflux terms. Iterative calculation using this method was found to readily converge, because parameters are estimated stepwise. Importantly, smoothing curves are determined so that metabolite concentrations satisfy mass balances. Furthermore, the slopes of smoothing curves are useful in estimating parameters, because they are probably close to their true behaviors regardless of errors that may be present in the actual data. Finally, calculations for each differential equation were found to converge in much less than one second if initial parameters are set at appropriate (guessed) values.

PASMet: a web-based platform for prediction, modelling and analyses of metabolic systems.

PASMet (Prediction, Analysis and Simulation of Metabolic networks) is a web-based platform for proposing and verifying mathematical models to understand the dynamics of metabolism. The advantages of PASMet include user-friendliness and accessibility, which enable biologists and biochemists to easily perform mathematical modelling. PASMet offers a series of user-functions to handle the time-series data of metabolite concentrations. The functions are organised into four steps: (i) Prediction of a probable metabolic pathway and its regulation; (ii) Construction of mathematical models; (iii) Simulation of metabolic behaviours; and (iv) Analysis of metabolic system characteristics. Each function contains various statistical and mathematical methods that can be used independently. Users who may not have enough knowledge of computing or programming can easily and quickly analyse their local data without software downloads, updates or installations. Users only need to upload their files in comma-separated values (CSV) format or enter their model equations directly into the website. Once the time-series data or mathematical equations are uploaded, PASMet automatically performs computation on server-side. Then, users can interactively view their results and directly download them to their local computers. PASMet is freely available with no login requirement at http://pasmet.riken.jp/ from major web browsers on Windows, Mac and Linux operating systems.

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.

Using dynamic sensitivities to characterize metabolic reaction systems.

Metabolite concentrations in cells are governed by enzyme kinetics in the metabolic reaction system. One can analyze how these concentrations depend on system variables such as enzyme activities by computing system sensitivities, which generally vary with time. Dynamic sensitivities, i.e., time-varying sensitivities, reflect the time-dependent response of the metabolic reaction network to perturbations. Unfortunately, dynamic sensitivity profiles are not commonly used in the analysis of metabolic reaction systems. In the present study, we demonstrate the use of dynamic logarithmic gains, i.e., normalized time-varying sensitivities, to gain insights into the dynamic behavior of metabolic networks. A biosynthetic reaction model of aromatic amino acids proposed by other researchers is used as a case study. The model system is analyzed using the dynamic logarithmic gains in parallel with simulations of the time-transient behavior of metabolite concentrations and metabolic fluxes. The result indicates that the influences of independent variables are most pronounced just after perturbations and the effects of perturbations on metabolite concentration at early times can be larger than those at steady state. These findings suggest that it is important to perform dynamic sensitivity analysis in addition to steady-state analysis. Furthermore, the rankings of the bottleneck ranking indicators, defined as the product of dynamic logarithmic gain and metabolite concentration, for three desired amino acids reveal that the degree of bottleneck for each enzyme changes with time. In conclusion, the dynamic logarithmic gains are not only useful for analyzing metabolic reaction systems but also can offer additional insights on the transient behavior of the system over steady state sensitivities, leading to a proper design of metabolic systems.

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 (Parameter Estimation in a N on- DImensionalized S-system with Constraints), 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.

Estimation of kinetic parameters in an S-system equation model for a metabolic reaction system using the Newton-Raphson method.

Metabolic reaction systems can be modeled easily in terms of S-system type equations if their metabolic maps are available. This study therefore proposes a method for estimating parameters in decoupled S-system equations on the basis of the Newton-Raphson method and elucidates the performance of this estimation method. Parameter estimation from the time-course data of metabolite concentrations reveals that the parameters estimated are highly accurate, indicating that the estimation algorithm has been constructed correctly. The number of iterations is small and the calculation converges in a very short time (usually less than 1s). The method is also applied to time course data with noise and found to estimate parameters efficiently. Results indicate that the present method has the potential to be extended to a method for estimating parameters in large-scale metabolic reaction systems.

A U-system approach for predicting metabolic behaviors and responses based on an alleged metabolic reaction network.

BACKGROUND: Progress in systems biology offers sophisticated approaches toward a comprehensive understanding of biological systems. Yet, computational analyses are held back due to difficulties in determining suitable model parameter values from experimental data which naturally are subject to biological fluctuations. The data may also be corrupted by experimental uncertainties and sometimes do not contain all information regarding variables that cannot be measured for technical reasons. RESULTS: We show here a streamlined approach for the construction of a coarse model that allows us to set up dynamic models with minimal input information. The approach uses a hybrid between a pure mass action system and a generalized mass action (GMA) system in the framework of biochemical systems theory (BST) with rate constants of 1, normal kinetic orders of 1, and -0.5 and 0.5 for inhibitory and activating effects, named Unity (U)-system. The U-system model does not necessarily fit all data well but is often sufficient for predicting metabolic behavior of metabolites which cannot be simultaneously measured, identifying inconsistencies between experimental data and the assumed underlying pathway structure, as well as predicting system responses to a modification of gene or enzyme. The U-system approach was validated with small, generic systems and implemented to model a large-scale metabolic reaction network of a higher plant, Arabidopsis. The dynamic behaviors obtained by predictive simulations agreed with actually available metabolomic time-series data, identified probable errors in the experimental datasets, and estimated probable behavior of unmeasurable metabolites in a qualitative manner. The model could also predict metabolic responses of Arabidopsis with altered network structures due to genetic modification. CONCLUSIONS: The U-system approach can effectively predict metabolic behaviors and responses based on structures of an alleged metabolic reaction network. Thus, it can be a useful first-line tool of data analysis, model diagnostics and aid the design of next-step experiments.

SS-mPMG and SS-GA: tools for finding pathways and dynamic simulation of metabolic networks.

Metabolomics analysis tools can provide quantitative information on the concentration of metabolites in an organism. In this paper, we propose the minimum pathway model generator tool for simulating the dynamics of metabolite concentrations (SS-mPMG) and a tool for parameter estimation by genetic algorithm (SS-GA). SS-mPMG can extract a subsystem of the metabolic network from the genome-scale pathway maps to reduce the complexity of the simulation model and automatically construct a dynamic simulator to evaluate the experimentally observed behavior of metabolites. Using this tool, we show that stochastic simulation can reproduce experimentally observed dynamics of amino acid biosynthesis in Arabidopsis thaliana. In this simulation, SS-mPMG extracts the metabolic network subsystem from published databases. The parameters needed for the simulation are determined using a genetic algorithm to fit the simulation results to the experimental data. We expect that SS-mPMG and SS-GA will help researchers to create relevant metabolic networks and carry out simulations of metabolic reactions derived from metabolomics data.

Identification of a metabolic reaction network from time-series data of metabolite concentrations.

Recent development of high-throughput analytical techniques has made it possible to qualitatively identify a number of metabolites simultaneously. Correlation and multivariate analyses such as principal component analysis have been widely used to analyse those data and evaluate correlations among the metabolic profiles. However, these analyses cannot simultaneously carry out identification of metabolic reaction networks and prediction of dynamic behaviour of metabolites in the networks. The present study, therefore, proposes a new approach consisting of a combination of statistical technique and mathematical modelling approach to identify and predict a probable metabolic reaction network from time-series data of metabolite concentrations and simultaneously construct its mathematical model. Firstly, regression functions are fitted to experimental data by the locally estimated scatter plot smoothing method. Secondly, the fitted result is analysed by the bivariate Granger causality test to determine which metabolites cause the change in other metabolite concentrations and remove less related metabolites. Thirdly, S-system equations are formed by using the remaining metabolites within the framework of biochemical systems theory. Finally, parameters including rate constants and kinetic orders are estimated by the Levenberg-Marquardt algorithm. The estimation is iterated by setting insignificant kinetic orders at zero, i.e., removing insignificant metabolites. Consequently, a reaction network structure is identified and its mathematical model is obtained. Our approach is validated using a generic inhibition and activation model and its practical application is tested using a simplified model of the glycolysis of Lactococcus lactis MG1363, for which actual time-series data of metabolite concentrations are available. The results indicate the usefulness of our approach and suggest a probable pathway for the production of lactate and acetate. The results also indicate that the approach pinpoints a probable strong inhibition of lactate on the glycolysis pathway.

Investigation of the performance of fermentation processes using a mathematical model including effects of metabolic bottleneck and toxic product on cells.

A number of recent research studies have focused on theoretical and experimental investigation of a bottleneck in a metabolic reaction network. However, there is no study on how the bottleneck affects the performance of a fermentation process when a product is highly toxic and remarkably influences the growth and death of cells. The present work therefore studies the effect of bottleneck on product concentrations under different product toxicity conditions. A generalized bottleneck model in a fed-batch fermentation is constructed including both the bottleneck and the product influences on cell growth and death. The simulation result reveals that when the toxic product strongly influences the cell growth and death, the final product concentration is hardly changed even if the bottleneck is removed, whereas it is markedly changed by the degree of product toxicity. The performance of an ethanol fermentation process is also discussed as a case example to validate this result. In conclusion, when the product is highly toxic, one cannot expect a significant increase in the final product concentration even if removing the bottleneck; rather, it may be more effective to somehow protect the cells so that they can continuously produce the product.

Calculation errors of time-varying flux control coefficients obtained from elasticity coefficients by means of summation and connectivity theorems in metabolic control analysis.

This paper investigates the accuracy of a matrix method proposed by other researchers to calculate time-varying flux control coefficients (dynamic FCCs) from elasticity coefficients by means of summation and connectivity theorems in the framework of metabolic control analysis. A mathematical model for the fed-batch penicillin V fermentation process is used as a case example for discussion. Calculated results reveal that this method produces significant calculation errors because the theorems are essentially valid only in steady state, although it may provide rough time-transient behaviors of FCCs. Strictly, therefore, dynamic FCCs should be directly calculated from the differential equations for metabolite concentrations and sensitivities.