An official website of the United States government
Models of reaction networks within interacting compartments (RNIC) are a generalization of stochastic reaction networks. It is most natural to think of the interacting compartments as “cells” that can appear, degrade, split, and even merge, with each cell containing an evolving copy of the underlying stochastic reaction network. Such models have a number of parameters, including those associated with the internal chemical model and those associated with the compartment interactions, and it is natural to want efficient computational methods for the numerical estimation of sensitivities of model statistics with respect to these parameters. Motivated by the extensive work on computational methods for parametric sensitivity analysis in the context of stochastic reaction networks over the past few decades, we provide a number of methods in the basic RNIC setting. Provided methods include the (unbiased) Girsanov transformation method (also called the likelihood ratio method) and a number of coupling methods for the implementation of finite differences, each motivated by methods from previous work related to stochastic reaction networks. We provide several numerical examples comparing the various methods in the new setting. We find that the relative performance of each method is in line with its analog in the “standard” stochastic reaction network setting. We have made all the MATLAB codes used to implement the various methods freely available for download.
more »
« less
Stochastic Reaction Networks Within Interacting Compartments
Stochastic reaction networks, which are usually modeled as continuous-time Markov chains on $$\ZZ^d_{\ge 0}$$, and simulated via a version of the ``Gillespie algorithm,'' have proven to be a useful tool for the understanding of processes, chemical and otherwise, in homogeneous environments. There are multiple avenues for generalizing away from the assumption that the environment is homogeneous, with the proper modeling choice dependent upon the context of the problem being considered. One such generalization was recently introduced in \cite{Duso_Zechner_2020}, where the proposed model includes a varying number of interacting compartments, or cells, each of which contains an evolving copy of the stochastic reaction system. The novelty of the model is that these compartments also interact via the merging of two compartments (including their contents), the splitting of one compartment into two, and the appearance and destruction of compartments. In this paper we begin a systematic exploration of the mathematical properties of this model. We (i) obtain basic/foundational results pertaining to explosivity, transience, recurrence, and positive recurrence of the model, (ii) explore a number of examples demonstrating some possible non-intuitive behaviors of the model, and (iii) identify the limiting distribution of the model in a special case that generalizes three formulas from an example in \cite{Duso_Zechner_2020}.
more »
« less
- Award ID(s):
- 2051498
- PAR ID:
- 10533274
- Publisher / Repository:
- Springer
- Date Published:
- Journal Name:
- Bulletin of Mathematical Biology
- Volume:
- 85
- Issue:
- 10
- ISSN:
- 0092-8240
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Spreading and mixing are complementary processes that promote reaction of two reactive aqueous solutes present in contiguous plumes in groundwater. Spreading reconfigures the plume geometry, elongating the interface between the plumes, while mixing increases the volume of aquifer occupied by each plume, bringing the solute molecules together to react. Since reaction only occurs where the two solute plumes are in contact with each other, local mechanisms that drive flow and transport near the interface between the plumes control the amount of reaction. This work uses local characteristics of the plumes and the flow field near the plume interface to analyze the relative contributions of pore‐scale mixing and mechanical dispersion to instantaneous, irreversible, bimolecular reaction in a homogeneous aquifer with active spreading caused by radial flow from a well. Two solutes are introduced in sequence at the well, creating concentric circular plumes. We allow for incomplete mixing of the solutes in the pore space, by modeling the pore space as a segregated compartment and a mixed compartment with first‐order mass transfer between the two compartments. We develop semi‐analytical expressions for concentrations of the solutes in both compartments. We found that the relative contribution of mechanical dispersion to reaction increases over time and also increases due to increases in the Peclet number, in the relative source concentration of the chasing solute, and in the mass transfer rate from the segregated compartment to the mixed compartment of the pore space.more » « less
-
Compartmental exchange regulates steady states and stochastic switching of a phosphorylation networkThe phosphoregulation of proteins with multiple phosphorylation sites is governed by biochemical reaction networks that can exhibit multistable behavior. However, the behavior of such networks is typically studied in a single reaction volume, while cells are spatially organized into compartments that can exchange proteins. In this work, we use stochastic simulations to study the impact of compartmentalization on a two-site phosphorylation network. We characterize steady states and fluctuation-driven transitions between them as a function of the rate of protein exchange between two compartments. Surprisingly, the average time spent in a state before stochastically switching to another depends nonmonotonically on the protein exchange rate, with the most frequent switching occurring at intermediate exchange rates. At sufficiently small exchange rates, the state of the system and mean switching time are controlled largely by fluctuations in the balance of enzymes in each compartment. This leads to negatively correlated states in the compartments. For large exchange rates, the two compartments behave as a single effective compartment. However, when the compartmental volumes are unequal, the behavior differs from a single compartment with the same total volume. These results demonstrate that exchange of proteins between distinct compartments can regulate the emergent behavior of a common signaling motif.more » « less
-
Stochastic second-order methods are known to achieve fast local convergence in strongly convex optimization by relying on noisy Hessian estimates to precondition the gradient. Yet, most of these methods achieve superlinear convergence only when the stochastic Hessian noise diminishes, requiring an increase in the per-iteration cost as time progresses. Recent work in \cite{na2022hessian} addressed this issue via a Hessian averaging scheme that achieves a superlinear convergence rate without increasing the per-iteration cost. However, the considered method exhibits a slow global convergence rate, requiring up to ~O(κ^2) iterations to reach the superlinear rate of ~O((1/t)^{t/2}), where κ is the problem's condition number. In this paper, we propose a novel stochastic Newton proximal extragradient method that significantly improves these bounds, achieving a faster global linear rate and reaching the same fast superlinear rate in ~O(κ) iterations. We achieve this by developing a novel extension of the Hybrid Proximal Extragradient (HPE) framework, which simultaneously achieves fast global and local convergence rates for strongly convex functions with access to a noisy Hessian oracle.more » « less
-
The basic reproduction number R0 is a concept which originated in population dynamics, mathematical epidemiology, and ecology and is closely related to the mean number of children in branching processes (reflecting the fact that the phenomena of interest are well approximated via branching processes, at their inception). Despite the very extensive literature around R0 for deterministic epidemic models, we believe there are still aspects which are not fully understood. Foremost is the fact that R0 is not a function of the original ODE model, unless we also include in it a certain (F,V) gradient decomposition, which is not unique. This is related to the specification of the “infected compartments”, which is also not unique. A second interesting question is whether the extinction probabilities of the natural continuous time Markovian chain approximation of an ODE model around boundary points (disease-free equilibrium and invasion points) are also related to the (F,V) gradient decomposition. We offer below several new contributions to the literature: (1) A universal algorithmic definition of a (F,V) gradient decomposition (and hence of the resulting R0). (2) A fixed point equation for the extinction probabilities of a stochastic model associated to a deterministic ODE model, which may be expressed in terms of the (F,V) decomposition. Last but not least, we offer Mathematica scripts and implement them for a large variety of examples, which illustrate that our recipe offers always reasonable results, but that sometimes other reasonable (F,V) decompositions are available as well.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1007/s11538-023-01185-w
