We address the behavior of the Dirac equation with the Killingbeck radial potential including the external magnetic and Aharonov-Bohm (AB) flux fields. The spin and pseudo-spin symmetries are considered. The correct bound state spectra and their corresponding wave functions are obtained. We seek such a solution using the biconfluent Heun's differential equation method. Further, we give some of our results at the end of this study. Our final results can be reduced to their non-relativistic forms by simply using some appropriate transformations. The spectra, in the spin and pseudo-spin symmetries, are very similar with a slight difference in energy spacing between different states.