Gene regulatory networks: a coarse-grained, equation-free approach to multiscale computation

Abstract

We present computer-assisted methods for analyzing stochastic models of gene regulatory networks. The main idea that underlies this equation-free analysis is the design and execution of appropriately-initialized short bursts of stochastic simulations; the results of these are processed to estimate coarse-grained quantities of interest, such as mesoscopic transport coefficients. In particular, using a simple model of a genetic toggle switch, we illustrate the computation of an effective free energy and of a state-dependent effective diffusion coefficient that characterize an unavailable effective Fokker-Planck equation. Additionally we illustrate the linking of equation-free techniques with continuation methods for performing a form of stochastic "bifurcation analysis"; estimation of mean switching times in the case of a bistable switch is also implemented in this equation-free context. The accuracy of our methods is tested by direct comparison with long-time stochastic simulations. This type of equation-free analysis appears to be a promising approach to computing features of the long-time, coarse-grained behavior of certain classes of complex stochastic models of gene regulatory networks, circumventing the need for long Monte Carlo simulations.

arXiv comment33 pages, submitted to The Journal of Chemical Physics

References in corpus (7)

• Coarse Molecular Dynamics of a Peptide Fragment: Free Energy, Kinetics, and Long-Time Dynamics Computations
• Coarse Stability and Bifurcation Analysis Using Stochastic Simulators: Kinetic Monte Carlo Examples
• Coarse Brownian Dynamics for Nematic Liquid Crystals: Bifurcation Diagrams via Stochastic Simulation
• Coarse-Grained Kinetic Computations for Rare Events: Application to Micelle Formation
• Coarse Projective kMC Integration: Forward/Reverse Initial and Boundary Value Problems
• Coarse Nonlinear Dynamics and Metastability of Filling-Emptying Transitions: Water in Carbon Nanotubes
• Apparent hysteresis in a driven system with self-organized drag

Cited by in corpus (13)

• Reduction of dynamical biochemical reaction networks in computational biology
• Variable-free exploration of stochastic models: a gene regulatory network example
• Coarse-grained dynamics of an activity bump in a neural field model
• Perfect Sampling of the Master Equation for Gene Regulatory Networks
• Model reduction for stochastic chemical systems with abundant species
• Onset of collective motion in locusts is captured by a minimal model
• Deterministic continutation of stochastic metastable equilibria via Lyapunov equations and ellipsoids
• Ergodic directional switching in mobile insect groups
• Spatially Distributed Stochastic Systems: equation-free and equation-assisted preconditioned computation
• Statistically consistent coarse-grained simulations for critical phenomena in complex networks
• Constrained Approximation of Effective Generators for Multiscale Stochastic Reaction Networks and Application to Conditioned Path Sampling
• Coarse-grained Monte Carlo simulations of the phase transition of Potts model on weighted networks
• Synchronized molecular dynamics method for thin-layer flows of complex fluids