Effects of periodic input on the quasi-steady state assumptions for enzyme-catalysed reactions

Balanced truncation for model reduction of biological oscillators

Model reduction is a central problem in mathematical biology. Reduced order models enable modeling of a biological system at different levels of complexity and the quantitative analysis of its properties, like sensitivity to parameter variations and resilience to exogenous perturbations. However, available model reduction methods often fail to capture a diverse range of nonlinear behaviors observed in biology, such as multistability and limit cycle oscillations. The paper addresses this need using differential analysis. This approach leads to a nonlinear enhancement of classical balanced truncation for biological systems whose behavior is not restricted to the stability of a single equilibrium. Numerical results suggest that the proposed framework may be relevant to the approximation of classical models of biological systems.

Bogdanov-Takens bifurcations in the enzyme-catalyzed reaction comprising a branched network

There have been some results on bifurcations of codimension one (such as saddle-node, transcritical, pitchfork) and degenerate Hopf bifurcations for an enzyme-catalyzed reaction system comprising a branched network but no further discussion for bifurcations at its cusp. In this paper we give conditions for the existence of a cusp and compute the parameter curves for the Bogdanov-Takens bifurcation, which induces the appearance of homoclinic orbits and periodic orbits, indicating the tendency to steady-states or a rise of periodic oscillations for the concentrations of the substrate and the product.

A method for zooming of nonlinear models of biochemical systems

BACKGROUND

Models of biochemical systems are typically complex, which may complicate the discovery of cardinal biochemical principles. It is therefore important to single out the parts of a model that are essential for the function of the system, so that the remaining non-essential parts can be eliminated. However, each component of a mechanistic model has a clear biochemical interpretation, and it is desirable to conserve as much of this interpretability as possible in the reduction process. Furthermore, it is of great advantage if we can translate predictions from the reduced model to the original model.

RESULTS

In this paper we present a novel method for model reduction that generates reduced models with a clear biochemical interpretation. Unlike conventional methods for model reduction our method enables the mapping of predictions by the reduced model to the corresponding detailed predictions by the original model. The method is based on proper lumping of state variables interacting on short time scales and on the computation of fraction parameters, which serve as the link between the reduced model and the original model. We illustrate the advantages of the proposed method by applying it to two biochemical models. The first model is of modest size and is commonly occurring as a part of larger models. The second model describes glucose transport across the cell membrane in baker's yeast. Both models can be significantly reduced with the proposed method, at the same time as the interpretability is conserved.

CONCLUSIONS

We introduce a novel method for reduction of biochemical models that is compatible with the concept of zooming. Zooming allows the modeler to work on different levels of model granularity, and enables a direct interpretation of how modifications to the model on one level affect the model on other levels in the hierarchy. The method extends the applicability of the method that was previously developed for zooming of linear biochemical models to nonlinear models.

How molecular should your molecular model be? On the level of molecular detail required to simulate biological networks in systems and synthetic biology

The recent advance of genetic studies and the rapid accumulation of molecular data, together with the increasing performance of computers, led researchers to design more and more detailed mathematical models of biological systems. Many modeling approaches rely on ordinary differential equations (ODE) which are based on standard enzyme kinetics. Michaelis-Menten and Hill functions are indeed commonly used in dynamical models in systems and synthetic biology because they provide the necessary nonlinearity to make the dynamics nontrivial (i.e., limit-cycle oscillations or multistability). For most of the systems modeled, the actual molecular mechanism is unknown, and the enzyme equations should be regarded as phenomenological. In this chapter, we discuss the validity and accuracy of these approximations. In particular, we focus on the validity of the Michaelis-Menten function for open systems and on the use of Hill kinetics to describe transcription rates of regulated genes. Our discussion is illustrated by numerical simulations of prototype systems, including the Repressilator (a genetic oscillator) and the Toggle Switch model (a bistable system). We systematically compare the results obtained with the compact version (based on Michaelis-Menten and Hill functions) with its corresponding developed versions (based on "elementary" reaction steps and mass action laws). We also discuss the use of compact approaches to perform stochastic simulations (Gillespie algorithm). On the basis of these results, we argue that using compact models is suitable to model qualitatively biological systems.

Use and abuse of the quasi-steady-state approximation

The transient kinetic behaviour of an open single enzyme, single substrate reaction is examined. The reaction follows the Van Slyke-Cullen mechanism, a spacial case of the Michaelis-Menten reaction. The analysis is performed both with and without applying the quasi-steady-state approximation. The analysis of the full system shows conditions for biochemical pathway coupling, which yield sustained oscillatory behaviour in the enzyme reaction. The reduced model does not demonstrate this behaviour. The results have important implications in the analysis of open biochemical reactions and the modelling of metabolic systems.

Kinetic constraints for formation of steady states in biochemical networks

The constraint-based analysis has emerged as a useful tool for analysis of biochemical networks. This work introduces the concept of kinetic constraints. It is shown that maximal reaction rates are appropriate constraints only for isolated enzymatic reactions. For biochemical networks, it is revealed that constraints for formation of a steady state require specific relationships between maximal reaction rates of all enzymes. The constraints for a branched network are significantly different from those for a cyclic network. Moreover, the constraints do not require Michaelis-Menten constants for most enzymes, and they only require the constants for the enzymes at the branching or cyclic point. Reversibility of reactions at system boundary or branching point may significantly impact on kinetic constraints. When enzymes are regulated, regulations may impose severe kinetic constraints for the formation of steady states. As the complexity of a network increases, kinetic constraints become more severe. In addition, it is demonstrated that kinetic constraints for networks with co-regulation can be analyzed using the approach. In general, co-regulation enhances the constraints and therefore larger fluctuations in fluxes can be accommodated in the networks with co-regulation. As a first example of the application, we derive the kinetic constraints for an actual network that describes sucrose accumulation in the sugar cane culm, and confirm their validity using numerical simulations.