Using Molecular Dynamics to Predict Factors Affecting Binding Strength and Magnetic Relaxivity of MRI Contrast Agents

Abstract

We demonstrate the use of molecular dynamics and molecular mechanics methods to calculate properties and behavior of metal-chelate complexes that can be used as MRI contrast agents. Static and dynamic properties of several known agents were calculated and compared with experiment. We calculated the static properties such as the q-values (number of inner shell waters) and binding distances of chelate atoms to the metal ion for a set of chelates with known X-ray structure. The dynamic flexibility of the chelate arms was also calculated. These computations were extended to a series of exploratory chelate structures in order to estimate their potential as MRI contrast agents. We have also calculated for the first time the NMR relaxivity of an MRI contrast agent using a long (5 nsec) molecular dynamics simulation. Our predictions are promising enough that the method should prove useful for evaluating novel candidate compounds before they are synthesized. One novel static property, the projected area of chelate atoms onto a virtual surface centered on the metal ion (gnomonic projection), was found to give an effective measure of how well the chelate atoms use the free space around the metal ion.