for participant in rxn.participants:
                part = participant.species
                if participant.coefficient < 0:
                    # reactant
                    digraph.add_edge(part, rxn)