reaction.participants.append(metabolites['cmp'].species_coefficients.get_or_create(
                    coefficient=ntp_count['C']))