Force Fields, Quantum-Mechanical- and Molecular-Dynamics-Based Descriptors of Radiometal-Chelator Complexes

Radiopharmaceuticals are currently a key tool in cancer diagnosis and therapy. Metal-based radiopharmaceuticals are characterized by a radiometal-chelator moiety linked to a bio-vector that binds the biological target (e.g., a protein overexpressed in a particular tumor). The right match between radiometal and chelator influences the stability of the complex and the drug's efficacy. Therefore, the coupling of the radioactive element to the correct chelator requires consideration of several features of the radiometal, such as its oxidation state, ionic radius, and coordination geometry. In this work, we systematically investigated about 120 radiometal-chelator complexes taken from the Cambridge Structural Database. We considered 25 radiometals and about 30 chelators, featuring both cyclic and acyclic geometries. We used quantum mechanics methods at the density functional theoretical level to generate the general AMBER force field parameters and to perform 1 µs-long all-atom molecular dynamics simulations in explicit water solution. From these calculations, we extracted several key molecular descriptors accounting for both electronic- and dynamical-based properties. The whole workflow was carefully validated, and selected test-cases were investigated in detail. Molecular descriptors and force field parameters for the complexes considered in this study are made freely available, thus enabling their use in predictive models, molecular modelling, and molecular dynamics investigations of the interaction of compounds with macromolecular targets. Our work provides new insights in understanding the properties of radiometal-chelator complexes, with a direct impact for rational drug design of this important class of drugs.