Solid-state nuclear magnetic resonance (NMR) of paramagnetic samples has the potential to provide a detailed insight into the environments and processes occurring in a wide range of technologically-relevant phases, but the acquisition and interpretation of spectra is typically not straightforward. Structural complexity and/or the occurrence of charge or orbital ordering further compound such difficulties. In response to such challenges, the present article outlines how the total Fermi contact (FC) shifts of NMR observed centers (OCs) may be decomposed into sets of pairwise metal–OC bond pathway contributions via solid-state hybrid density functional theory calculations. A generally applicable “spin flipping” approach is outlined wherein bond pathway contributions are obtained by the reversal of spin moments at selected metal sites. The applications of such pathway contributions in interpreting the NMR spectra of structurally and electronically complex phases are demonstrated in a range of paramagnetic Li-...