cyto_deg_rate = cyto_deg_reaction.rate_laws[0].expression._parsed_expression.eval({
                    wc_lang.Species: cyto_species_counts,
                    wc_lang.Compartment: {
                        cytosol.id: cytosol.init_volume.mean * \
                            cytosol.init_density.value}