The electrostatic energy can be computed, at least the difference between two conformations:
the charge at one place creates an electrostatic potential everywhere. Where the other(s) charge is, multiply by this charge (normally less than q) to get an energy. Move the other charge between both conformations, get an energy difference.
Then, convert the electrostatic energy into what thermodynamic value you want.
I don't believe the electrostatic energy is defined by the sole electric dipole, because moving 1*q over 2nm to create the dipole takes less energy than moving 2*q over 1nm, and because the same charge spread over a bigger group of atoms needs less energy - telling why a phenyl can more easily accept or give some charge.