There exist many deterministic models for signaling pathways in systems biology. However the models do not consider the stochastic properties of the pathways, which means the models fit well with experimental data in certain situations but poorly in others. Incorporating stochasticity into deterministic models is one way to handle this problem. In this paper the way is used to produce stochastic models based on the deterministic differential equations for the published extracellular signal-regulated kinase (ERK) pathway. We consider strong convergence and stability of the numerical approximations for the stochastic models.