Hyperfine interactions (HFI) on the nuclei of the first coordination sphere water molecules in a model [Gd(H2O)(8)](3+) aqua complex and in the magnetic resonance imaging contrast agent [Gd(DOTA)(H2O)](-) were studied theoretically. Density functional theory (DFT) calculations combined with classical molecular dynamics (MD) simulations have been used in order to take into account dynamic effects in aqueous solution. DFT relativistic calculations show a strong spin-polarization of the first coordination sphere water molecules. This spin-polarization leads to a positive O-17 isotropic hyperfine coupling constant (A(iso)(O-17) = 0.58 +/- 0.11 MHz) and to a significant increase of the effective distance (< r(eff)(Gd-O)> = 2.72 +/- 0.06 angstrom) of dipolar interaction compared to the mean internuclear distance (< r(Gd-O)> = 2.56 +/- 0.06 angstrom) obtained from the MD trajectory of [Gd(DOTA)(H2O)](-) in aqueous solution. The point-dipole model for anisotropic hyperfine interaction overestimates therefore the longitudinal relaxation rate of the O-17 nucleus by similar to 45%. The H-1 isotropic hyperfine coupling constant of the bound water molecule is predicted to be very small (A(iso)(H-1) = 0.03 +/- 0.02 MHz), and the point-dipole approximation for first coordination sphere water protons holds. The calculated hyperfine parameters are in good agreement with available experimental data.