The availability of a huge amount of molecular data concerning various biochemical reactions provoked numerous attempts to study the dynamics of cellular processes by means of kinetic models and computer simulations. Biochemical processes frequently involve small numbers of molecules (e.g. a few molecules of a transcriptional regulator binding to one molecule of a DNA regulatory region). Such reactions are subject to significant stochastic fluctuations. Monte Carlo methods must be employed to study the functional consequences of the fluctuations and simulate processes that cannot be modelled by continuous fluxes of matter. This provides the motivation to develop software dedicated to Monte Carlo simulations of cellular processes with the rigorously proven Gillespie algorithm.