Dataset Information

A generalised enzyme kinetic model for predicting the behaviour of complex biochemical systems.

ABSTRACT: Quasi steady-state enzyme kinetic models are increasingly used in systems modelling. The Michaelis Menten model is popular due to its reduced parameter dimensionality, but its low-enzyme and irreversibility assumption may not always be valid in the in vivo context. Whilst the total quasi-steady state assumption (tQSSA) model eliminates the reactant stationary assumptions, its mathematical complexity is increased. Here, we propose the differential quasi-steady state approximation (dQSSA) kinetic model, which expresses the differential equations as a linear algebraic equation. It eliminates the reactant stationary assumptions of the Michaelis Menten model without increasing model dimensionality. The dQSSA was found to be easily adaptable for reversible enzyme kinetic systems with complex topologies and to predict behaviour consistent with mass action kinetics in silico. Additionally, the dQSSA was able to predict coenzyme inhibition in the reversible lactate dehydrogenase enzyme, which the Michaelis Menten model failed to do. Whilst the dQSSA does not account for the physical and thermodynamic interactions of all intermediate enzyme-substrate complex states, it is proposed to be suitable for modelling complex enzyme mediated biochemical systems. This is due to its simpler application, reduced parameter dimensionality and improved accuracy.

SUBMITTER: Wong MK

PROVIDER: S-EPMC4383669 | biostudies-literature | 2015

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

A generalised enzyme kinetic model for predicting the behaviour of complex biochemical systems.

Wong Martin Kin Lok MK Krycer James Robert JR Burchfield James Geoffrey JG James David Ernest DE Kuncic Zdenka Z

FEBS open bio 20150309

Quasi steady-state enzyme kinetic models are increasingly used in systems modelling. The Michaelis Menten model is popular due to its reduced parameter dimensionality, but its low-enzyme and irreversibility assumption may not always be valid in the in vivo context. Whilst the total quasi-steady state assumption (tQSSA) model eliminates the reactant stationary assumptions, its mathematical complexity is increased. Here, we propose the differential quasi-steady state approximation (dQSSA) kinetic ...[more]

PMID: 25859426

Similar Datasets

Project description:Dealing with a system of first-order reactions is a recurrent issue in chemometrics, especially in the analysis of data obtained by spectroscopic methods applied on complex biological systems. We argue that global multiexponential fitting, the still common way to solve such problems, has serious weaknesses compared to contemporary methods of sparse modeling. Combining the advantages of group lasso and elastic net-the statistical methods proven to be very powerful in other areas-we created an optimization problem tunable from very sparse to very dense distribution over a large pre-defined grid of time constants, fitting both simulated and experimental multiwavelength spectroscopic data with high computational efficiency. We found that the optimal values of the tuning hyperparameters can be selected by a machine-learning algorithm based on a Bayesian optimization procedure, utilizing widely used or novel versions of cross-validation. The derived algorithm accurately recovered the true sparse kinetic parameters of an extremely complex simulated model of the bacteriorhodopsin photocycle, as well as the wide peak of hypothetical distributed kinetics in the presence of different noise levels. It also performed well in the analysis of the ultrafast experimental fluorescence kinetics data detected on the coenzyme FAD in a very wide logarithmic time window. We conclude that the primary application of the presented algorithms-implemented in available software-covers a wide area of studies on light-induced physical, chemical, and biological processes carried out with different spectroscopic methods. The demand for this kind of analysis is expected to soar due to the emerging ultrafast multidimensional infrared and electronic spectroscopic techniques that provide very large and complex datasets. In addition, simulations based on our methods could help in designing the technical parameters of future experiments for the verification of particular hypothetical models.

Project description:A persistent obstacle for constructing kinetic models of metabolism is uncertainty in the kinetic properties of enzymes. Currently, available methods for building kinetic models can cope indirectly with uncertainties by integrating data from different biological levels and origins into models. In this study, we use the recently proposed computational approach iSCHRUNK (in Silico Approach to Characterization and Reduction of Uncertainty in the Kinetic Models), which combines Monte Carlo parameter sampling methods and machine learning techniques, in the context of Bayesian inference. Monte Carlo parameter sampling methods allow us to exploit synergies between different data sources and generate a population of kinetic models that are consistent with the available data and physicochemical laws. The machine learning allows us to data-mine the a priori generated kinetic parameters together with the integrated datasets and derive posterior distributions of kinetic parameters consistent with the observed physiology. In this work, we used iSCHRUNK to address a design question: can we identify which are the kinetic parameters and what are their values that give rise to a desired metabolic behavior? Such information is important for a wide variety of studies ranging from biotechnology to medicine. To illustrate the proposed methodology, we performed Metabolic Control Analysis, computed the flux control coefficients of the xylose uptake (XTR), and identified parameters that ensure a rate improvement of XTR in a glucose-xylose co-utilizing S. cerevisiae strain. Our results indicate that only three kinetic parameters need to be accurately characterized to describe the studied physiology, and ultimately to design and control the desired responses of the metabolism. This framework paves the way for a new generation of methods that will systematically integrate the wealth of available omics data and efficiently extract the information necessary for metabolic engineering and synthetic biology decisions.

Project description:BACKGROUND:Systems Biology continues to produce increasingly large models of complex biochemical reaction networks. In applications requiring, for example, parameter estimation, the use of agent-based modelling approaches, or real-time simulation, this growing model complexity can present a significant hurdle. Often, however, not all portions of a model are of equal interest in a given setting. In such situations methods of model reduction offer one possible approach for addressing the issue of complexity by seeking to eliminate those portions of a pathway that can be shown to have the least effect upon the properties of interest. METHODS:In this paper a model reduction algorithm bringing together the complementary aspects of proper lumping and empirical balanced truncation is presented. Additional contributions include the development of a criterion for the selection of state-variable elimination via conservation analysis and use of an 'averaged' lumping inverse. This combined algorithm is highly automatable and of particular applicability in the context of 'controlled' biochemical networks. RESULTS:The algorithm is demonstrated here via application to two examples; an 11 dimensional model of bacterial chemotaxis in Escherichia coli and a 99 dimensional model of extracellular regulatory kinase activation (ERK) mediated via the epidermal growth factor (EGF) and nerve growth factor (NGF) receptor pathways. In the case of the chemotaxis model the algorithm was able to reduce the model to 2 state-variables producing a maximal relative error between the dynamics of the original and reduced models of only 2.8% whilst yielding a 26 fold speed up in simulation time. For the ERK activation model the algorithm was able to reduce the system to 7 state-variables, incurring a maximal relative error of 4.8%, and producing an approximately 10 fold speed up in the rate of simulation. Indices of controllability and observability are additionally developed and demonstrated throughout the paper. These provide insight into the relative importance of individual reactants in mediating a biochemical system's input-output response even for highly complex networks. CONCLUSIONS:Through application, this paper demonstrates that combined model reduction methods can produce a significant simplification of complex Systems Biology models whilst retaining a high degree of predictive accuracy. In particular, it is shown that by combining the methods of proper lumping and empirical balanced truncation it is often possible to produce more accurate reductions than can be obtained by the use of either method in isolation.

Dataset Information

A generalised enzyme kinetic model for predicting the behaviour of complex biochemical systems.

Publications

A generalised enzyme kinetic model for predicting the behaviour of complex biochemical systems.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets