Abstract:BackgroundSystems 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.MethodsIn 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.ResultsThe 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.ConclusionsThrough 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.

SIMILARITY REDUCTIONS AND INTEGRABILITY FOR THE BRUSSELATOR REACTION–DIFFUSION MODEL WITH SYMBOLIC COMPUTATION

Coexistence of Positive Steady-state Solutions for a Class of Reaction-diffusion Model with Holling-Tanner Term

Model Reducted-order of Linear Discrete Time Systems Based on Bilinear Transformation

Modified Hermite Wavelet Discrete Matrix Approach for the Brusselator Chemical Model

A combined model reduction algorithm for controlled biochemical systems

Analytic Investigation of a Reaction-diffusion Brusselator Model with the Time-space Fractional Derivative

Turing instability and attractor bifurcation for the general Brusselator model

Intrinsic Fluctuations of Chemical Reactions with Different Approaches

Stability and Bifurcation Analysis in a Diffusive Brusselator System with Delayed Feedback Control.

Spatial patterns of the Brusselator model with the asymmetric Lévy diffusion

A novel conformable fractional approach to the Brusselator system with numerical simulation

Amplitude equation for a diffusion–reaction system in presence of complexing reaction with the activator species: the Brusselator model

Multiple Bifurcations and Spatiotemporal Patterns for A Coupled Two-Cell Brusselator Model

The Shigesada-Kawasaki-Teramoto model: conditional symmetries, exact solutions and their properties

A finite element approach to capture Turing patterns of autocatalytic Brusselator model

Bifurcation of a three molecular saturated reaction with impulsive input

A New Mechanism Revealed by Cross-Diffusion-Driven Instability and Double-Hopf Bifurcation in the Brusselator System

Graph-based, dynamics-preserving reduction of (bio)chemical systems

A HYBRID APPROACH FOR SYNCHRONIZING BETWEEN TWO REACTION-DIFFUSION SYSTEMS OF INTEGER- AND FRACTIONAL-ORDER APPLIED ON CERTAIN CHEMICAL MODELS

Kinetic Pattern Formation with Intermolecular Interactions: A Modified Brusselator Model

Constrained Broyden Minimization Combined with the Dimer Method for Locating Transition State of Complex Reactions