The free energy profiles were computed from the distribution densities of pore waters and pore Na + s. Gaussian kernel density estimates were generated with the gaussian kde function under the scipy.stats module in SCIPY from the Na + and water positions along the z-axis of all the conformations. The densities were then converted to free energy: