Parallelization of tau-leap coarse-grained Monte Carlo simulations on GPUs