In this work we have combined a Wang-Landau sampling scheme [F. Wang and D. Landau, Phys. Rev. Lett. 86, 2050 (2001)] with an expanded ensemble formalism to yield a simple and powerful method for computing potentials of mean force. The new method is implemented to investigate the mechanical deformation of proteins. Comparisons are made with analytical results for simple model systems such as harmonic springs and Rouse chains. The method is then illustrated on a model 15-residue alanine molecule in an implicit solvent. Results for mechanical unfolding of this oligopeptide are compared to those of steered molecular dynamics calculations. (C) 2004 American Institute of Physics.