In a recent article published by Bastug and Kuyucak [J. Chem. Phys.126, 105103 (2007)] investigated the microscopic factors affecting double ion occupancy in the gramicidin channel. The analysis relied largely on the one-dimensional potential of mean force of ions along the axis of the channel (the so-called free energy profile of the ion along the channel axis), as well as on the calculation of the equilibrium association constant of the ions in the channel binding sites. It is the purpose of this communication to clarify this issue.