Using data from equilibrium molecular dynamics computer simulations we have built up a catalog of response functions for the Coulomb one-component plasma over a wide range of Γ coupling values, including the strongly coupled Γ>1 liquid regime. We focus on the domain of negative compressibility (Γ>3), where the proper response displays an acausal behavior, implying a modification of the relation between its real and imaginary parts in the Kramers-Kronig relations. We give a description of the details of this acausal feature, in both the frequency and time domains. We show that the viscoelastic pole of the density response function morphs into an imaginary pole in the upper ω half-plane that is responsible for the anomalous behavior of the response in this coupling range. By examining the plasmon dispersion relation through the dielectric response function, rather than via the peaks of the dynamical structure function, we obtain a more reliable representation for the dispersion. We demonstrate that there is an intimate link between the formation of the roton minimum in the dispersion and the negative compressibility of the system. The feasibility of the extension of our analysis to systems with a short-range interaction is explored.