This paper addresses probabilistic shaping (PS) which has been a latest key technique to approach capacity of fiber-optic channels. We investigate the impact of PS on nonlinear interference (NLI), including self channel interference (SCI), cross channel interference (XCI), and multi channel interference (MCI) for a polarization multiplexed 16-ary quadrature amplitude modulation format in a wavelength division multiplexed (WDM) system. To this end, we consider performing PS in two scenarios: (i) Solely on the channel of interest and (ii) over all C-band WDM channels of a fiber-optic link by analyzing the effective signal to noise ratio and symbol error rate. It is demonstrated that using the enhanced Gaussian noise model with merely 10% overhead in the first scenario, the applied PS scheme increases the SCI and the total experienced NLI by about 19.23%, and 6.6%, respectively. Interestingly, despite enhancing the NLI in this scenario, the simulated PS technique leads to about 47.6% increase in the transmission reach. In the second scenario, the numerical results show increase of the SCI, XCI, and total NLI around 19.8%, 23.34%, and 20.2%, respectively, but resulting in an increase of 32.3% in the transmission reach.