We prove the asymptotic normality of nonparametric estimator of pairwise interaction function for a stationary pairwise interaction point process characterized by the Papangelou conditional intensity and observed in a bounded window of a sequence of cubes growing up to $\mathbb{R}^d$. Formula for the variance of the resulting estimator can be obtained using Papangelou conditional intensity of the point process. This is a random function satisfying the counterpart of the Georgii-Nguyen-Zessin formula. The proof of the asymptotic normality of the resulting estimator is based on the $m_n$-approximation method in the setting of dependent random fields indexed by $\mathbb{Z}^d$ where $d$ is a positive integer.