It is an implicit assumption in almost all papers other than our own that the distribution of RBFEs obtained from different replicas in an ensemble is Gaussian. It is why, even if just 2 or 3 replica simulations are performed, those are deemed sufficient. Paliwal et al. claimed that hydration free energy distributions are Gaussian and deviate from this behaviour at the 95% confidence level only when the ensemble sizes are increased above 110. However, in addition to the fact that their study is based on toy models with simple interactions, the dependence of normality on the ensemble size that they report raises questions as to the reliability of their interpretation of the general nature of free energy distributions. We have repeatedly shown that such distributions deviate from normality for complex protein-ligand systems. Moreover, this is not just confined to ligand-protein interactions but valid for many non-linear systems with long-range interactions. Furthermore, a recent paper of ours demonstrates that the influence of the random number seed used in each MD simulation completely dominates the uncertainty accruing from the uncertainty in the parameters used to perform such simulations. The nature of free energy distributions associated with MD simulations and their important consequences are simply not touched upon by any of the publications discussed in the Introduction.