The role of the Hulthén potential on the spin and pseudospin symmetry solutions is investigated systematically by solving the Dirac equation with attractive scalar
and repulsive vector
potentials. The spin and pseudospin symmetry along with orbital dependency (pseudospin–orbit and spin–orbit dependent couplings) of the Dirac equation are included to the solution by introducing the Hulthén‐square approximation. This effective approach is based on forming the spin and pseudo‐centrifugal kinetic energy term from the square of the Hulthén potential. The analytical solutions of the Dirac equation for the Hulthén potential with the spin–orbit and pseudospin–orbit‐dependent couplings are obtained by using the Nikiforov–Uvarov (NU) method. The energy eigenvalue equations and wave functions for various degenerate states are presented for several spin–orbital, pseudospin–orbital and radial quantum numbers under the condition of the spin and pseudospin symmetry.