Spin-orbit coupling (SOC) is crucial for correct electronic structure analysis in molecules and materials, for example, in large molecular systems such as superatoms, for understanding the role of transition metals in enzymes, and when investigating the energy transfer processes in metal-organic frameworks. We extend the GFN-xTB method, popular to treat extended systems, by including SOC into the hamiltonian operator. We followed the same approach as previously reported for the density-functional tight-binding method and provide and validate the necessary parameters for all elements throughout the Periodic Table. The parameters have been obtained consistently from atomic SOC calculations using the density-functional theory. We tested them for reference structures where SOC is decisive, as in the transition metal containing heme moiety, chromophores in metal-organic frameworks, and in superatoms. Our parameterization paves the path for incorporation of SOC in the GFN-xTB based electronic structure calculations of computationally expensive molecular systems.