The Schrodinger equation incorporating the long-range Coulomb potential takes the form of a Fredholm equation whose kernel is singular on its diagonal when represented by a basis bearing a continuum of states, such as in a Fourier-Bessel transform. Several methods have been devised to tackle this difficulty, from simply removing the infinite-range of the Coulomb potential with a screening or cut function to using discretizing schemes which take advantage of the integrable character of Coulomb kernel singularities. However, they have never been tested in the context of Berggren bases, which allow many-body nuclear wave functions to be expanded, with halo or resonant properties within a shell model framework. It is thus the object of this paper to test different discretization schemes of the Coulomb potential kernel in the framework of complex-energy nuclear physics. For that, the Berggren basis expansion of proton states pertaining to the sd-shell arising in the A ~ 20 region, being typically resonant, will be effected. Apart from standard frameworks involving a cut function or analytical integration of singularities, a new method will be presented, which replaces diagonal singularities by finite off-diagonal terms. It will be shown that this methodology surpasses in precision the two former techniques.