Determination of Henry's law constants through transition matrix Monte Carlo simulation