For further reducing these error, we plan to combine the scheme presented here with the density-based many-body expansion , developed in our group for molecular clusters . Moreover, the efficiency of the calculations can likely be further improved by more refined strategies for screening the two-and three-body contributions. This will open the door to the application of fragmentation schemes in combination with highlyaccurate wavefunction-based quantum chemical methods to proteins and proteinligand complexes.