def METHOD_NAME(expression):
    expr_sps = []
    for a in expression.expr.atoms():
        if isinstance(a, pysb.core.Observable):
            sps = a.species
            expr_sps += sps
    return set(expr_sps)