Endpoint bilinear estimates and applications to the two-dimensional Poisson–Nernst–Planck system