In this paper, we use Krein regularization to study certain standard computations in quantum chromodynamics (QCD). In this method, the auxiliary modes -those with negative norms -are employed to calculate the quark self-energy, vacuum polarizations and vertex functions. We explicitly show that after making use of these modes and by taking into account the quantum metric fluctuation for the problems at hand, the conventional results can indeed be reproduced; but with the advantage of finite answers which require fewer mathematical procedures. An obvious merit of this approach is that the theory is naturally renormalized. The ultraviolet (UV) divergences disappear due to the presence of negative norm state, similar to the Pauli-Villars regularization method. We compare the answers of Krein regularization with the results of calculations which have been done in Hilbert space.