The relative binding energies of lanthanide–EDTA complexes, at varying solution pH, were predicted from electronic structure calculations using the solution structures of Ln–EDTA complexes resolved withab initiomolecular dynamics simulations.