Abstract:Background: Stochastic simulation has become a useful tool to both study natural biological systems and design new synthetic ones. By capturing the intrinsic molecular fluctuations of "small" systems, these simulations produce a more accurate picture of single cell dynamics, including interesting phenomena missed by deterministic methods, such as noise-induced oscillations and transitions between stable states. However, the computational cost of the original stochastic simulation algorithm can be high, motivating the use of hybrid stochastic methods. Hybrid stochastic methods partition the system into multiple subsets and describe each subset as a different representation, such as a jump Markov, Poisson, continuous Markov, or deterministic process. By applying valid approximations and self-consistently merging disparate descriptions, a method can be considerably faster, while retaining accuracy. In this paper, we describe Hy3S, a collection of multiscale simulation programs. Results: Building on our previous work on developing novel hybrid stochastic algorithms, we have created the Hy3S software package to enable scientists and engineers to both study and design extremely large well-mixed biological systems with many thousands of reactions and chemical species. We have added adaptive stochastic numerical integrators to permit the robust simulation of dynamically stiff biological systems. In addition, Hy3S has many useful features, including embarrassingly parallelized simulations with MPI; special discrete events, such as transcriptional and translation elongation and cell division; mid-simulation perturbations in both the number of molecules of species and reaction kinetic parameters; combinatorial variation of both initial conditions and kinetic parameters to enable sensitivity analysis; use of NetCDF optimized binary format to quickly read and write large datasets; and a simple graphical user interface, written in Matlab, to help users create biological systems and analyze data. We demonstrate the accuracy and efficiency of Hy3S with examples, including a large-scale system benchmark and a complex bistable biochemical network with positive feedback. The software itself is open-sourced under the GPL license and is modular, allowing users to modify it for their own purposes. Conclusion: Hy3S is a powerful suite of simulation programs for simulating the stochastic dynamics of networks of biochemical reactions. Its first public version enables computational biologists to more efficiently investigate the dynamics of realistic biological systems.

Snoopy’s Hybrid Simulator: a Tool to Construct and Simulate Hybrid Biological Models

Biochemical Reaction Networks with Fuzzy Kinetic Parameters in Snoopy

Petri Net Models for the Semi-Automatic Construction of Large Scale Biological Networks

Additional File 3 of Snoopy’s Hybrid Simulator: a Tool to Construct and Simulate Hybrid Biological Models

Hybrid Simulation of Biochemical Systems Using Hybrid Adaptive Petri Nets

A Systematic Petri Net Approach for Multiple-Scale Modeling and Simulation of Biochemical Processes

Multiscale Hy3S: hybrid stochastic simulation for supercomputers

An advanced environment for hybrid modeling of biological systems based on modelica

Snoopy --- a unifying petri net tool

A Workflow for the Hybrid Modelling and Simulation of Multi-Timescale Biological Systems

Modelling and Simulation of Metabolic Networks Petri Nets Approach and Perspective

Colouring fuzziness for systems biology

BioNetSim: a Petri Net-Based Modeling Tool for Simulations of Biochemical Processes

Modeling The Krebs Cycle Based On Hybrid Petri Net

Exact hybrid particle/population simulation of rule-based models of biochemical systems

DynaSim: A MATLAB Toolbox for Neural Modeling and Simulation

Colored Petri Nets for Multiscale Systems Biology - Current Modeling and Analysis Capabilities in Snoopy

Coloured Hybrid Petri Nets: An adaptable modelling approach for multi-scale biological networks

Constructing Hybrid Software Process Simulation Models.

Hybrid Models and Biological Model Reduction with PyDSTool

Computer simulation for biological dynamic networks