A systematic theoretical investigation on the interaction energies of halogen-ionic bridges formed between halide ions and the polar H atoms bonded to N of protein moieties has been carried out by employing a variety of density functional methods. In this procedure, full geometry optimizations are performed at the Møller-Plesset second-order perturbation (MP2) level of theory in conjunction with the Dunning's augmented correlation-consistent basis set, aug-cc-pVDZ. Subsequently, two distinct basis sets, i.e. 6-311++G(df,pd) and aug-cc-pVTZ, are employed in the following single-point calculations so as to check the stability of the results obtained at the different levels of DFT. The performance of DFT methods has been evaluated by comparing the results with those obtained from the rigorous MP2 theory. It is shown that the B98, B97-1, and M05 give the lowest root-mean-square error (RMSE) for predicting fluoride-binding energies, M05-2X, MPW1B95, and MPW1PW91 have the best performance in reproducing chloride-binding energies, B97-1, PBEKCIS, and PBE1KCIS present the optimal result for bromide-binding energies, while B97-1, MPW1PW91, and TPSS perform most well on iodide-binding energies. The popular B3LYP functional seems to be quite modest for studying halide-protein moiety interactions. In addition, the PBE1KCIS functional provide accuracies close to the computationally expensive MP2 method for the calculation of interaction energies of all halide-binding systems.