Force-field based grand-canonical Monte Carlo simulations are employed to predict the hydrogen adsorption properties of seven structurally different MOFs. The performance of different parameter sets is assessed by comparison with experimental data, and the capabilities and limitations of the methodology are critically discussed, with a particular emphasis on systems with unsaturated metal sites. In addition to adsorption isotherms and isosteric heats of adsorption, the preferred adsorption sites are obtained from a detailed analysis of the calculated hydrogen density fields. Where possible, these positions are compared to the results of neutron diffraction experiments. This study highlights the capabilities of computational methods to identify the structural features which are most favourable for hydrogen adsorption, providing valuable implications for the synthesis of novel MOFs.