The solvation energy was calculated using the free energy perturbation method implemented in the GROMACS packages. The system consisted of one solute and 512 To circumvent numerical instability in the calculation of LJ interactions, a softcore LJ potential function was employed between the solute and solvent, with the corresponding parameter alpha and sigma taking values of 0.5 and 0.3, respectively.