def bounding_box(points):
    """
    Computes the bounding box of a list of 2D points.

    :param points: list of (x, y) tuples, representing the points
    :return: list of four (x, y) tuples, representing the vertices of the bounding box
    """
    xs, ys = zip(*points)
    x_min = min(xs)
    x_max = max(xs)
    y_min = min(ys)
    y_max = max(ys)
    return [(x_min, y_min), (x_max, y_min), (x_max, y_max), (x_min, y_max)]