Ligand-protein binding kinetic rates are growing in importance as parameters to consider in drug discovery and lead optimization. In this study we analysed using surface plasmon resonance (SPR) the transition state (TS) properties of a set of six adenosine A2Areceptor inhibitors, belonging to both the xanthine and the triazolo-triazine scaffolds. SPR highlighted interesting differences among the ligands in the enthalpic and entropic components of the TS energy barriers for the binding and unbinding events. To better understand at a molecular level these differences, we developed suMetaD, a novel molecular dynamics (MD)-based approach combining supervised MD and metadynamics. This method allows simulation of the ligand unbinding and binding events. It also provides the system conformation corresponding to the highest energy barrier the ligand is required to overcome to reach the final state. For the six ligands evaluated in this study their TS thermodynamic properties were linked in particular to the role of water molecules in solvating/desolvating the pocket and the small molecules. suMetaD identified kinetic bottleneck conformations near the bound state position or in the vestibule area. In the first case the barrier is mainly enthalpic, requiring the breaking of strong interactions with the protein. In the vestibule TS location the kinetic bottleneck is instead mainly of entropic nature, linked to the solvent behaviour.
