for rnap_kb in cell.observables.get_one(id='rna_polymerase_obs').expression.species:
                rnap_species_type_model = model.species_types.get_one(id=rnap_kb.species_type.id)
                rnap_model = rnap_species_type_model.species.get_one(compartment=cytosol)

                reaction.participants.add(rnap_model.species_coefficients.get_or_create(coefficient=-1))