Introducing and benchmarking the accuracy of cayenne: a Python package for stochastic simulations

Biological systems are intrinsically noisy and this noise may determine the qualitative outcome of the system. In the absence of analytical solutions to mathematical models incorporating noise, stochastic simulation algorithms are useful to explore the possible trajectories of these systems. Algorithms used for such stochastic simulations include the Gillespie algorithm and its approximations. In this study we introduce cayenne, an easy to use Python package containing accurate and fast implementations of the Gillespie algorithm (direct method), the tau-leaping algorithm and a tau-adaptive algorithm. We compare the accuracy of cayenne with other stochastic simulation libraries (BioSimulator.jl, GillespieSSA and Tellurium) and find that cayenne offers the best trade-off between accuracy and speed. Additionally, we highlight the importance of performing accuracy tests for stochastic simulation libraries, and hope that it becomes standard practice when developing the same. The cayenne package can be found at https://github.com/Heuro-labs/cayenne while the bench-marks can be found at https://github.com/Heuro-labs/cayenne-benchmarks