The one-step perturbation approach is an efficient means to calculate many relative free energies from a common reference compound. Combining lessons learned in previous studies, an application of the method is presented that allows for the calculation of relative binding free energies for structurally rather diverse compounds from only a few simulations. Based on the well known statistical-mechanical perturbation formula, the results do not require any empirical parameters, or training sets, only limited knowledge of the binding characteristics of the ligands suffices to design appropriate reference compounds. Depending on the choice of reference compound, relative free energies of binding rigid ligands to the ligand-binding domain of the estrogen receptor can be obtained that show good agreement with the experimental values. The approach presented here can easily be applied to many rigid ligands, and it should be relatively easy to extend the method to account for ligand flexibility. The free-energy calculations can be straightforwardly parallelized, allowing for an efficient means to understand and predict relative binding free energies. free-energy calculation ͉ molecular dynamics simulation ͉ estrogen receptor ͉ one-step perturbation T he calculation of relative binding free energies for several ligands to a common receptor is of relevance for drug design purposes, and for obtaining a better understanding of the molecular interactions of proteins with small compounds. Despite the increased availability of computational power, calculating free energies from molecular dynamics simulations is still time-consuming, often requiring several extensive simulations to obtain a single relative free energy. The approach in which several free energies can be obtained from a single simulation of a, not necessarily physically meaningful, reference state (1) has been applied successfully in recent years (2-6). The idea behind the method is to simulate a judiciously chosen reference compound R, generating an ensemble of structures that contains conformations representative for several physically relevant compounds. The free-energy difference between any real ligand A and the reference compound can be calculated from the perturbation formula (7)where the angle brackets indicate the ensemble average over the structures generated in a simulation of R. H A and H R are the Hamiltonians for the real compound (A) and the reference compound (R), respectively. Because this expression only involves the difference between the two Hamiltonians, only interactions that differ between compounds A and R need to be reevaluated over the ensemble. k B is the Boltzmann constant, and T is the temperature. As was demonstrated before, one can split up the process of changing the reference compound into a real ligand into insertion of a real ligand and the removal of the reference compound (6). Assuming that the removal of the reference compound is independent of the real ligand that was added, one can estimate the free-energy difference aswith G...