The electron paramagnetic resonance (EPR) investigation of mononuclear cis- and trans-(L1O)MoOCl complexes [L1OH = bis(3,5-dimethylpyrazolyl)-3-tert-butyl-2-hydroxy-5-methylphenyl)methane] reveals a significant difference in their spin Hamiltonian parameters which reflect different equatorial and axial ligand fields created by the heteroscorpionate donor atoms. Density functional theory (DFT) was used to calculate the values of principal components and relative orientations of the g and A tensors, and the molecular framework in four pairs of isomeric mononuclear oxo‑molybdenum(V) complexes (cis- and trans-(L1O)MoOCl, cis,cis- and cis,trans-(L-NS)MoOCl [L-NSH = N,N'-dimethyl-N,N'-bis(mercaptophenyl)ethylenediamine], cis,cis- and cis,trans-(L-NS)MoO(SCN), and cis- and trans-[(dt)MoO(OMe)] [dtH = 2,3-dimercapto-2-butene]). Scalar relativistic DFT calculations were conducted using three different exchange-correlation functionals. It was found that the use of hybrid exchange-correlation functional with 25% of the Hartree-Fock exchange leads to the best quantitative agreement between theory and experiment. A simplified ligand-field approach was used to analyze the influence of the ligand fields in all cis- and trans-isomers on energies and contributions of molybdenum d-orbital manifold to g and A tensors and relative orientations. Specifically, contributions that originated from the spin-orbit coupling of the d, d, and d orbitals into the ground state have been discussed. The new findings are discussed in the context of the experimental data of mononuclear molybdoenzyme, DMSO reductase.