Magnetic shielding and indirect spin-spin coupling phenomena are tensorial properties and both their isotropic and anisotropic parts do affect NMR spectra. The involved interaction tensors, (sigma) over cap and (J) over cap, can nowadays be theoretically calculated, although the reliability of such methods in the case of anisotropic parameters, Delta sigma and Delta J, in systems involving heavy nuclei, yet demands testing. In this communication the results of the experimental and theoretical investigations of bis(phenylethynyl)mercury (I) labeled with C-13 isotope at positions neighboring Hg are reported. The theoretical calculations of molecular geometry and values of NMR parameters for I have been performed by the ZORA/DFT method, including the relativistic scalar and spin-orbit coupling contributions, using the PBE0 functional and TZP (or jcpl) basis set. These values have been confronted with the experimentally measured ones. The isotropic parameters have been measured by the standard C-13 and Hg-199 NMR spectra. The shielding anisotropies for the atoms in the central part of molecule I have been determined in a liquid sample using magnetic relaxation measurements. The relaxation data have been interpreted within the rotational diffusion theory, assuming the symmetrical top reorientation model. The anisotropies of one-bond C-13-Hg-199 and two-bond C-13-Hg-C-13 spin-spin couplings have been determined exploiting the temperature-dependent C-13 NMR spectra of I in the ZLI1167 liquid-crystal phase. We have found that our theoretical calculations reproduce experimental values of both isotropic and anisotropic NMR parameters very well.