Linlog Kinetics Research Articles

Construction of kinetic models for metabolic reaction networks: Lessons learned in analysing short-term stimulus response data

Construction of dynamic models of large-scale metabolic networks is one of the central issues in the engineering of living cells. However, construction of such models is often hampered by a number of challenges, for example, data availability, compartmentalization and parameter identification coupled with design of in vivo perturbations. As a solution to the latter, short-term perturbation experiments are proposed and are proven to be a useful experimental method to obtain insights into the in vivo kinetic properties of the metabolic pathways. The aim of this work is to construct a kinetic model using the available experimental data obtained by short-term perturbation experiments, where the steady state of a glucose-limited anaerobic chemostat culture of Saccharomyces cerevisiae was perturbed. In constructing the model, we first determined the steady-state flux distribution using the data before the glucose pulse and the known stoichiometry. For the rate expressions, we used approximative linlog kinetics, which allows the enzyme–metabolite kinetic interactions to be represented by an elasticity matrix. We performed a priori model reduction based on timescale analysis and parameter identifiability analysis allowing the information content of the experimental data to be assessed. The final values of the elasticities are estimated by fitting the model to the available short-term kinetic response data. The final model consists of 16 metabolites and 14 reactions. With 25 parameters, the model adequately describes the short-term response of the cells to the glucose perturbation, pointing to the fact that the assumed kinetic interactions in the model are sufficient to account for the observed response.

Hybrid dynamic modeling of Escherichia coli central metabolic network combining Michaelis–Menten and approximate kinetic equations

The construction of dynamic metabolic models at reaction network level requires the use of mechanistic enzymatic rate equations that comprise a large number of parameters. The lack of knowledge on these equations and the difficulty in the experimental identification of their associated parameters, represent nowadays the limiting factor in the construction of such models. In this study, we compare four alternative modeling approaches based on Michaelis–Menten kinetics for the bi-molecular reactions and different types of simplified rate equations for the remaining reactions (generalized mass action, convenience kinetics, lin-log and power-law). Using the mechanistic model for Escherichia coli central carbon metabolism as a benchmark, we investigate the alternative modeling approaches through comparative simulations analyses. The good dynamic behavior and the powerful predictive capabilities obtained using the hybrid model composed of Michaelis–Menten and the approximate lin-log kinetics indicate that this is a possible suitable approach to model complex large-scale networks where the exact rate laws are unknown.

Towards a genome-scale kinetic model of cellular metabolism

BackgroundAdvances in bioinformatic techniques and analyses have led to the availability of genome-scale metabolic reconstructions. The size and complexity of such networks often means that their potential behaviour can only be analysed with constraint-based methods. Whilst requiring minimal experimental data, such methods are unable to give insight into cellular substrate concentrations. Instead, the long-term goal of systems biology is to use kinetic modelling to characterize fully the mechanics of each enzymatic reaction, and to combine such knowledge to predict system behaviour.ResultsWe describe a method for building a parameterized genome-scale kinetic model of a metabolic network. Simplified linlog kinetics are used and the parameters are extracted from a kinetic model repository. We demonstrate our methodology by applying it to yeast metabolism. The resultant model has 956 metabolic reactions involving 820 metabolites, and, whilst approximative, has considerably broader remit than any existing models of its type. Control analysis is used to identify key steps within the system.ConclusionsOur modelling framework may be considered a stepping-stone toward the long-term goal of a fully-parameterized model of yeast metabolism. The model is available in SBML format from the BioModels database (BioModels ID: MODEL1001200000) and at http://www.mcisb.org/resources/genomescale/.

Evaluating the integration of proteomic data for the prediction of intracellular fluxes after knockout experiments

So far, few large scale kinetic models of metabolic networks have been successfully constructed. The main reasons for this are not only the associated mathematical complexity, but also the large number of unknown kinetic parameters required in the rate equations to define the system. In contrast to kinetic models, the constraint-based modelling approach bypasses these difficulties by using basically only stoichiometric information with certain physicochemical constraints to delimit the solution space without large fitted parameter sets. Although these constraint-based models are highly relevant to predict feasible steady-state fluxes under a diverse range of genetic and environmental conditions, the steady-state assumption may oversimplify cellular behaviour and cannot predict time-course profiles. To overcome these problems, combining these two approaches appears as a reasonable alternative to modelling large-scale metabolic networks. On the other hand, several of the experimental data required for model construction are often rare and in this way it is usually assumed that the enzyme concentrations are constant.In this work, we used a central carbon metabolic network of E. coli to investigate whether including high throughput enzyme concentration data into a kinetic model allows improved predictions of metabolic flux distributions in response to single knockouts perturbations. For this purpose, an E. coli model, based on results obtained from flux balance analysis (FBA) and approximate lin-log kinetics was constructed. The intracellular fluxes distributions, obtained using this model, were compared with published in vivo measurements.

Model reduction and a priori kinetic parameter identifiability analysis using metabolome time series for metabolic reaction networks with linlog kinetics

In this work, we present a time-scale analysis based model reduction and parameter identifiability analysis method for metabolic reaction networks. The method uses the information obtained from short term chemostat perturbation experiments. We approximate the time constant of each metabolite pool by their turn-over time and classify the pools accordingly into two groups: fast and slow pools. We performed a priori model reduction, neglecting the dynamic term of the fast pools. By making use of the linlog approximative kinetics, we obtained a general explicit solution for the fast pools in terms of the slow pools by elaborating the degenerate algebraic system resulting from model reduction. The obtained relations yielded also analytical relations between a subset of kinetic parameters. These relations also allow to realize an analytical model reduction using lumped reaction kinetics. After solving these theoretical identifiability problems and performing model reduction, we carried out a Monte Carlo approach to study the practical identifiability problems. We illustrated the methodology on model reduction and theoretical/practical identifiability analysis on an example system representing the glycolysis in Saccharomyces cerevisiae cells.

Something from nothing − bridging the gap between constraint‐based and kinetic modelling

Two divergent modelling methodologies have been adopted to increase our understanding of metabolism and its regulation. Constraint-based modelling highlights the optimal path through a stoichiometric network within certain physicochemical constraints. Such an approach requires minimal biological data to make quantitative inferences about network behaviour; however, constraint-based modelling is unable to give an insight into cellular substrate concentrations. In contrast, kinetic modelling aims to characterize fully the mechanics of each enzymatic reaction. This approach suffers because parameterizing mechanistic models is both costly and time-consuming. In this paper, we outline a method for developing a kinetic model for a metabolic network, based solely on the knowledge of reaction stoichiometries. Fluxes through the system, estimated by flux balance analysis, are allowed to vary dynamically according to linlog kinetics. Elasticities are estimated from stoichiometric considerations. When compared to a popular branched model of yeast glycolysis, we observe an excellent agreement between the real and approximate models, despite the absence of (and indeed the requirement for) experimental data for kinetic constants. Moreover, using this particular methodology affords us analytical forms for steady state determination, stability analyses and studies of dynamical behaviour.

Determination of elasticities, concentration and flux control coefficients from transient metabolite data using linlog kinetics

This paper presents a practical approach to estimate the kinetic parameters of a metabolic network from in vivo kinetics experiments. This method is based on the linlog kinetics format (Visser and Heijnen, 2003, Metab. Eng. 5(3), 164–176; Wu et al., 2004, Eur. J. Biochem. 271, 3348–3359), of which the kinetic parameters, called elasticities, are estimated by an iterative linear optimization followed by non-linear optimization, from transient metabolite concentration data which are directly obtainable from rapid pulse experiments. In this way, not only the parameters are estimated but also a full kinetic model, based on linlog kinetics, is developed. The obtained elasticities also allow immediate calculation of all control coefficients. As an in silico case study, the estimation of elasticities of a linear pathway is presented. The method is shown to be able to estimate the elasticities quite accurately and to be robust toward errors in the metabolite data originating from sampling and measurement inaccuracy. The method allows experimental redesign to get more accurate estimated parameters and accommodates various types of experimentally applied disturbances in the pathway: changes in independent metabolites, dependent metabolites or enzyme levels/activities.

Optimal re-design of primary metabolism in Escherichia coli using linlog kinetics

This paper examines the validity of the linlog approach, which was recently developed in our laboratory, by comparison of two different kinetic models for the metabolic network of Escherichia coli. The first model is a complete mechanistic model; the second is an approximative model in which linlog kinetics are applied. The parameters of the linlog model (elasticities) are derived from the mechanistic model.Three different optimization cases are examined. In all cases, the objective is to calculate the enzyme levels that maximize a certain flux while keeping the total amount of enzyme constant and preventing large changes of metabolite concentrations. For an average variation of metabolite levels of 10% and individual changes of a factor 2, the predicted enzyme levels, metabolite concentrations and fluxes of both models are highly similar. This similarity holds for changes in enzyme level of a factor 4–6 and for changes in fluxes up to a factor 6.In all three cases, the predicted optimal enzyme levels could neither have been found by intuition-based approaches, nor on basis of flux control coefficients. This demonstrates that kinetic models are essential tools in Metabolic Engineering. In this respect, the linlog approach is a valuable extension of MCA, since it allows construction of kinetic models, based on MCA parameters, that can be used for constrained optimization problems and are valid for large changes of metabolite and enzyme levels.

Dynamic simulation and metabolic re-design of a branched pathway using linlog kinetics

This paper presents a new mathematical framework for modeling of in vivo dynamics and for metabolic re-design: the linlog approach. This approach is an extension of metabolic control analysis (MCA), valid for large changes of enzyme and metabolite levels. Furthermore, the presented framework combines MCA with kinetic modeling, thereby also combining the merits of both approaches. The linlog framework includes general expressions giving the steady-state fluxes and metabolite concentrations as a function of enzyme levels and extracellular concentrations, and a metabolic design equation that allows direct calculation of required enzyme levels for a desired steady state when control and response coefficients are available. Expressions giving control coefficients as a function of the enzyme levels are also derived. The validity of the linlog approximation in metabolic modeling is demonstrated by application of linlog kinetics to a branched pathway with moiety conservation, reversible reactions and allosteric interactions. Results show that the linlog approximation is able to describe the non-linear dynamics of this pathway very well for concentration changes up to a factor 20. Also the metabolic design equation was tested successfully.