Abstract:Background: The temporal progression of many fundamental processes in cells and organisms, including homeostasis, differentiation and development, are governed by gene regulatory networks (GRNs). GRNs balance fluctuations in the output of their genes, which trace back to the stochasticity of molecular interactions. Although highly desirable to understand life processes, predicting the temporal progression of gene products within a GRN is challenging when considering stochastic events such as transcription factor-DNA interactions or protein production and degradation. Results: We report a method to simulate and infer GRNs including genes and biochemical reactions at molecular detail. In our approach, we consider each network element to be isolated from other elements during small time intervals, after which we synchronize molecule numbers across all network elements. Thereby, the temporal behaviour of network elements is decoupled and can be treated by local stochastic or deterministic solutions. We demonstrate the working principle of this modular approach with a repressive gene cascade comprising four genes. By considering a deterministic time evolution within each time interval for all elements, our method approaches the solution of the system of deterministic differential equations associated with the GRN. By allowing genes to stochastically switch between on and off states or by considering stochastic production of gene outputs, we are able to include increasing levels of stochastic detail and approximate the solution of a Gillespie simulation. Thereby, CaiNet is able to reproduce noise-induced bi-stability and oscillations in dynamically complex GRNs. Notably, our modular approach further allows for a simple consideration of deterministic delays. We further infer relevant regulatory connections and steady-state parameters of a GRN of up to ten genes from steady-state measurements by identifying each gene of the network with a single perceptron in an artificial neuronal network and using a gradient decent method originally designed to train recurrent neural networks. To facilitate setting up GRNs and using our simulation and inference method, we provide a fast computer-aided interactive network simulation environment, CaiNet. Conclusion: We developed a method to simulate GRNs at molecular detail and to infer the topology and steady-state parameters of GRNs. Our method and associated user-friendly framework CaiNet should prove helpful to analyze or predict the temporal progression of reaction networks or GRNs in cellular and organismic biology. CaiNet is freely available at https://gitlab.com/GebhardtLab/CaiNet .

Exact derivation and practical application of a hybrid stochastic simulation algorithm for large gene regulatory networks

Exact Parallelization of the Stochastic Simulation Algorithm for Scalable Simulation of Large Biochemical Networks

A Scalable Moment-Closure Approximation for Large-Scale Biochemical Reaction Networks

Stochastic Differential Equations for Practical Simulation of Gene Circuits

Multiscale Hy3S: hybrid stochastic simulation for supercomputers

Constant-complexity Stochastic Simulation Algorithm with Optimal Binning

From Homogeneity to Heterogeneity: Refining Stochastic Simulations of Gene Regulation

Holimap: an accurate and efficient method for solving stochastic gene network dynamics

Periodic synchronization of isolated network elements facilitates simulating and inferring gene regulatory networks including stochastic molecular kinetics

Inferring Gene Regulatory Networks Based on a Hybrid Parallel Genetic Algorithm and the Threshold Restriction Method

Effective computational methods for hybrid stochastic gene networks

Coarse-graining stochastic biochemical networks: quasi-stationary approximation and fast simulations using a stochastic path integral technique

A multiscale compartment-based model of stochastic gene regulatory networks using hitting-time analysis

Green's Function Reaction Dynamics: a new approach to simulate biochemical networks at the particle level and in time and space

A Hierarchical Exact Accelerated Stochastic Simulation Algorithm

BIOCOMP ’ 11 Stability Analysis of Hybrid Stochastic Gene Regulatory Networks

High Dimensional ODEs Coupled with Mixed-Effects Modeling Techniques for Dynamic Gene Regulatory Network Identification

Coarse-Grained Stochastic Particle-based Reaction-Diffusion Simulation Algorithm

Stochastic Simulation of Biomolecular Networks in Dynamic Environments

Experimental design trade-offs for gene regulatory network inference: an in silico study of the yeast Saccharomyces cerevisiae cell cycle

Inferring Large-Scale Gene Regulatory Networks Using a Randomized Algorithm Based on Singular Value Decomposition.