Numerical simulations often involve 3D objects with spherical topology, e.g. rigid particles, drops, vesicles. When the underlying numerical method is based on boundary integral equations, standard quadrature rules can yield large errors in computing the layer potentials if the distance between the objects is too small and the associated integrals become nearly singular. I will present numerical and analytical approaches to efficiently evaluate the quadrature error estimates for these situations.