The usual approach to the study of signaling pathways in biological systems is to assume that high numbers of cells and of perfectly mixed molecules within cells are involved. To study the temporal evolution of the system averaged over the cell population, ordinary differential equations are usually used. However, this approach has been shown to be inadequate if few copies of molecules and/or cells are present. In such situation, a stochastic or a hybrid stochastic/deterministic approach needs to be used. Moreover, considering a perfectly mixed system in cases where spatial effects are present can be an over-simplifying assumption. This can be corrected by adding diffusion terms to the ordinary differential equations describing chemical reactions and proliferation kinetics. However, there exist cases in which both stochastic and spatial effects have to be considered. We study the relevance of differential equations, stochastic Gillespie algorithm, and deterministic and stochastic reaction-diffusion models for the study of important biological processes, such as viral infection and early carcinogenesis. To that end we have developed two optimized libraries of C functions for R (r-project.org) to simulate biological systems using Petri Nets, in a pure deterministic, pure stochastic, or hybrid deterministic/stochastic fashion, with and without spatial effects. We discuss our findings in the terms of specific biological systems including signaling in innate immune response, early carcinogenesis and spatial spread of viral infection.