We present a derivation and computations of the paramagnetic enhancement of the nuclear magnetic resonance (NMR) spin-spin coupling, which may be expressed in terms of the hyperfine coupling (HFC) and (for systems with multiple unpaired electrons) zero-field splitting (ZFS) tensors. This enhancement is formally analogous to the hyperfine contributions to the NMR shielding tensor as formulated by Kurland and McGarvey. The significance of the spinspin coupling enhancement is demonstrated by using a combination of density-functional theory and correlated ab initio calculations, to determine the HFC and ZFS tensors, respectively, for two paramagnetic 3d metallocenes, a Cr(II)(acac) 2 complex, a Co(II) pyrazolylborate complex, and a lanthanide system, Gd-DOTA. Particular attention is paid to relativistic effects in HFC tensors, which are calculated using two methods: a nonrelativistic method supplemented by perturbational spin-orbit coupling corrections, and a fully relativistic, 4-component matrixDirac-Kohn-Sham approach. The paramagnetic enhancement lacks a direct dependence on the