We propose a novel Foldy-Wouthuysen transformation including a vector potential A, which can be used to introduce restricted magnetic balance in the Douglas-Kroll-Hess (DKH) method. We also demonstrate that the DKH method can be used in combination with the gauge-including atomic orbital method. For the numerical examination, we calculate the NMR shielding constants and anisotropies of noble gas atoms (Ne, Ar, Kr, Xe), halogens (F, Cl, Br, I), and chalcogens (O, S, Se, Te) in molecules without using extremely large basis sets. (C) 2013 Elsevier B.V. All rights reserved.