Density functional theory was applied to study the interaction of group llb transition-metal cations (Zn2+, Cd2+, and Hg2+) with one and two fully or partially deprotonated 3-mercaptopropionic acid ligands. In this investigation, we determined the geometries of all possible complexes resulting from the coordination of the metal ions with the ligands at different binding sites selected on each ligand. The relative energies of the complexes, metal-ion affinities, free energies, and entropies were also determined. The natures of the bonds were critically analyzed by natural bond orbital (NBO) analysis and clarified further using the atoms-in-molecules (MM) approach. The substantial influence of the solvent (water) polarization on the energetics, geometries, and bonding of the molecular complexes was also investigated by the conductor-like screening solvation model (COSMO). In an attempt to simulate the complexes in an aqueous environment, water molecules were added explicitly to complete the coordination sphere of the metal cations, and the corresponding metal-ion affinities were calculated to study the effect of microhydration.