While the above test cases only consider a single point on the potential energy surface (i.e., a fixed protein structure), molecular dynamics simulations require a balanced description of the low-energy regions of the potential energy surface. This is particularly relevant if quantum-chemical fragmentation methods are to be used for the parametrization of machine learning potentials with the aim to replace classical biomolecular force fields.