It is shown that the paramagnetic properties of lanthanides can be exploited to obtain information on specific parts of a protein surface. Owing to the high affinity of coordinatively unsaturated lanthanide complexes for oxygen donors, carboxylate groups can be used as preferential targets for the interaction. The DO3A ligand is particularly useful in these studies, as it coordinates lanthanides in a heptadentate fashion, leaving two sites available for exogenous donors. A solution of a (15)N-labeled sample protein, calbindin D(9k) (75 residues), was titrated with up to 200% of Gd(III)-DO3A complex, and an inversion recovery (15)N-(1)H HSQC experiment was used to measure the paramagnetic contributions to the longitudinal relaxation rates of the amide protons. Relaxation data were used as distance constraints to estimate the number of interacting complexes and the occupancies of their binding sites. Four preferential interaction sites on the protein surface are found. Inspection of the various carboxylate side chains on the surface of the protein indicates that Gd(III)-DO3A interacts preferentially with carboxylate-rich regions, rather than with isolated carboxylates, suggesting the possibility of chelation of one Gd(III)-DO3A molecule by two carboxylate groups. Gd(III)-DO3A is thus a valuable semi-selective probe for clusters of negative charges on the protein surface.