Systems approaches to modelling pathways and networks.

Peer reviewedPreprin

Model simplification of signal transduction pathway networks via a hybrid inference strategy

A full-scale mathematical model of cellular networks normally involves a large number of variables and parameters. How to effectively develop manageable and reliable models is crucial for effective computation, analysis and design of such systems. The aim of model simplification is to eliminate parts of a model that are unimportant for the properties of interest. In this work, a model reduction strategy via hybrid inference is proposed for signal pathway networks. It integrates multiple techniques including conservation analysis, local sensitivity analysis, principal component analysis and flux analysis to identify the reactions and variables that can be considered to be eliminated from the full-scale model. Using an I·B-NF-·B signalling pathway model as an example, simulation analysis demonstrates that the simplified model quantitatively predicts the dynamic behaviours of the network

Formulating genome-scale kinetic models in the post-genome era.

The biological community is now awash in high-throughput data sets and is grappling with the challenge of integrating disparate data sets. Such integration has taken the form of statistical analysis of large data sets, or through the bottom-up reconstruction of reaction networks. While progress has been made with statistical and structural methods, large-scale systems have remained refractory to dynamic model building by traditional approaches. The availability of annotated genomes enabled the reconstruction of genome-scale networks, and now the availability of high-throughput metabolomic and fluxomic data along with thermodynamic information opens the possibility to build genome-scale kinetic models. We describe here a framework for building and analyzing such models. The mathematical analysis challenges are reflected in four foundational properties, (i) the decomposition of the Jacobian matrix into chemical, kinetic and thermodynamic information, (ii) the structural similarity between the stoichiometric matrix and the transpose of the gradient matrix, (iii) the duality transformations enabling either fluxes or concentrations to serve as the independent variables and (iv) the timescale hierarchy in biological networks. Recognition and appreciation of these properties highlight notable and challenging new in silico analysis issues

'BioNessie(G) - a grid enabled biochemical networks simulation environment

The simulation of biochemical networks provides insight and understanding about the underlying biochemical processes and pathways used by cells and organisms. BioNessie is a biochemical network simulator which has been developed at the University of Glasgow. This paper describes the simulator and focuses in particular on how it has been extended to benefit from a wide variety of high performance compute resources across the UK through Grid technologies to support larger scale simulations

BioNessie - a grid enabled biochemical networks simulation environment

The simulation of biochemical networks provides insight and understanding about the underlying biochemical processes and pathways used by cells and organisms. BioNessie is a biochemical network simulator which has been developed at the University of Glasgow. This paper describes the simulator and focuses in particular on how it has been extended to benefit from a wide variety of high performance compute resources across the UK through Grid technologies to support larger scale simulations

Inference of nonlinear state-space models for sandwich-type lateral flow immunoassay using extended Kalman filtering

Copyright [2011] IEEE. This material is posted here with permission of the IEEE. Such permission of the IEEE does not in any way imply IEEE endorsement of any of Brunel University's products or services. Internal or personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or redistribution must be obtained from the IEEE by writing to [email protected]. By choosing to view this document, you agree to all provisions of the copyright laws protecting it.In this paper, a mathematical model for sandwichtype lateral flow immunoassay is developed via short available time series. A nonlinear dynamic stochastic model is considered that consists of the biochemical reaction system equations and the observation equation. After specifying the model structure, we apply the extend Kalman filter (EKF) algorithm for identifying both the states and parameters of the nonlinear state-space model. It is shown that the EKF algorithm can accurately identify the parameters and also predict the system states in the nonlinear dynamic stochastic model through an iterative procedure by using a small number of observations. The identified mathematical model provides a powerful tool for testing the system hypotheses and also inspecting the effects from various design parameters in a both rapid and inexpensive way. Furthermore, by means of the established model, the dynamic changes of the concentration of antigens and antibodies can be predicted, thereby making it possible for us to analyze, optimize and design the properties of lateral flow immunoassay devices.This work was supported in part by the International Science and Technology Cooperation Project of China under Grant 2009DFA32050, Natural Science Foundation of Fujian Province of China under Grants 2009J01280 and 2009J01281

Reduction of dynamical biochemical reaction networks in computational biology

Biochemical networks are used in computational biology, to model the static and dynamical details of systems involved in cell signaling, metabolism, and regulation of gene expression. Parametric and structural uncertainty, as well as combinatorial explosion are strong obstacles against analyzing the dynamics of large models of this type. Multi-scaleness is another property of these networks, that can be used to get past some of these obstacles. Networks with many well separated time scales, can be reduced to simpler networks, in a way that depends only on the orders of magnitude and not on the exact values of the kinetic parameters. The main idea used for such robust simplifications of networks is the concept of dominance among model elements, allowing hierarchical organization of these elements according to their effects on the network dynamics. This concept finds a natural formulation in tropical geometry. We revisit, in the light of these new ideas, the main approaches to model reduction of reaction networks, such as quasi-steady state and quasi-equilibrium approximations, and provide practical recipes for model reduction of linear and nonlinear networks. We also discuss the application of model reduction to backward pruning machine learning techniques

Comment on "Regularizing capacity of metabolic networks"

In a recent paper, Marr, Muller-Linow and Hutt [Phys. Rev. E 75, 041917 (2007)] investigate an artificial dynamic system on metabolic networks. They find a less complex time evolution of this dynamic system in real networks, compared to networks of reference models. The authors argue that this suggests that metabolic network structure is a major factor behind the stability of biochemical steady states. We reanalyze the same kind of data using a dynamic system modeling actual reaction kinetics. The conclusions about stability, from our analysis, are inconsistent with those of Marr et al. We argue that this issue calls for a more detailed type of modeling

Elasticity sampling links thermodynamics to metabolic control

Metabolic networks can be turned into kinetic models in a predefined steady state by sampling the reaction elasticities in this state. Elasticities for many reversible rate laws can be computed from the reaction Gibbs free energies, which are determined by the state, and from physically unconstrained saturation values. Starting from a network structure with allosteric regulation and consistent metabolic fluxes and concentrations, one can sample the elasticities, compute the control coefficients, and reconstruct a kinetic model with consistent reversible rate laws. Some of the model variables are manually chosen, fitted to data, or optimised, while the others are computed from them. The resulting model ensemble allows for probabilistic predictions, for instance, about possible dynamic behaviour. By adding more data or tighter constraints, the predictions can be made more precise. Model variants differing in network structure, flux distributions, thermodynamic forces, regulation, or rate laws can be realised by different model ensembles and compared by significance tests. The thermodynamic forces have specific effects on flux control, on the synergisms between enzymes, and on the emergence and propagation of metabolite fluctuations. Large kinetic models could help to simulate global metabolic dynamics and to predict the effects of enzyme inhibition, differential expression, genetic modifications, and their combinations on metabolic fluxes. MATLAB code for elasticity sampling is freely available