reaction.participants.append(metabolites['ump'].species_coefficients.get_or_create(
                    coefficient=ntp_count['U']))