We introduce a new approach to model the behavior of neuronal signal transduction networks using stochastic differential equations. We present first a mathematical formulation for a stochastic model of protein kinase C pathway. Different kinds of numerical integration methods, including the explicit and implicit Euler