“…The obtained file was employed for docking simulations performed by Grid-based ligand docking with energetics (GLIDE) [34] , [35] , [36] , [37] , [38] . In order to properly take into account putative conformational rearrangements of the protein binding site during molecular recognition, Induced Fit Docking (IFD) [39] , [40] simulations were performed using the SP mode (all the used parameters are available in Table S1 - Supporting information ) and all the default settings, building a cubic grid centered on the center of mass of H52 (predicted to be protonated during the protein preparation step) and H65 (predicted to be neutral and with the δ -N atom bound to an hydrogen atom), namely the residues mostly affected by the interaction with BPH03 as indicated by NMR data [29] . Furthermore, an inner box of 10 Å × 10 Å × 10 Å and an outer box of 30 Å × 30 Å × 30 Å were employed.…”