Stochasic simulation has limited applicability to whole-cell or multi-cell biochemical systems due to immense number of reaction channels. We have employed parallel exact stochastic simulation in modeling large-scale biochemical systems and tested on parallel shared memory supercomputers. Efficient model subdivision techniques has been utilized to significantly enhance the performance. Our results demonstrated our method is applicable and effective in a variety of biochemical models.