Motivation: Many biochemical systems require stochastic descriptions. Unfortunately these can only be solved for the simplest cases and their direct simulation can become prohibitively expensive, precluding thorough analysis. As an alternative, moment closure approximation methods generate equations for the time-evolution of the system’s moments and apply a closure ansatz to obtain a closed set of differential equations; that can become the basis for the deterministic analysis of the moments of the outputs of stochastic systems.Results: We present a free, user-friendly tool implementing an efficient moment expansion approximation with parametric closures that integrates well with the IPython interactive environment. Our package enables the analysis of complex stochastic systems without any constraints on the number of species and moments studied and the type of rate laws in the system. In addition to the approximation method our package provides numerous tools to help non-expert users in stochastic analysis.Availability and implementation: https://github.com/theosysbio/meansContacts: m.stumpf@imperial.ac.uk or e.lakatos13@imperial.ac.ukSupplementary information: Supplementary data are available at Bioinformatics online.

Chemotaxis is vital cellular movement in response to environmental chemicals. Unlike the canonical chemotactic pathway in Escherichia coli , Rhodobacter sphaeroides has both transmembrane and cytoplasmic sensory clusters, with the latter possibly interacting with essential components in the electron transport system. However, the effect of the cytoplasmic sensor and the mechanism of signal integration from both sensory clusters remain unclear. Based on a minimal model of the chemotaxis pathway in this species, we show that signal integration at the motor level produces realistic chemotactic behaviour in line with experimental observations. Our model also suggests that the core pathway of R. sphaeroides , at least its ancestor, may represent a metabolism-dependent selective stopping strategy, which alone can steer cells to favourable environments. Our results not only clarify the potential roles of the two sensory clusters but also put in question the current definitions of attractants and repellents.