Approximate analytical solutions of the Dirac equation with Tietz-Hua (TH) potential are obtained for arbitrary spin-orbit quantum number using the Pekeris approximation scheme to deal with the spin-orbit coupling terms 2 ( 1) . r In the presence of exact spin and pseudo-spin (pspin) symmetric limitation, the bound state energy eigenvalues and associated two-component wave functions of the Dirac particle moving in the field of attractive and repulsive TH potential are obtained using the parametric generalization of the Nikiforov-Uvarov (NU) method. The cases of the Morse potential, the generalized Morse potential and non-relativistic limits are studied.The spin and the pseudo-spin (pspin) symmetries of the Dirac Hamiltonian had been discovered many years ago, however, these symmetries have recently been recognized empirically in nuclear and hadronic spectroscopes [1]. Within the framework of Dirac equation, pspin symmetry used to feature deformed nuclei and superdeformation to establish an effective shell-model [2-4], whereas spin symmetry is relevant for mesons [5]. The spin symmetry occurs when the difference of the scalar Sr and vector Vrpotentials are a constant, i.e., s rC and Over the past years, the Nikiforov-Uvarov (NU) method [18] has shown to be a powerful tool in solving second-order differential equations. It was applied successfully to a large number of potential models [19-25]. This method has also been used to solve the spinless (spin-0 ) Schrödinger [26-30] and Klein-Gordon (KG) [31-35] equations and also relativistic spin-2 / 1 Dirac equation [36-40] with different potential models.