We simulate with state-of-the-art enhanced sampling techniques the binding of
Benzamidine to Trypsin which is a much studied and paradigmatic ligand-protein
system. We use machine learning methods and in particular Time-lagged
Independent Component Analysis to determine efficient collective coordinates.
These coordinates are used to perform On-the-fly Probability Enhanced Sampling
simulations, which we adapt to calculate also the ligand residence time. Our
results, both static and dynamic, are in good agreement with experiments. We
underline the role of water in the unbinding process and find that the presence
of a water molecule located at the bottom of the binding pocket allows via a
network of hydrogen bonds the ligand to be released into the solution. On a
finer scale, even when unbinding is allowed, another water molecule further
modulates the exit time