Stability of Dirichlet heat kernel estimates for non-local operators under Feynman-Kac perturbation