The stability of DOTA-chelated radiopharmaceuticals within $^{225}$Ac decay pathway studied with density functional theory.

Ligand structures encapsulating metal ions play a central role as contrast agents in Magnetic Resonance Imaging (MRI) or as agents delivering toxic cargo directly to tumor cells in targeted cancer therapy. The structural stability and interaction with solutions of such complexes are the key elements in understanding the foundation of delivery process. We present a comparative study for the 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid (DOTA) chelated to radioactive isotopes of $^{225}$Ac, $^{221}$Fr, $^{217}$At, $^{213}$Bi and a control $^{68}$Gd. Using density functional theory methods we investigate the structural stability of complexes for cancer therapy including binding energies, charge transfer, electron densities. The van der Waals interactions are included in the simulations to take into account weak dispersion forces present in such structures. Our results reveal that Ac-DOTA, Bi-DOTA and Gd-DOTA are the most stable complexes in the group. We also show that the water environment is a key ingredient for the structural coordination of the DOTA structures.