We have implemented and tested the method we have recently proposed [J. Phys. Chem. Lett. 10, 1537] to treat dispersion interactions derived from a supramolecular wavefunction constrained to leave the diagonal of the many-body reduced density matrix of each monomer unchanged. The corresponding variational optimization leads to expressions for the dispersion coefficients in terms of the ground-state pair densities of the isolated monomers only. We have used three levels of theory for the groundstate monomer pair densities: Hartree-Fock, MP2 and CCSD, looking at both isotropic and anisotropic dispersion coeffcients. For closed-shell systems, CCSD monomer pair densities yield the best results with smaller variance, with mean average percent error on isotropic dispersion coefficents of about 7%, and similar accuracy for anisotropies.The performance for open shell systems is less satisfactory, with CCSD not always providing the best result.