Basis set superposition error (BSSE) correction with counterpoise (CP) procedure under the environmental electrostatic potential is newly introduced to interfragment interaction energy (IFIE), which is important for interaction analysis in the fragment molecular orbital method. The CP correction for IFIE is applied to a stacked dimer of base pair and a protein-ligand complex of estrogen receptor and 17 beta-estradiol with scaled third-order Moller-Plesset perturbation theory. The BSSEs amount to about quarter of IFIE for hydrogen-bonding and electrostatic interactions and half or even more for dispersion interactions. Estimation of IFIE with the CP correction is therefore preferred for the quantitative discussion. (C) 2011 Elsevier B.V. All rights reserved.