1
|
Morgado G, Nowakowski B, Lemarchand A. Elimination of fast variables in stochastic nonlinear kinetics. Phys Chem Chem Phys 2020; 22:20801-20814. [PMID: 32909563 DOI: 10.1039/d0cp02785e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A reduced chemical scheme involving a small number of variables is often sufficient to account for the deterministic evolution of the concentration of the main species contributing to a reaction. However, its predictions are questionable in small systems used, for example in fluorescence correlation spectroscopy (FCS) or in explosive systems involving strong nonlinearities such as autocatalytic steps. We make precise dynamical criteria defining the validity domain of the quasi-steady-state approximation and the elimination of a fast concentration in deterministic dynamics. Designing two different three-variable models converging toward the same two-variable model, we show that the variances and covariance of the fluctuations of the slow variables are not correctly predicted using the two-variable model, even in the limit of a large system size. The more striking weaknesses of the reduced scheme are figured out in mesoscaled systems containing a small number of molecules. The results of two stochastic approaches are compared and the shortcomings of the Langevin equations with respect to the master equation are pointed out. We conclude that the description of the fluctuations and their coupling with nonlinearities of deterministic dynamics escape reduced chemical schemes.
Collapse
Affiliation(s)
- Gabriel Morgado
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland and Laboratoire de Physique Théorique de la Matière Condensée, Sorbonne Université, CNRS, 4 Place Jussieu, Case Courrier 121, 75252 Paris CEDEX 05, France.
| | - Bogdan Nowakowski
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland and SGGW, Warsaw University of Life Sciences, Nowoursynowska 159, 02-776 Warsaw, Poland
| | - Annie Lemarchand
- Laboratoire de Physique Théorique de la Matière Condensée, Sorbonne Université, CNRS, 4 Place Jussieu, Case Courrier 121, 75252 Paris CEDEX 05, France.
| |
Collapse
|
2
|
Heineken, Tsuchiya and Aris on the mathematical status of the pseudo-steady state hypothesis: A classic from volume 1 of Mathematical Biosciences. Math Biosci 2019; 318:108274. [PMID: 31697965 DOI: 10.1016/j.mbs.2019.108274] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2019] [Revised: 10/26/2019] [Accepted: 10/26/2019] [Indexed: 01/25/2023]
Abstract
Volume 1, Issue 1 of Mathematical Biosciences was the venue for a now-classic paper on the application of singular perturbation theory in enzyme kinetics, "On the mathematical status of the pseudo-steady state hypothesis of biochemical kinetics" by F. G. Heineken, H. M. Tsuchiya and R. Aris. More than 50 years have passed, and yet this paper continues to be studied and mined for insights. This perspective discusses both the strengths and weaknesses of the work presented in this paper. For many, the justification of the pseudo-steady-state approximation using singular perturbation theory is the main achievement of this paper. However, there is so much more material here, which laid the foundation for a great deal of research in mathematical biochemistry in the intervening decades. The parameterization of the equations, construction of the first-order uniform singular-perturbation solution, and an attempt to apply similar principles to the pseudo-equilibrium approximation are discussed in particular detail.
Collapse
|
3
|
Vigelius M, Meyer B, Pascoe G. Multiscale modelling and analysis of collective decision making in swarm robotics. PLoS One 2014; 9:e111542. [PMID: 25369026 PMCID: PMC4219724 DOI: 10.1371/journal.pone.0111542] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2014] [Accepted: 10/05/2014] [Indexed: 11/18/2022] Open
Abstract
We present a unified approach to describing certain types of collective decision making in swarm robotics that bridges from a microscopic individual-based description to aggregate properties. Our approach encompasses robot swarm experiments, microscopic and probabilistic macroscopic-discrete simulations as well as an analytic mathematical model. Following up on previous work, we identify the symmetry parameter, a measure of the progress of the swarm towards a decision, as a fundamental integrated swarm property and formulate its time evolution as a continuous-time Markov process. Contrary to previous work, which justified this approach only empirically and a posteriori, we justify it from first principles and derive hard limits on the parameter regime in which it is applicable.
Collapse
Affiliation(s)
- Matthias Vigelius
- FIT Centre for Research in Intelligent Systems, Monash University, Melbourne, Australia
| | - Bernd Meyer
- FIT Centre for Research in Intelligent Systems, Monash University, Melbourne, Australia
- * E-mail:
| | - Geoffrey Pascoe
- FIT Centre for Research in Intelligent Systems, Monash University, Melbourne, Australia
| |
Collapse
|
4
|
Dóka É, Lente G. Stochastic mapping of the Michaelis-Menten mechanism. J Chem Phys 2012; 136:054111. [DOI: 10.1063/1.3681942] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
5
|
Dóka É, Lente G. Mechanism-Based Chemical Understanding of Chiral Symmetry Breaking in the Soai Reaction. A Combined Probabilistic and Deterministic Description of Chemical Reactions. J Am Chem Soc 2011; 133:17878-81. [DOI: 10.1021/ja207408y] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Éva Dóka
- Department of Inorganic and Analytical Chemistry, University of Debrecen, POB 21, Debrecen H-4010, Hungary
| | - Gábor Lente
- Department of Inorganic and Analytical Chemistry, University of Debrecen, POB 21, Debrecen H-4010, Hungary
| |
Collapse
|
6
|
Gómez-Uribe CA, Verghese GC, Tzafriri AR. Enhanced identification and exploitation of time scales for model reduction in stochastic chemical kinetics. J Chem Phys 2009; 129:244112. [PMID: 19123500 DOI: 10.1063/1.3050350] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Widely different time scales are common in systems of chemical reactions and can be exploited to obtain reduced models applicable to the time scales of interest. These reduced models enable more efficient computation and simplify analysis. A classic example is the irreversible enzymatic reaction, for which separation of time scales in a deterministic mass action kinetics model results in approximate rate laws for the slow dynamics, such as that of Michaelis-Menten. Recently, several methods have been developed for separation of slow and fast time scales in chemical master equation (CME) descriptions of stochastic chemical kinetics, yielding separate reduced CMEs for the slow variables and the fast variables. The paper begins by systematizing the preliminary step of identifying slow and fast variables in a chemical system from a specification of the slow and fast reactions in the system. The authors then present an enhanced time-scale-separation method that can extend the validity and improve the accuracy of existing methods by better accounting for slow reactions when equilibrating the fast subsystem. The resulting method is particularly accurate in systems such as enzymatic and protein interaction networks, where the rates of the slow reactions that modify the slow variables are not a function of the slow variables. The authors apply their methodology to the case of an irreversible enzymatic reaction and show that the resulting improvements in accuracy and validity are analogous to those obtained in the deterministic case by using the total quasi-steady-state approximation rather than the classical Michaelis-Menten. The other main contribution of this paper is to show how mass fluctuation kinetics models, which give approximate evolution equations for the means, variances, and covariances of the concentrations in a chemical system, can feed into time-scale-separation methods at a variety of stages.
Collapse
|
7
|
Pigolotti S, Vulpiani A. Coarse graining of master equations with fast and slow states. J Chem Phys 2008; 128:154114. [DOI: 10.1063/1.2907242] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
8
|
Fiasconaro A, Spagnolo B, Ochab-Marcinek A, Gudowska-Nowak E. Co-occurrence of resonant activation and noise-enhanced stability in a model of cancer growth in the presence of immune response. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2006; 74:041904. [PMID: 17155093 DOI: 10.1103/physreve.74.041904] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/01/2005] [Revised: 08/18/2006] [Indexed: 05/12/2023]
Abstract
We investigate a stochastic version of a simple enzymatic reaction which follows the generic Michaelis-Menten kinetics. At sufficiently high concentrations of reacting species, that represent here populations of cells involved in cancerous proliferation and cytotoxic response of the immune system, the overall kinetics can be approximated by a one-dimensional overdamped Langevin equation. The modulating activity of the immune response is here modeled as a dichotomous random process of the relative rate of neoplastic cell destruction. We discuss physical aspects of environmental noises acting in such a system, pointing out the possibility of coexistence of dynamical regimes where noise-enhanced stability and resonant activation phenomena can be observed together. We explain the underlying mechanisms by analyzing the behavior of the variance of first passage times as a function of the noise intensity.
Collapse
Affiliation(s)
- Alessandro Fiasconaro
- Dipartimento di Fisica e Tecnologie Relative and CNISM, Group of Interdisciplinary Physics, Università di Palermo, Viale delle Scienze, I-90128 Palermo, Italy
| | | | | | | |
Collapse
|
9
|
Pineda M, Imbihl R, Schimansky-Geier L, Zülicke C. Theoretical analysis of internal fluctuations and bistability in CO oxidation on nanoscale surfaces. J Chem Phys 2006; 124:044701. [PMID: 16460194 DOI: 10.1063/1.2140709] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The bistable CO oxidation on a nanoscale surface is characterized by a limited number of reacting molecules on the catalytic area. Internal fluctuations due to finite-size effects are studied by the master equation with a Langmuir-Hinshelwood mechanism for CO oxidation. Analytical solutions can be found in a reduced one-component model after the adiabatic elimination of one variable which in our case is the oxygen coverage. It is shown that near the critical point, with decreasing surface area, one cannot distinguish between two macroscopically stable stationary states. This is a consequence of the large fluctuations in the coverage which occur on a fast time scale. Under these conditions, the transition times between the macroscopic states also are no longer separated from the short-time scale of the coverage fluctuations as is the case for large surface areas and far away from the critical point. The corresponding stationary solutions of the probability distribution and the mean first passage times calculated in the reduced model are supported by numerics of the full two-component model.
Collapse
Affiliation(s)
- M Pineda
- Institut für Physikalische Chemie und Elektrochemie, Universität Hannover, Callinstrasse 3-3a, D-30167 Hannover, Germany.
| | | | | | | |
Collapse
|
10
|
Roussel MR, Zhu R. Reducing a chemical master equation by invariant manifold methods. J Chem Phys 2004; 121:8716-30. [PMID: 15527335 DOI: 10.1063/1.1802495] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We study methods for reducing chemical master equations using the Michaelis-Menten mechanism as an example. The master equation consists of a set of linear ordinary differential equations whose variables are probabilities that the realizable states exist. For a master equation with s(0) initial substrate molecules and e(0) initial enzyme molecules, the manifold can be parametrized by s(0) of the probability variables. Fraser's functional iteration method is found to be difficult to use for master equations of high dimension. Building on the insights gained from Fraser's method, techniques are developed to produce s(0)-dimensional manifolds of larger systems directly from the eigenvectors. We also develop a simple, but surprisingly effective way to generate initial conditions for the reduced models.
Collapse
Affiliation(s)
- Marc R Roussel
- Department of Chemistry and Biochemistry, University of Lethbridge, Lethbridge, Alberta T1K 3M4, Canada.
| | | |
Collapse
|