Simulating gene regulatory networks using reaction systems