Modified perturbation-variation method for calculating spin-spin coupling constants