The energy landscape of proteins is thought to have an intricate, corrugated structure. Such roughness should have important consequences on the folding and binding kinetics of proteins, as well as on their equilibrium fluctuations. So far, no direct measurement of protein energy landscape roughness has been made. Here, we combined a recent theory with single-molecule dynamic force spectroscopy experiments to extract the overall energy scale of roughness epsilon for a complex consisting of the small GTPase Ran and the nuclear transport receptor importin-beta. The results gave epsilon > 5k(B)T, indicating a bumpy energy surface, which is consistent with the ability of importin-beta to accommodate multiple conformations and to interact with different, structurally distinct ligands.