In a recent paper, Lucco Castello et al. (arXiv:2107.03537) provided an accurate parameterization of classical one-component plasma bridge functions that was embedded in a novel dielectric scheme for strongly coupled electron liquids. Here, this approach is rigorously formulated, its set of equations is formally derived, and its numerical algorithm is scrutinized. A systematic comparison with available and new path integral Monte Carlo simulations reveals a rather unprecedented agreement especially in terms of the interaction energy and the long wavelength limit of the static local field correction.