eval_rate_law = reaction.rate_laws[0].expression._parsed_expression.eval({
                    wc_lang.Species: init_species_counts,
                    wc_lang.Compartment: {
                        transcription_compartment.id: transcription_compartment.init_volume.mean * \
                            transcription_compartment.init_density.value},