The direct quantum-chemical treatment of large biomolecules such as proteins with accurate methods is generally not feasible because of the steep scaling of the computational effort with system size. However, there is a need for such calculations for numerous applications, ranging from the prediction of protein-ligand interaction energies in drug discovery over studies of reaction mechanisms is biocatalysis to computational enzyme design .