reaction.participants.append(metabolites['gmp'].species_coefficients.get_or_create(
                    coefficient=ntp_count['G']))