Molecular determinants and bottlenecks in the unbinding dynamics of biotin-streptavidin

Biotin-streptavidin is a very popular system used to gain insight into protein-ligand interactions. In its tetrameric form, it is well-known for its extremely long residence times, being one of the strongest known non-covalent interactions in nature, and is heavily used across the biotechnological industry. In this work we gain understanding into the molecular determinants and bottlenecks in the unbinding of the dimeric biotinstreptavidin system in its wild type and with N23A mutation. Using new enhanced sampling methods with full atomistic resolution, we reproduce the variation caused by N23A mutation in experimentally reported residence time. We also answer a longstanding question regarding cause/effect in the coupled events of bond stretching and bond hydration during unbinding and establish that in this system, it is the bond stretching and not hydration which forms the bottleneck in the early parts of the unbinding. We believe these calculations represent a step forward in the use of atomistic simulations to study pharmacodynamics. An improved understanding of biotin-streptavidin unbinding dynamics should also have direct benefits in biotechnological and nanobiotechnological applications.