Molecular dynamics simulation of [Gd(egta)(H(2)O)](-) in aqueous solution: internal motions of the poly(amino carboxylate) and water ligands, and rotational correlation times

Molecular dynamics simulations of Gd(egta)(H(2)O)](-) (egta(4-)=3,12-bis(carboxymethyl)-6,9-dioxa-3,12-diazatetradecanedioate(4-)) have been performed without any artificial constraint on the first coordination sphere, such as covalent bonds between the Gd(3+) and the coordination sites. Two new crystallographic structures were determined for this gadolinium chelate and used to start two molecular dynamics simulations. Gd(egta)(H(2)O)](-) and Gd(egta)](-) were both observed during the simulations, with a mean volume for the reaction of dissociation Gd(egta)(H(2)O)](-)-->Gd(egta)](-)+H(2)O of +7.2 cm(3)mol(-1), which corroborates the previously published experimental value of +10.5 cm(3)mol(-1). Changes in the conformation of the complex with the inversion of several dihedral angles are observed in the simulations independently from the water dissociation. Very fast changes of the third-order rotation axis direction of the Gd(3+) coordination polyhedron (of symmetry D(3h)) are observed during the simulations and are related to the mechanism of electronic relaxation of the complex. Different rotational correlation times (tau(R)) were calculated from the simulations on various observables of the complex. Protons of the inner sphere have different tau(R). The mean tau(R) of the two Gd-HW(HW=hydrogen of water molecule) vectors is 72% lower than tau(R) of the complex, and 75% lower than tau(R) of the vector Gd-OW (OW=oxygen of water molecule). This discrimination of the tumbling rates should be taken into account in future global (17)O NMR, EPR and NMRD (nuclear magnetic relaxation dispersion) data analysis.