Interactions of COOCH3 functional group of methyl benzoate (MB) and its derivatives can substantially be affected by substituents and heteroatoms. The substituent and heteroatom effects on the intermolecular interactions in the complexes of MB and its derivatives with the water molecule have been investigated by the quantum mechanical calculations in the present study. Electron donating (ED)/electron withdrawing (EW) substituents increase/decrease the binding energy (ΔEbin) between mentioned functional group and hydrogen bond donor, and the N heteroatom in the six membered ring weakens the effect. The effects of the N heteroatoms depend on their positions in the ring, and the effects of substituents depend on the number and the positions of heteroatoms. Some relations were obtained between the dipole moment of carbonyl group, results obtained from the natural bond orbital (NBO) and the atoms in molecules (AIM) analysis, binding energy and Hammett constants of substituents. The results can be used to predict the factors affect the intermolecular interactions, and to improve the ligand-protein bindings in biological systems.