In the past decade, there has been a growth in using Zirconium-89 (89Zr) as a radionuclide in nuclear medicine for cancer diagnostic imaging and drug discovery processes. Although one of the most popular chelators for 89Zr, desferrioxamine (DFO) is typically presented as a hexadentate ligand, our work suggests a different scenario. The coordination structure of the Zr4+–DFO complex has primarily been informed by DFT-based calculations, which typically ignore temperature and therefore entropic and dynamic solvent effects. In this work, free energy calculations using molecular dynamics simulations, where the conformational fluctuations of both the ligand and the solvent are explicitly included, are used to compare the binding of Zr4+ cations with two different chelators, DFO and 4HMS, the latter of which is an octadentate ligand that has been recently proposed as a better chelator due to the presence of four hydroxymate groups. We find that thermally induced disorder leads to an open hexadentate chelate structure of the Zr4+–DFO complex, leaving the Zr4+ metal exposed to the solvent. A stable coordination of Zr4+ with 4HMS, however, is formed by involving both hydroxamate groups and water molecules in a more closely packed structure.
