| 1. |
|
|
| 2. |
- Coulier, Adrien, et al.
(författare)
-
A multiscale compartment-based model of stochastic gene regulatory networks using hitting-time analysis
- 2021
-
Ingår i: Journal of Chemical Physics. - : American Institute of Physics (AIP). - 0021-9606 .- 1089-7690. ; 154:18
-
Tidskriftsartikel (refereegranskat)abstract
- Spatial stochastic models of single cell kinetics are capable of capturing both fluctuations in molecular numbers and the spatial dependencies of the key steps of intracellular regulatory networks. The spatial stochastic model can be simulated both on a detailed microscopic level using particle tracking and on a mesoscopic level using the reaction–diffusion master equation. However, despite substantial progress on simulation efficiency for spatial models in the last years, the computational cost quickly becomes prohibitively expensive for tasks that require repeated simulation of thousands or millions of realizations of the model. This limits the use of spatial models in applications such as multicellular simulations, likelihood-free parameter inference, and robustness analysis. Further approximation of the spatial dynamics is needed to accelerate such computational engineering tasks. We here propose a multiscale model where a compartment-based model approximates a detailed spatial stochastic model. The compartment model is constructed via a first-exit time analysis on the spatial model, thus capturing critical spatial aspects of the fine-grained simulations, at a cost close to the simple well-mixed model. We apply the multiscale model to a canonical model of negative-feedback gene regulation, assess its accuracy over a range of parameters, and demonstrate that the approximation can yield substantial speedups for likelihood-free parameter inference.
|
|
| 3. |
|
|
| 4. |
- Hellander, Andreas, et al.
(författare)
-
Coupled mesoscopic and microscopic simulation of stochastic reaction-diffusion processes in mixed dimensions
- 2011
-
Rapport (övrigt vetenskapligt/konstnärligt)abstract
- We present a new simulation algorithm that allows for dynamic switching between a mesoscopic and a microscopic modeling framework for stochastic reaction-diffusion kinetics. The more expensive and more accurate microscopic model is used only for those species and in those regions in space where there is reason to believe that a microscopic model is needed to capture the dynamics correctly. The microscopic algorithm is extended to simulation on curved surfaces in order to model reaction and diffusion on membranes. The accuracy of the method on and near a spherical membrane is analyzed and evaluated in a numerical experiment. Two biologically motivated examples are simulated in which the need for microscopic simulation of parts of the system arises for different reasons. First, we apply the method to a model of the phosphorylation reactions in a MAPK signaling cascade where microscale methods are necessary to resolve fast rebinding events. Then a model is considered for transport of a species over a membrane coupled to reactions in the bulk. The new algorithm attains an accuracy similar to a full microscopic simulation by handling critical interactions on the microscale, but at a significantly reduced cost by using the mesoscale framework for most parts of the biological model.
|
|
| 5. |
|
|
| 6. |
|
|
| 7. |
- Hellander, Stefan, et al.
(författare)
-
Mesoscopic-microscopic spatial stochastic simulation with automatic system partitioning
- 2017
-
Ingår i: Journal of Chemical Physics. - : AIP Publishing. - 0021-9606 .- 1089-7690. ; 147:23
-
Tidskriftsartikel (refereegranskat)abstract
- The reaction-diffusion master equation (RDME) is a model that allows for efficient on-lattice simulation of spatially resolved stochastic chemical kinetics. Compared to off-lattice hard-sphere simulations with Brownian dynamics or Green's function reaction dynamics, the RDME can be orders of magnitude faster if the lattice spacing can be chosen coarse enough. However, strongly diffusion-controlled reactions mandate a very fine mesh resolution for acceptable accuracy. It is common that reactions in the same model differ in their degree of diffusion control and therefore require different degrees of mesh resolution. This renders mesoscopic simulation inefficient for systems with multiscale properties. Mesoscopic-microscopic hybrid methods address this problem by resolving the most challenging reactions with a microscale, off-lattice simulation. However, all methods to date require manual partitioning of a system, effectively limiting their usefulness as "black-box" simulation codes. In this paper, we propose a hybrid simulation algorithm with automatic system partitioning based on indirect a priori error estimates. We demonstrate the accuracy and efficiency of the method on models of diffusion-controlled networks in 3D.
|
|
| 8. |
|
|
| 9. |
|
|
| 10. |
- Bani-Hashemian, Hossein, et al.
(författare)
-
Efficient sampling in event-driven algorithms for reaction-diffusion processes
- 2011
-
Rapport (övrigt vetenskapligt/konstnärligt)abstract
- In event-driven algorithms for simulation of diffusing, colliding, and reacting particles, new positions and events are sampled from the cumulative distribution function (CDF) of a probability distribution. The distribution is sampled frequently and it is important for the efficiency of the algorithm that the sampling is fast. The CDF is known analytically or computed numerically. Analytical formulas are sometimes rather complicated making them difficult to evaluate. The CDF may be stored in a table for interpolation or computed directly when it is needed. Different alternatives are compared for chemically reacting molecules moving by Brownian diffusion in two and three dimensions. The best strategy depends on the dimension of the problem, the length of the time interval, the density of the particles, and the number of different reactions.
|
|
