Conformational changes modulate macromolecular function by promoting the specific binding of ligands (such as in antigen recognition) or the stabilization of transition states in enzymatic reactions. However, quantitative characterization of the energetics underlying dynamic structural interconversions is still challenging and lacks a unified method. Here, we introduce a novel in silico approach based on the combined use of essential dynamics sampling and nonequilibrium free-energy calculations to obtain quantitative data on conformational energy landscapes. This technique allows the unbiased investigation of highly complex rearrangements, and does not require the crucial definition of user-defined collective variables. We show that free-energy values derived from profiles connecting the unliganded and ligand-bound X-ray structures of a bacterial nucleoside hydrolase match the experimental binding constant. This approach also provides first evidence for a rate-limiting character of the conformational transition in this enzyme, and an unexpected role of the protonation state of a single residue in regulating substrate binding and product release.