Coarse graining of master equations with fast and slow states

We propose a general method for simplifying master equations by eliminating from the description rapidly evolving states. The physical recipe we impose is the suppression of these states and a renormalization of the rates of all the surviving states. In some cases, this decimation procedure can be analytically carried out and is consistent with other analytical approaches, like in the problem of the random walk in a double-well potential. We discuss the application of our method to nontrivial examples: diffusion in a lattice with defects and a model of an enzymatic reaction outside the steady state regime.Comment: 9 pages, 9 figures, final version (new subsection and many minor improvements

Generalized Haldane Equation and Fluctuation Theorem in the Steady State Cycle Kinetics of Single Enzymes

Enyzme kinetics are cyclic. We study a Markov renewal process model of single-enzyme turnover in nonequilibrium steady-state (NESS) with sustained concentrations for substrates and products. We show that the forward and backward cycle times have idential non-exponential distributions: \QQ_+(t)=\QQ_-(t). This equation generalizes the Haldane relation in reversible enzyme kinetics. In terms of the probabilities for the forward ($p_+$) and backward ($p_-$) cycles, $k_BT\ln(p_+/p_-)$ is shown to be the chemical driving force of the NESS, $\Delta\mu$. More interestingly, the moment generating function of the stochastic number of substrate cycle $\nu(t)$, follows the fluctuation theorem in the form of Kurchan-Lebowitz-Spohn-type symmetry. When $\lambda$ = $\Delta\mu/k_BT$, we obtain the Jarzynski-Hatano-Sasa-type equality: $\equiv$ 1 for all $t$, where $\nu\Delta\mu$ is the fluctuating chemical work done for sustaining the NESS. This theory suggests possible methods to experimentally determine the nonequilibrium driving force {\it in situ} from turnover data via single-molecule enzymology.Comment: 4 pages, 3 figure

Rigorous elimination of fast stochastic variables from the linear noise approximation using projection operators

The linear noise approximation (LNA) offers a simple means by which one can study intrinsic noise in monostable biochemical networks. Using simple physical arguments, we have recently introduced the slow-scale LNA (ssLNA) which is a reduced version of the LNA under conditions of timescale separation. In this paper, we present the first rigorous derivation of the ssLNA using the projection operator technique and show that the ssLNA follows uniquely from the standard LNA under the same conditions of timescale separation as those required for the deterministic quasi-steady state approximation. We also show that the large molecule number limit of several common stochastic model reduction techniques under timescale separation conditions constitutes a special case of the ssLNA.Comment: 10 pages, 1 figure, submitted to Physical Review E; see also BMC Systems Biology 6, 39 (2012

Michaelis-Menten Relations for Complex Enzymatic Networks

All biological processes are controlled by complex systems of enzymatic chemical reactions. Although the majority of enzymatic networks have very elaborate structures, there are many experimental observations indicating that some turnover rates still follow a simple Michaelis-Menten relation with a hyperbolic dependence on a substrate concentration. The original Michaelis-Menten mechanism has been derived as a steady-state approximation for a single-pathway enzymatic chain. The validity of this mechanism for many complex enzymatic systems is surprising. To determine general conditions when this relation might be observed in experiments, enzymatic networks consisting of coupled parallel pathways are investigated theoretically. It is found that the Michaelis-Menten equation is satisfied for specific relations between chemical rates, and it also corresponds to the situation with no fluxes between parallel pathways. Our results are illustrated for simple models. The importance of the Michaelis-Menten relationship and derived criteria for single-molecule experimental studies of enzymatic processes are discussed.Comment: 10 pages, 4 figure

Blood coagulation dynamics: mathematical modeling and stability results

The hemostatic system is a highly complex multicomponent biosystem that under normal physiologic conditions maintains the fluidity of blood. Coagulation is initiated in response to endothelial surface vascular injury or certain biochemical stimuli, by the exposure of plasma to Tissue Factor (TF), that activates platelets and the coagulation cascade, inducing clot formation, growth and lysis. In recent years considerable advances have contributed to understand this highly complex process and some mathematical and numerical models have been developed. However, mathematical models that are both rigorous and comprehensive in terms of meaningful experimental data, are not available yet. In this paper a mathematical model of coagulation and fibrinolysis in flowing blood that integrates biochemical, physiologic and rheological factors, is revisited. Three-dimensional numerical simulations are performed in an idealized stenosed blood vessel where clot formation and growth are initialized through appropriate boundary conditions on a prescribed region of the vessel wall. Stability results are obtained for a simplified version of the clot model in quiescent plasma, involving some of the most relevant enzymatic reactions that follow Michaelis-Menten kinetics, and having a continuum of equilibria.CEMAT/IST through FCT [PTDC/MAT/68166/2006]; Czech Science Foundation [201/09/0917]; Grant Agency of the Academy of Sciences of the CR [IAA100190804]; Ministry of Education of Czech Republic [6840770010]info:eu-repo/semantics/publishedVersio

Moment Closure - A Brief Review

Moment closure methods appear in myriad scientific disciplines in the modelling of complex systems. The goal is to achieve a closed form of a large, usually even infinite, set of coupled differential (or difference) equations. Each equation describes the evolution of one "moment", a suitable coarse-grained quantity computable from the full state space. If the system is too large for analytical and/or numerical methods, then one aims to reduce it by finding a moment closure relation expressing "higher-order moments" in terms of "lower-order moments". In this brief review, we focus on highlighting how moment closure methods occur in different contexts. We also conjecture via a geometric explanation why it has been difficult to rigorously justify many moment closure approximations although they work very well in practice.Comment: short survey paper (max 20 pages) for a broad audience in mathematics, physics, chemistry and quantitative biolog

Overview of mathematical approaches used to model bacterial chemotaxis II: bacterial populations

We review the application of mathematical modeling to understanding the behavior of populations of chemotactic bacteria. The application of continuum mathematical models, in particular generalized Keller–Segel models, is discussed along with attempts to incorporate the microscale (individual) behavior on the macroscale, modeling the interaction between different species of bacteria, the interaction of bacteria with their environment, and methods used to obtain experimentally verified parameter values. We allude briefly to the role of modeling pattern formation in understanding collective behavior within bacterial populations. Various aspects of each model are discussed and areas for possible future research are postulated

SBMLsqueezer: A CellDesigner plug-in to generate kinetic rate equations for biochemical networks

<p>Abstract</p> <p>Background</p> <p>The development of complex biochemical models has been facilitated through the standardization of machine-readable representations like SBML (Systems Biology Markup Language). This effort is accompanied by the ongoing development of the human-readable diagrammatic representation SBGN (Systems Biology Graphical Notation). The graphical SBML editor CellDesigner allows direct translation of SBGN into SBML, and vice versa. For the assignment of kinetic rate laws, however, this process is not straightforward, as it often requires manual assembly and specific knowledge of kinetic equations.</p> <p>Results</p> <p>SBMLsqueezer facilitates exactly this modeling step via automated equation generation, overcoming the highly error-prone and cumbersome process of manually assigning kinetic equations. For each reaction the kinetic equation is derived from the stoichiometry, the participating species (e.g., proteins, mRNA or simple molecules) as well as the regulatory relations (activation, inhibition or other modulations) of the SBGN diagram. Such information allows distinctions between, for example, translation, phosphorylation or state transitions. The types of kinetics considered are numerous, for instance generalized mass-action, Hill, convenience and several Michaelis-Menten-based kinetics, each including activation and inhibition. These kinetics allow SBMLsqueezer to cover metabolic, gene regulatory, signal transduction and mixed networks. Whenever multiple kinetics are applicable to one reaction, parameter settings allow for user-defined specifications. After invoking SBMLsqueezer, the kinetic formulas are generated and assigned to the model, which can then be simulated in CellDesigner or with external ODE solvers. Furthermore, the equations can be exported to SBML, LaTeX or plain text format.</p> <p>Conclusion</p> <p>SBMLsqueezer considers the annotation of all participating reactants, products and regulators when generating rate laws for reactions. Thus, for each reaction, only applicable kinetic formulas are considered. This modeling scheme creates kinetics in accordance with the diagrammatic representation. In contrast most previously published tools have relied on the stoichiometry and generic modulators of a reaction, thus ignoring and potentially conflicting with the information expressed through the process diagram. Additional material and the source code can be found at the project homepage (URL found in the Availability and requirements section).</p