On the rejection-based algorithm for simulation and analysis of large-scale reaction networks