Efficient stochastic simulation of gene regulatory networks using hybrid models of transcriptional bursting