The choice of truncated basis sets and their optimization for MBPT calculations of molecular properties are discussed. It is pointed out that computing the correlation corrections to some kth order property by using the MBPT approach requires the knowledge of accurate perturbed orbitals through the kth order. Hence, it is argued that the basis set functions can be optimized with respect to the perturbed energies calculated within the coupled Hartree-Fock method. The proposed procedure is illustrated by MBPT calculations of quadrupole moments of H2 and FH. Additionally, also some estimates of the quadrupole polarizability tensor components for these molecules are obtained.