reaction.participants.append(metabolites['amp'].species_coefficients.get_or_create(
                    coefficient=ntp_count['A']))