On the variational perturbation approach to the Dirac equation