To circumvent inaccuracies of approximate NAKEs for strongly interacting subsystems during molecular dynamics simulations in the liquid phase, Pavanello and co-workers presented a workflow in which subsystem boundaries are dynamically adjusted between different simulation steps , thus making sure that subsystems which interact too strongly are combined into one fragment. As a simple criterion for the interaction strength, the distances between subsystems are employed in that work.