Actinide +IV complexes with six nitrates [AnIV(NO3)6]2- (An = Th, U, Np, and Pu) have been studied by 15N and 17O NMR spectroscopy in solution and first-principles calculations. Magnetic susceptibilities were evaluated experimentally using the Evans method and are in good agreement with the ab initio values. The evolution in the series of the crystal field parameters deduced from ab initio calculations is discussed. The NMR paramagnetic shifts are analyzed based on ab initio calculations. Because the cubic symmetry of the complex quenches the dipolar contribution, they are only of Fermi contact origin. They are evaluated from first-principles based on a complete active space/density functional theory (DFT) strategy, in good accordance with the experimental one. The ligand hyperfine coupling constants are deduced from paramagnetic shifts and calculated using unrestricted DFT. The latter are decomposed in terms of the contribution of molecular orbitals. It highlights two pathways for the delocalization of the spin density from the metallic open-shell 5f orbitals to the NMR active nuclei, either through the valence 5f hybridized with 6d to the valence 2p molecular orbitals of the ligands, or by spin polarization of the metallic 6p orbitals which interact with the 2s-based molecular orbitals of the ligands.