Molecular dynamics simulations make it possible to study the structural dynamics of biomolecules in solution, and provide valuable insights into the emergence and mechanisms of biological function . Conventionally, such simulations rely on classical force fields to model the potential energy surface of biomolecules such as proteins . In the past decade, considerable efforts were undertaken to extend the time scales of classical molecular dynamics simulations, which made it possible to directly simulate protein folding processes .